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I.  INTRODUCTION 


In  the  field  of  quasi-one-dimensional  conductors  much  interest 
has  been  generated  by  the  recent  discovery! >2  of  anomalous  electronic 
transport  properties  in  the  linear-chain  metal  NbSe3.  This  interest  is 
two-fold. 

First,  from  the  theoretical  point  of  view  the  observed  properties 
give  every  indication  that  a  new  many-body  effect  involving  coherent 
electronic  charge  transport  is  being  excited  in  NbSe3  below  the  phase  tran¬ 
sition  where  charge-density -waves  (CDW)  form.  The  theory  that  the  coupled 
electron-and-lattice  condensate  should  display  conduction  processes  such  as 
sliding  conductivity  and  phase  oscillations  has  been  discussed  by  many  workers3 
following  the  original  work  of  FrohlicM*.  Attempts  to  apply  these  ideas  to 
real  solids  have  been  controversial  in  the  recent  past7  because  the  experi¬ 
mental  evidence  is  often  marginal  or  not  reproducible.  In  most  cases 
conventional  single-particle  theory  can  accommodate  the  observed  transport 
properties.  In  NbSe3  many-body  effects  on  the  electronic  transport  are  so 
predominant  that  they  were  inferred  to  be  the  signature  of  the  Frohlich 
mechanism  long  before  it  was  experimentally  established  that  the  phase 
transitions  were,  in  fact,  caused  by  CDW  formation.  These  anomalous 
transport  properties  are  easily  observed,  reproducible  from  sample 
to  sample,  and  indisputably  at  variance  with  (by  orders  of  magnitude) 
the  predictions  of  single-particle  theory.  The  dynamics  of  CDW  metals 
provide  a  third  situation  where  many-body  effects  dominate  the  electronic 
transport.  (The  other  two  are  superconductivity  and  the  Kondo  effect.) 
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The  second  viewpoint  is  from  the  applications  aspect.  Much  of 
the  impetus  behind  activity  in  this  field  has  come  from  the  hope  that 
the  Frohlich  mechanism^  may  lead  to  high  temperature  superconductivity. 

This  hope  has  dimmed  considerably,  partly  because  none  of  the  one¬ 
dimensional  conductors  studied  exceed  3,000  (ftcm)'^  in  conductivity 
at  room  temperature  or  show  great  conductivity  enhancement  in  the  CDW 
phase.  However,  study  of  the  (non-Ohmic)  mechanism  in  NbSej  (which  is 
by  far  the  strongest  candidate  for  the  Frohlich  mechanism)  will  increase 
our  understanding  of  the  Frohlich  mechanism  (especially  on  the  question 
of  its  pinning  to  the  host  lattice  or  suppression  in  most  CDW  systems.) 

Quite  apart  from  such  considerations  the  non-linear  I-V  characteristics 
in  the  CDW  phase  may  be  exploited  for  device  applications.  The 
strong  dispersion  of  the  Ohmic  conductivity  in  NbSe3  in  the  GHz  region 
together  with  the  I-V  non-linearities  may  also  find  useful  applica¬ 
tions  in  signal  detection  at  these  frequencies.  (Doping  studies 
show  that  the  dispersion  can  be  pushed  into  the  100GHz  region.) 


U.  PHYSICS  OF  PHASE  TRAMSITIONS 
AND  ANOMALOUS  TRANSPORT  IN  NbSe3 

As  a  background  for  the  following  report  in  III  on  the  work  per¬ 
formed  on  NbSe3  we  present  here  a  brief  discussion  of  the  physics  under¬ 
lying  the  various  experiments. 

PEIERLS  DISTORTION^:  A  strictly  one-dimensional  electron 
gas  confined  to  a  chain  of  atoms  of  spacing  a  is  inherently  unstable 
to  the  following  distortion.  The  ions  move  to  new  equilibrium  posi¬ 
tions  which  are  described  by  a  static,  longitudinal  sinusoidal  wave. 

This  superlattice  distortion  presents  a  new  periodic  potential  which 
induces  new  gaps  in  the  electronic  band.  If  the  superlattice  periodi¬ 
city  is  such  that  the  new  gaps  occur  at  the  Fermi  level  of  the 
original  band  then  the  electron  gas  gains  energy,  because  the 
opening  of  the  gaps  lowers  the  energy  of  the  single-particle  states 
below  the  gap.  A  phase  transition9  into  the  distorted  phase  occurs 
when  the  gain  in  electronic  energy  offsets  the  cost  in  strain  energy 
coming  from  the  periodic  lattice  distortion.  In  most  systems  the 
transition  occurs  as  a  function  of  temperature  because  a  decrease  in 
temperature  sharpens  the  Fermi -Dirac  distribution  and  enhances  the 
gain  in  electronic  energy.  Conversely,  the  effect  of  pressure  is  to 
increase  the  strain  energy  cost  of  the  lattice  which  favors  the 
normal  phase.  This  leads  to  a  decrease  of  the  transition  temperature 
with  increasing  pressure. 
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In  the  condensed  phase  the  static  lattice  distortion  appears  as 
faint  superlattice  Bragg  spots10  in  x-ray  and  electron  scattering 
experiments.  The  charge  density  of  the  condensed  electrons  is 
also  modulated  by  the  same  periodicity  (hence  the  name  charge-density- 
wave  ( CDW ) ) -  Since  the  periodicity  is  tied  to  the  Fermi  Surface 
dimension,  the  COW  can  either  be  commensurate  with  the  original  host 
lattice  (CCDW)  or  incommensurate  (ICDW).  The  conditions  under  which 
the  CDW  is  commensurate  or  incommensurate,  or  changes  from  one  to 
the  other,  is  at  present  an  active  area  of  research11. 

ELECTRON  TRANSPORT  IN  CDW  STATE.  In  the  foregoing  picture  the 
CDW  transition  converts  a  partially  filled  band  (metal)  to  a  fill ed- 
band  (below  the  CDW  gap).  Thus  in  a  one-dimensional  system  the 
transition  is  from  metal  to  semi-conductor,  as  is  in  fact  observed 
in  KCP,  TTF-TCNQ  and  its  analogs  and  in  TaS3-  However,  in  the  two- 
dimensional  layered-compounds10  (TaSe^,  NbSe^  etc.)  and  in  NbSe3 
the  transition  is  from  metal -to-metal .  Presumably,  other  electronic 
bands  which  do  not  participate  in  the  transition  survive  to  provide 
metallic  conduction  in  the  low  temperature  phase.  In  the  higher¬ 
dimensional  systems  there  exist  flat  portions  of  the  Fermi  surface 
(FS)  which  are  sufficiently  parallel  ("nesting"  criterion)  to  drive 
the  CDW  instability.  The  superlattice  periodicity  is  then  determined 
by  the  wave  vector  ("spanning  vector")  which  --  by  translation  in  k- 
space  —  brings  one  portion  of  the  FS  into  coincidence  with  its  nesting 
partner.  Gaps  occur  over  these  nesting  portions  only.  Other  pockets 
of  the  FS  are  not  affected  by  the  CDW  transition.  To  clarify  the  dis¬ 
cussion  of  the  conduction  mechanisms  in  the  CDW  state  we  shall  distinguish 
the  uncondensed  (free)  carriers  from  the  carriers  in  the  bands  that 


develop  a  gap.  The  latter  will  be  called  collectively  the  condensate. 


Note  that  very  near  the  transition  temperature  free  carriers  can 
be  thermally  generated  from  the  condensate  as  well  if  kT  >  acdw. 

Further,  to  distinguish  the  anomalous  transport  properties  next  to  be 
discussed,  we  shall  call  the  transport  by  the  uncondensed  and  thermally 
generated  carriers  "normal"  conductivity  (even  though  we  are  in  the 
CDW  phase).  Thus,  in  the  normal  phase  (T  >  T^)  current  is  carried 
by  all  portions  of  the  FS.  Below  T^q^  gaps  occur  over  the  nesting 
portions.  The  gaps  grow  from  0  to  their  maximum  value  in  analogy  to  the 
case  in  superconductivity.  Thermal  generation  of  quasi-particles  across 
the  growing  gap  leads  to  a  semiconducting-type  conductivity  contri¬ 
bution  from  the  nesting  portions.  At  very  low  temperatures  when  the 
gaps  are  fully  developed  these  quasi-particles  are  frozen  out 
leaving  only  the  uncondensed  carriers  to  carry  the  current.  The  normal 
conductivity  decreases  to  a  value  appropriate  to  the  reduced  FS  area, 
but  remains  metallic  down  to  the  lowest  temperature. 


MANY-BODY  EFFECTS:  FR0HL1CH  NODE,  DISCOMMENSURATIONS.  If  the  CDW 
is  commensurate  with  the  host  lattice,  the  phase  of  the  sinusoidal 
electronic  charge  density  (relative  to  the  host)  can  be  adjusted  so 
that  peaks  in  the  electronic  density  occur  between  adjacent  metal  ions. 
This  enhances  the  gain  in  bonding  energy.  The  CCDW  is  said  to  be 
"locked-in"  to  the  host.  In  the  incommensurate  case  "lock-in"  is  not 
possible  if  the  ICDW  is  rigid  because  the  gain  in  bonding  energy 
averages  to  zero  over  the  whole  crystal.  This  implies  that  the  phase 
4>  has  no  preferred  value  in  the  ground  state  if  the  host  lattice  is 
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free  of  defects.  In  particular  <<>  can  be  made  to  vary  linearly  in  time 
with  no  cost  in  energy.  This  would  correspond  to  the  condensate  sliding 
relative  to  the  host  lattice,  and  is  the  basis  for  charge  transport 
in  the  Frohlich  model^’®. 

In  a  real  solid  this  sliding  mode  is  impeded  by  pinning  to  im- 

3  13  5 

purities  and  lattice  defects  .  Fukuyama  and  Lee  ,  and  Lee  and  Rice 

have  studied  the  pinning  effects  of  impurities  on  an  ICDW.  They 

distinguish  between  the  two  cases  of  strong  and  weak  pinning.  In 

the  strong  pinning  case  the  impurity  potential  is  sufficiently  strong 

that  <j>  adjusts  itself  locally  at  each  impurity  site.  In  the  weak 

pinning  case  the  pinning  energy  gained  by  local  phase  deformation  at 

each  site  is  not  sufficient  to  overcome  the  cost  in  strain  energy. 

However,  pinning  still  occurs  as  follows.  The  ICDW  breaks  up  into 

finite  domains.  Within  each  domain  the  phase  is  more  or  less  uniform 

and  assumes  an  optimum  value  determined  by  all  the  weak  impurities 

within  the  domain.  The  total  pinning  energy  gained  from  this  collection 

of  impurities  now  becomes  sufficient  to  overcome  the  strain  energy  in 

the  domain  walls  (where  <j>  changes  rapidly).  The  ICDW  is  therefore 

pinned  by  the  collection  of  impurities.  Moreover,  the  pinning  energy 

is  expected  to  increase  with  the  concentration  of  weak  impurities. 

When  an  electric  field  is  applied  the  condensate  does  not  move  until 

the  field  exceeds  some  threshold  value  required  to  dislodge  the  domain. 

Such  a  stick-slip  process  is  expected  to  give  rise  to  strong  non- 
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linearities  in  the  conductivity.  Lee  and  Rice  have  studied  this 
problem  and  predict  that  the  threshold  field  Eg  should  increase  as 

E0  -  c2 


(D 


where  c  is  the  concentration  of  weak  impurities. 

The  impurities  also  affect  the  ac  response  of  the  condensate. 

The  pinned  condensate  will  in  general  have  a  characteristic  pinning 
frequency  so  that  an  applied  weak  ac  field  near  the  pinning  frequency 
will  induce  segments  of  the  condensate  to  undergo  forced  oscillations. 
The  pinning  frequency  may  be  quite  low  if  the  pinning  is  weak.  Again, 
since  the  pinning  energy  increases  with  impurity  concentration  the 
pinning  frequency  should  also  increase  with  c. 

In  summary,  over  and  above  the  normal  conductivity  carried  by 
the  uncondensed  carriers  the  CDW  state  should  display  two  conductivity 
mechanisms  which  are  manifestations  of  the  Frohlich  mode.  These  are 
non-Ohmicity  in  the  dc  conductivity  at  relatively  low  fields,  and 
excess  absorption  at  a  relatively  low  rf  frequency.  Since  these  are 
mechanisms  involving  collective  motion  of  the  charged  condensate  they 
cannot  be  described  by  the  single-particle  models,  and  their  effects 
on  the  transport  may  manifest  themselves  as  spectacular  departures 
from  conventional  Drude  behavior. 
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III.  SUMMARY  OF  EXPERIMENTS  DONE  ON  NbSe3  AT  USC 


The  report  on  the  experiments  performed  on  NbSe^  at  USC  in  the 
last  two  years  (and  funded  by  the  Office  of  Naval  Research)  follows 
more  or  less  the  chronological  order  in  which  the  experiments  were 
done.  This  sequence  of  presentation  will  show  in  a  more  transparent 
way  the  evolution  of  the  physical  ideas  used  to  describe  NbSe^. 

Some  significant  work  performed  at  other  institutions  will  be  men¬ 
tioned  where  appropriate. 

1 )  Non-Ohmic  Conductivity 

Following  the  original  observation  by  Monceau,  Ong,  and  co- 

workers^  at  U.C.  Berkeley  that  NbSe^  is  non-Ohmic  below  the  two 

phase  transitions  at  T^  (142  K)  and  T£  (58  K),  as  well  as  anomalous 
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in  its  microwave  absorption  at  9.8  GHz,  we  set  out  to  measure  more 
systematically  the  temperature  dependence  of  the  non-Ohmic  para¬ 
meters.  Monceau,  et  alJ,  had  found  that  the  non-Ohmicity  obeyed 
the  relation 

-En(T)/E 

o(E,  T)  =  oa(T)  +  ab(T)  e  0  (2) 

The  temperature  dependence  of  ogS  ab  and  E^  for  both  transitions 
2  14 

were  measured  by  us  ’  and  interpreted  in  terms  of  the  CDW  model 
as  follows.  Arguing  from  the  separability  of  the  Ohmic  conduc¬ 
tivity  into  two  parts,  one  coning  from  the  non-nesting  FS  portions 

and  the  other  from  the  condensate  through  thermal  generation  of 

14 

carriers,  we  formed  the  parameter 
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i 


a(T)  s  ab(T)/(aa(T)  +  ob(T))  , 


(3) 


and  related  it  to  the  CDW  cjap  a  by 

o(T)  =  a(0)  tanh  ba/2  >  e  =  (kgT)'1  .  (4) 

Data  on  a  confirmed  that  it  rapidly  saturated  to  a  constant  value 
as  T  decreased,  as  required  by  Eq.  (4).  Moreover,  by  fitting 
Eq.  (4)  to  the  data  we  extracted  values  for  the  gap  at  both  the 
and  transition.  This  was  the  first  measurement  of  the  CDW 
gap  in  NbSe3,  and  lent  strong  support  to  the  unconfirmed  (at  that 
time)  CDW  model. 

2)  Young's  Modulus 

The  CDW  transition  may  also  be  regarded  as  the  softening  and 

condensation  of  the  longitudinal  phonon  with  q  =  2kp.  There  is 

therefore  some  likelihood  that  acoustic  or  mechanical  probes  may 

be  sensitive  to  the  transition.  Working  with  Dr.  Brill  at  Stanford 
1 5 

we  measured  the  Young's  Modulus  in  NbSe^  and  detected  an  elastic 
anomaly  at  the  T^  transition.  However,  no  anomaly  was  seen  to  the 
limits  of  our  uncertainty  (10~^)  at  the  T 2  transition.  The  modulus 
was  measured  using  a  vibrating  reed  technique. 

3)  Hall  Effect 

One  of  the  consequences  of  the  appearance  of  a  gap  over  por¬ 
tions  of  the  FS  is  the  loss  of  carriers.  If  the  CDW  model  is 
valid  in  NbSe3  the  loss  of  carriers  at  both  T-j  and  T2  should  have 
a  drastic  effect  on  the  Hall  constant.  We  succeeded  in  attaching 
transverse  probes  to  the  small  samples  and  were  able  to  measure 


the  Hall  constant  both  as  a  function  of  magnetic  field  and  tem¬ 
perature^.  The  Hall  constant  does  change  significantly  at  both 
T-|  and  T2.  However,  the  relative  change  is  much  larger  than  what 
one  might  have  inferred  (on  the  basis  of  a  single-band  model)  from 
the  conductivity  data.  The  Hall  constant  also  changes  sign  at 
low  temperature  indicating  the  presence  of  both  holes  and  elec¬ 
trons.  Clearly  a  full  two-band  analysis  would  be  required  to 
relate  the  measured  Hall  constant  directly  to  the  carrier  concen¬ 
trations. 

4)  Conductivity  Anisotropy,  Maqnetoresistance  and  Two-Band 
Model 

The  conductivity  anisotropy  and  transverse  magnetoresistance 
were  measured^  to  provide  enough  data  for  a  full  two-band  analysis. 
Together  with  the  Hall  data,  these  measurements  provided  six  ex¬ 
perimental  numbers  to  fix  the  six  parameters  of  the  two-band 
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model  at  each  temperature  below  T^.  Finally  we  were  able  to 
extract  the  number  of  electrons  and  holes  and  their  respective 
mobilities  along  the  chain  and  transverse  directions.  This  rather 
complete  analysis  provides  full  information  on  the  mobilities 
and  concentrations  in  the  CDW  phase  below  T2.  To  our  knowledge 
no  such  analysis  has  been  carried  out  for  any  other  CDW  system. 

5)  Superlattice  Measurement 

At  about  this  time  various  groups  reported  the  successful 
detection  of  superlattice  spots  associated  with  the  T^  and  T2 
transitions.  Tsutsumi  et  alJ®,  (Hokkaido  University)  using 


electron  diffraction  measured  the  spanning  vector  of  the  T-| 
CDW  to  be  (in  units  of  a*,  b*,  c*) 
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Q1  =  (0,  0.25,  0)  .  (5)  j 

I 

1 

I 
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Fleming  and  co-workers  (Bell  Labs)  using  diffuse  x-ray  scattering 
determined  (for  the  l^  CDW)  to  be 

(6) 

20 

In  addition,  Fleming  et  al.  ,  measured  the  superlattice  spot  inten¬ 
sity  in  the  presence  of  a  strong  electric  field  and  found  no 
decrease  in  the  spot  intensity.  In  other  words,  the  non-Ohmicity 
is  not  associated  with  a  break-up  of  the  CDW  condensate. 

6)  Effect  of  Impurities  on  the  Anomalous  Transport  Properties 
We  next  inquired  whether  the  Frohlich  mechanism,  which  can 
successfully  account  for  the  non-Ohmicity  as  well  as  the  excess 
microwave  absorption  mode,  was  capable  of  making  predictions  which 
can  be  tested.  Since  the  Frohlich  mode  is  sensitive  to  the  pre¬ 
sence  of  impurities  the  first  obvious  test  is  to  dope  the  nominally 
pure  compund  with  weak  impurities  and  observe  their  effect  on  both 

the  non-Ohmicity  and  microwave  absorption.  Specifically,  Lee  and 

5  .  2 

Rice  predict  that  the  threshold  field  Eg  should  increase  as  c 

(Eq.  (1)).  The  excess  microwave  absorption  should  also  decrease 

rapidly  as  the  impurities  stiffen  the  effective  spring  constant 

and  drive  up  the  pinning  frequency.  These  predictions  have  been 

verified  by  our  impurity  studies  .  The  nominally  pure  compound 

was  doped  with  Ta  in  the  1000  ppm  range  and  both  the  dc  Ohmic 

breakdown  and  the  excess  microwave  absorption  were  measured  as  a 

function  of  c,  the  Ta  concentration.  Preliminary  measurements 


Q2  =  (0.5,  0.26,  0.5)  . 


showed  that  both  effects  are  dramatically  suppressed  by  increasing 
c.  The  microwave  absorption  mode  at  9.3  GHz  decreases  monotonically 
with  increasing  c,  consistent  with  the  pinning  frequency  increasing 


with  c.  The  effect  on  the  depinning  field  Eg  in  the  dc  non-Ohmic 
measurements  is  even  more  pronounced.  Eg  rises  rapidly  as  c  in¬ 
creases.  Measurements  on  17  samples  showed  that  Eg  increases  as 
2 

c  in  striking  agreement  with  the  predictions  of  Lee  and  Rice. 

In  some  samples  Ti  was  used  as  dopant  instead  of  Ta.  Since 
Ti  is  group  IV  it  enters  the  lattice  as  a  charged  impurity  and 
is  expected  to  behave  as  a  strong  pinning  center.  Our  results 
show  that  its  effectiveness  in  suppressing  both  the  non-Ohmicity 
and  the  excess  microwave  absorption  is  ten  times  larger  than  for 
Ta  in  the  same  concentration.  This  confirms  that  charged  defects 
are  much  more  effective  in  pinning  the  condensate  than  uncharged 
ones. 

7)  Thermo power 

22 

In  collaboration  with  Drs.  Dee  and  Chaikin  at  UCLA  we 
have  measured  the  thermopower  in  the  regime  where  non-Ohmicity 
appears.  This  was  accomplished  by  applying  a  strong  low-frequency 
ac  field  to  the  sample.  In  the  absence  of  a  thermal  gradient 
across  the  length  of  the  sample  no  dc  voltage  develops  across  the 
sample.  When  a  temperature  gradient  &T  is  applied  a  dc  voltage 
proportional  to  aT  appears.  If  the  applied  ac  field  is  zero  this 
dc  voltage  gives  the  familiar  thermopower  S.  In  the  presence  of 
the  strong  ac  field  we  measure  an  enhancement  of  S  which  is  quite 
unexpected.  The  simple  Frohlich  model  would  predict  a  decrease 
in  S  when  the  condensate  is  made  to  slide.  At  present  the  con- 
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tradictory  results  are  not  understood  and  further  work  is  underway. 

8)  Far-Infrared  Magnetoabsorption 

Dr.  R.  Wagner  (NRL,  Wahsington,  D.C.)  in  collaboration  with 
23 

our  group  has  detected  cyclotron  resonance  lines  in  MbSe^.  The 
absorption  lines  occur  in  the  far-infrared  at  fields  of  the  order 
of  8  Tesla.  Sample  inhomogeneities  so  far  have  prevented  a  com¬ 
plete  analysis  of  the  line  shapes  despite  the  good  reproducibility 
of  the  data.  This  is  also  a  continuing  project. 

Present  Status  and  Future  Directions 

Two  and  a  half  years  ago,  the  nature  of  the  phase  transitions 
in  NbSe^  and  the  intriguing  anomalous  transport  properties 
associated  with  the  new  phases  was  unknown.  In  the  intervening 
period  much  work,  mainly  transport  measurements,  has  been  done  to 
study  the  transitions  as  well  as  the  anomalous  properties.  The 
CDW  model  has  been  placed  on  a  firm  footing  with  the  discovery 
of  two  incommensurate  superlattices.  Furthermore,  rather  complete 
measurements  on  the  galvanomagnetic  properties  in  the  Ohmic  (normal) 
regime  have  been  done,  and  analyzed  in  terms  of  a  two-band  model. 

All  the  transport  measurements  in  the  anomalous  (non-Ohmic,  high 
frequency)  regimes  support  the  Frohlich  model  of  a  sliding  CDW. 

This  model  has  been  shown  to  have  impressive  predictive  powers 
in  the  study  of  impurity  effects,  although  in  the  thermopower 
measurements  it  contradicts  experiment.  If  the  Frohlich  model 
turns  out  to  be  the  correct  one,  NbSe^  will  represent  the  only 
system  (so  far)  in  which  this  new  mechanism  of  charge  transport 
(sliding  conductivity)  can  be  conveniently  studied. 

The  questions  to  be  addressed  in  future  programs  are: 


1)  What  limits  the  conductivity  of  the  sliding  condensate? 

This  problem  clearly  requires  a  better  understanding  of  the  pinning 
mechanism  of  various  kinds  of  defects.  The  study  of  T  and  Ti 
impurities  is  a  start  in  this  direction.  A  related  issue  is  the 
viscous  coupling  between  the  uncondensed  carriers  and  the  condensate. 
This  may  ultimately  limit  the  maximum  conductivity. 

2)  What  is  the  galvanomagnetic  response  of  the  sliding  conden¬ 
sate?  Hall  and  magnetoresistance  measurements  are  feasible  in  the 
anomalous  regimes  (non-Ohmic  and  microwave  frequencies).  The 
response  of  the  condensate  current  to  high  magnetic  fields  should 

be  novel  and  interesting. 

3)  What  is  the  transient  response  of  the  condensate?  The 
anomalously  low  pinning  frequency  of  the  condensate  implies  that 
the  (low-current)  conductivity  has  a  measurable  frequency  depen¬ 
dence  in  the  MHz  region.  Such  long  time  constants  in  the  elec¬ 
tronic  response  have  not  been  previously  encountered  in  a  bulk, 
homogeneous  metal. 
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Using  the  Barmatz  vibrating  reed  technique  we  have  measured  the  Young’s 
modulus  of  NbSe3.  An  anomaly  is  observed  at  T,  (142  K.)  but  not  at  7) 

(58  K).  Taken  together  with  published  specific  heat  data  our  results 
provide  upper  (lower)  bounds  for  the  anisotropy  of  the  transition-temper¬ 
ature-stress-dependence  at  T\  (T2). 


1.  INTRODUCTION 

THE  TRANSITION-METAL-TR1CHALCOGENIDE 
NbSe3  has  aroused  much  interest  because  of  its  unusual 
transport  properties  [1,2].  The  resistivity  is  metallic  and 
shows  two  giant  anomalies  [3]  indicative  of  phase 
transitions  occuring  at  7”,  (142  K)  and  7j  (58  K).  At 
temperatures  below  142  K,  the  conductivity  is  highly 
non-Ohmic  and  has  been  shown  to  behave  [2]  as 

o  =  +  0bn  e~e°n/e,  n  =  1,2  (1) 

where  e  is  the  applied  electric  field  and  oan  and  obn  are 
field-independent  quantities.  In  equation  (1),  the  second 
subscript  n  labels  the  phase  transition.  The  activation 
field  e0„  diverges  close  to  the  transition  temperature  and 
the  quantity  a„  =  ob„/(oan  +  obn)  has  been  shown  [4] 
to  behave  like  an  order  parameter. 

NbSe3  undergoes  a  third  transition  to  the  supercon¬ 
ducting  state  [5]  when  hydrostatic  pressure  is  applied 
(Tc  is  6  K  at  10  K  bar)  and  recent  Hall  effect  [6] 
measurements  show  a  striking  loss  of  carriers  at  7)  and 
T2  which  correlates  with  the  rapid  rise  in  the  resistivity 
at  these  temperatures.  The  frequency  dependence  of 
the  resistivity  is  also  anomalous  [2] .  Both  anomalies 
are  sharply  reduced  at  10  GHz.  The  transitions  have 
been  tentatively  associated  [1—3]  with  charge-density- 
wave  [7]  (CDW)  formation;  the  non-Ohmic  behavior 
and  anomalous  frequency  dependence  of  the  conduc¬ 
tivity  have  been  described  [2, 8]  in  terms  of  excitations 
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in  the  CDW  condensate.  However,  this  hypothesis 
remains  unverified  by  direct  X-ray  evidence. 

Since  CDW  formation  involves  the  softening  and 
freezing  of  a  phonon  mode,  a  super-lattice  distortion 
corresponu.og  to  the  condensed  mode  appears  below 
the  CDW  transition.  It  is  natural  to  investigate  the 
elastic  properties  of  the  system  to  detect  anomalies  at 
these  transitions.  In  the  layered  dichalcogenides 
2H-TaSe2  and  2H-NbSe2,  where  CDW  formation  has 
been  experimentally  established  [7] ,  strong  anomalies 
in  the  Young’s  modulus  were  detected  [9, 10]  at  the 
CDW  transitions.  The  internal  friction  also  showed 
striking  changes  near  these  phase  transitions.  Anomalies 
have  also  been  observed  at  the  Peierl’s  distortion 
temperatures  in  the  one-dimensional  conductors 
TTF-TCQN  [11]  and  tri-methyl-ammonium-iodide- 
TCNQ  [12] .  To  acquire  more  information  on  the  tran¬ 
sitions  in  NbSe3,  we  have  performed  Young’s  modulus 
measurements  on  single  crystals,  utilizing  the  Barmatz 
[13]  vibrating  reed  technique. 

The  results  for  the  two  samples  investigated  show  a 
small  anomaly  in  the  modulus  at  142  K,  but  no  evidence 
of  the  58  K  transition  was  detected.  Using  Testardi’s 
theory  [14] ,  we  show  from  our  results  that  the  stress 
dependence  of  T2  is  highly  anisotropic.  The  7)  transition 
is  less  anisotropic,  but  a  numerical  estimate  is  not 
available  until  the  specific  heat  anomaly  at  7)  is 
measured. 

2.  EXPERIMENTAL 


The  flexural  resonant  frequency  of  a  cantilevered 
beam  of  length  /  and  thickness  t  is  given  by  [9,13] 
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Fig.  1.  The  change  in  Young’s  modulus  E  vs  temperature 
for  two  samples.  The  two  curves  for  sample  No.2  near 
room  temperature  indicate  hysteresis  between  cooling 
(solid  line)  and  warming  (dashed  line)  attributed  to 
creep  in  the  damp.  The  anomaly  near  141 K  is  ident¬ 
ified  with  the  large  resistivity  anomaly  observed  at 
r,  =  142  K.  £0(300  K)  ~  3  x  1012  dyne  cm’1 . 

where  tc  =  r/(12)l/2  is  the  radius  of  gyration,  the  k„  are 
constants  1 .875, 4.694, . . . ,  E  is  the  Young’s  modulus 
in  the  axial  direction,  and  p  is  the  density. 

Single  crystals  of  NbSe3  were  mounted  on  the  end 
of  copper  rods  using  low  temperature  solder  [15] .  The 
free  end  of  the  crystal  was  placed  between  two 
electrodes  ~  3  mm  apart.  An  a.c.  voltage  applied  to  one 
of  the  electrodes  excited  flexural  oscillations  of  the 
reed  at  twice  the  applied  frequency  [13] .  Typical 
strains  used  at  resonance  were  e  <  1CT6.  A  d.c.  bias  was 
applied  through  lOGft  to  the  second  electrode,  which 
picked  up  the  signal  due  to  the  modulated  capacitance 
to  the  grounded  reed.  A  MOSFET  mounted  directly  on 
the  pick  up  electrode  limited  the  static  capacitance  to 
ground  to  ~  3  pf.  The  signal  was  detected  by  a 
PAR  129A  lock-in  amplifier  operating  in  the  If  mode. 
The  output  of  the  lock-in  was  used  to  phase  lock  the 
drive  signal  to  the  resonance.  The  measurement  tech¬ 
nique  is  described  in  more  detail  elsewhere  [9, 12] . 

The  lock-in  output,  and  therefore  the  change  in 
resonant  frequency  (and  relative  change  in  modulus) 
and  temperature  are  monitored  continuously  on  an 
XY  recorder. 


Fig.  2.  XY  recorder  trace  for  four  runs  through  the  T, 
transition  in  sample  No.  2;  increasing  modulus  is  toward 
bottom  of  figure.  The  upper  two  curves  were  obtained 
using  a  d.c.  bias  of  50  V  and  driving  voltage  of  6  V.  The 
third  and  fourth  curves  have  a.c.  and  d.c.  voltages 
respectively  doubled.  The  arrows  indicate  warming  or 
cooling.  The  lack  of  hysteresis  is  consistent  with  the 
second  order  nature  of  transition. 

The  sample  is  cooled  with  a  closed-cycle  helium 
refrigerator  [16] .  The  temperature  of  the  sample  is 
measured  using  a  copper  wire  thermometer  [17] .  A 
thermostatted  heater  allows  the  sample  chamber 
temperature  to  be  held  constant  (fluctuations  <1/10  K) 
[13]  at  any  temperature  between  ambient  and  20  K. 

The  crystals  used  were  typically  a  few  microns 
thick,  50- 100  jam  wide,  and  a  few  mm  long.  Because 
the  width  >  thickness,  the  thickness  was  very  difficult 
to  measure,  making  absolute  values  of  the  modulus  very 
uncertain;  the  small  thickness  also  made  crystal 
mounting  difficult.  The  measured  Q's  of  the  resonances, 
typically  ~  1 000,  were  strain  dependent  and  probably 
limited  by  friction  at  the  clamp. 

For  crystal  No.  1 ,  the  length  equals  0.49  cm,  the 
thickness  (3.3  ±  1 .0)  jam,  and  the  room  temperature 
resonant  frequency  1 1 0  Hz,  giving  a  modulus  value 
E  =  (1 .5  ±  1.0)  x  1012  dyne  cm-2.  Crystal  No.  2  was 
0.35  cm  long  with  a  resonant  frequency  =  670  Hz;  the 
thickness  was  not  determined.  A  third  crystal  (which 
was  not  cooled)  was  0.16  cm  long,  2.8  ±  0.7  fim  thick, 
and  had  a  resonant  frequency  =  1 660  Hz,  giving  a 
modulus  value  E  =  (5.5  ±  2.8)  x  1012  dyne  cm'2.  The 
disparity  of  modulus  values  may  be  due  to  the  failure 
to  obtain  a  rigid  clamp  (which  would  lower  the  resonant 
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Fig.  3.  XY  recorder  trace  for  both  samples  near  the  7j 
transition;  increased  modulus  is  toward  bottom  of 
figure.  An  upper  limit  of  10'4  is  obtained  for  the 
anomaly  size  AE/E0. 

frequency),  or  slightly  warping  a  crystal  when  mounting 
(which  would  increase  the  frequency).  (These 
conditions  could  also  explain  the  slightly  different 
slopes,  din E/dT,  observed  in  crystals  Nos.  1  and  2 
when  cooled,  although  this  might  also  be  a  frequency 
dependent  effect  [9] .)  These  ambient  values  of  the 
modulus  should  be  compared  to  that  of  1.8  x  1012  dyne 
cm'2  for  2H-NbSe2  [9] . 

3.  RESULTS 

Figure  1  shows  the  change  in  E  as  a  function  of 
temperature  (taking  the  room  temperature  value  as  the 
base  line).  Although  an  anomaly  in  E  appears  at  142  K, 
no  anomaly  or  abrupt  change  in  slope  is  discernible 
near  58  K.  The  curve  for  sample  No.  2  shows  some 
hysteresis  (broken  curve)  near  300  K  between  warming 
and  cooling  runs,  probably  arising  from  creep  in  the 
clamp.  No  hysteresis  was  detected  for  sample  No.  1 . 

In  Fig.  2,  four  runs  through  the  7)  transition  in  sample 
2  are  displayed  (with  the  sign  of  A E  reversed  from  that 
in  Fig.  1).  The  upper  two  curves  were  taken  with  a  d.c. 
bias  =  50  V  and  an  a.c.  drive  of  6  V.  For  the  third  run, 
the  a.c.  voltage  was  increased  to  12  V,  and  for  the 
fourth  the  d.c.  was  increased  to  100  V.  The  arrows 
indicate  cooling  or  warming  runs.  The  four  curves 
demonstrate  that  the  anomaly  at  7)  is  not  affected  by 
the  small  strains  and  fields  used,  and  that  there  is  no 
hysteresis  effects  to  the  accuracy  of  our  measurements. 
This  is  consistent  with  non-Ohmic  measurements,  which 
show  that  the  transitions  at  7)  and  7)  are  both  2nd 
order.  [In  2H-TaSe2 ,  the  first-order  commensurate 
transition  near  90  K  shows  [9]  strong  elastic  hysteresis 
(>  5  K).]  The  curves  in  Fig.  2  were  obtained  at  a  drift 
rate  of  ~  0.5  K  min'1.  During  each  run  the  temperature 
was  allowed  to  stabilize  at  several  points  to  eliminate 
the  possibility  of  non-equilibrium  effects  inside  the 
sample  chamber.  Aside  from  the  elastic  anomaly,  the  7\ 
transition  manifests  itself  as  a  slight  change  in  the  slope 


of  E  vs  T.  For  sample  No.  1  we  measure  A(dE/Eo)/dT  = 
+  1 .26  x  1CT4  K'1  across  the  transition,  while  for  sample 
No.  2  the  change  in  relative  slope  is  +  1 .09  x  10'4  K*1. 
The  size  of  the  T  anomaly,  AE/E0  (defined  as  the 
difference  in  the  value  of  E  at  the  minimum  and  the 
extrapolation  of  the  pre-transition  behavior)  is  equal  to 
8.0  x  10'4  (sample  No.  1)  and  9.3  x  10"4  (sample  No.  2). 

In  Fig.  3  we  have  plotted  the  recorder  trace  of  the 
modulus  change  against  temperature  near  58  K  for  both 
samples.  To  the  level  of  uncertainty  in  the  measure¬ 
ments,  &E/E0  ~  10’4,  no  anomaly  or  slope  change  is 
detected. 

As  discussed  above,  the  measured  Q\  of  the 
resonance  were  dominated  by  losses  at  the  clamp.  No 
structure  in  the  Q  was  observed  at  either  transition  for 
either  sample.  An  upper  limit  on  the  internal  friction 
due  to  either  phase  transition  is  tan  6  <3  x  10~4. 


4.  DISCUSSION 

At  a  second-order  phase  transition,  the  change  in 
Young’s  modulus,  £)  (measured  along  the  ith  direction) 
is  related  to  the  specific  heat  anomaly,  ACp,  by  [14] 

A  F 

(6rc/6o,)2  =  --5-  TJACp  (3) 


where  Oj  is  the  uniaxial  stress  in  the  rth  direction  and  Tc 
is  the  transition  temperature.  Although  the  stress 
dependence  of  7)  and  7j  in  NbSe3  have  not  been 
measured  directly,  it  is  of  interest  to  compare  the  value 
obtained  from  equation  (3)  with  published  values  of 
5  TJ&p,  where  p  is  the  hydrostatic  pressure.  The 
pressure  dependence  is  related  to  the  stress  dependence 

(«£  +  *£+«£)  (4) 
,6  oa  8ob  8ac) 


8Tcl'&p  = 


_  _  'l 


where  a,  b  and  c  are  the  three  crystallographic 
directions.  Monceau  and  co-workers  [18]  have  measured 
bTc/8p  to  be  —  4  K  kbar"1  at  both  7)  and  7) .  The 
specific  heat  anomaly  at  T,  has  not  been  determined  but 
is  known  to  be  smaller  than  the  reported  value  of  6  x 
10s  erg  cm'5  K'1  for  the  T2  transition  [19],  Using  this 
value  and  the  value  of  AEb/Eb  reported  here,  we  obtain 
from  equation  (3)  a  lower  bound  for  |8T|/8ffb|.  (taking 
E  ~  3  x  1012  dyne  cm'2), 

\8TJ8ab\>  0.26  K  kbar'1.  (5) 


Assuming  67) /80(,  is  positive,  equation  (5)  places  an 
upper  limit  to  the  anisotropy  in  the  stress  dependence 
of  7), 

k*i0+K.c)<6.4,  (6) 

where 

s  |57’I/5oj|/|57’,/6a6|. 


(7) 
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More  accurate  measurements  to  obtain  ACp  at  7",  are 
necessary  to  get  specific  values  for  and  Af,c. 

The  present  measurements  detect  no  modulus 
anomaly  at  T2.  Thus  we  have  an  upper  bound  of 
A£J,//fb  <  10-4  at  the  lower  transition.  Using  the 
reported  specific  heat  anomaly  value  of  0.5  R,  we  may 
compute  an  upper  bound  for  &T2&ob : 

l6T2/8obl  <  0.05 1  K  kbar'1 .  (8) 

This  implies  that  the  stress  dependence  of  7j  is  highly 
anisotropic.  With  the  same  assumptions  as  those  that 
led  to  equation  (7),  we  have 

4(*a.+*ae)>  37.  (9) 

In  equation  (9),  K2a  is  the  Tj-analogue  of  the  anisotropy 
factor  defined  in  equation  (7). 

5.  CONCLUSION 

By  measuring  the  flexural  resonance  frequency  of 


two  NbSe3  crystals  clamped  at  one  end,  we  have 
observed  an  anomaly  in  the  Young’s  modulus  near  Tt. 
The  absence  of  elastic  hysteresis  is  consistent  with  the 
second  order  nature  of  the  transition.  From  the  observed 
AEf, /Eb,  we  have  observed  an  upper  bound  [equation 
(7)]  to  the  anisotropy  of  the  stress  dependence  of  7\. 
The  refining  of  this  value  must  await  more  accurate 
measurements  on  the  specific  heat  near  7j .  The  absence 
of  any  manifestation  of  the  7j-transition  has  been 
interpreted  as  a  lower  bound  to  the  anisotropy  of  the 
/^-stress  dependence.  It  appears  that  this  anisotropy  is 
an  order  of  magnitude  larger  than  that  at  7),  and  may 
correlate  with  the  larger  size  of  the  resistivity  anomaly 
at  T,. 
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NOTE  ADDED  IN  PROOF 

Superlattice  spots  have  been  recently  observed  in  NbSe3  by  electron  diffraction  [19],  Recent  efforts  to  measure  the 
specific  heat  show  that  the  T2  anomaly  is  much  smaller  than  the  reported  value  in  [18].  This  alters  our  conclusions 
in  equations  (8)  and  (9). 
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The  non-Ohmic  transport  properties  of  NbSe,  are  studied  near  the  two  phase  transitions  at  T,  and  7y 
From  the  non-Ohmic  data  a  parameter  a  is  obtained  which  provides  information  on  the  fraction  of  Fermi 
surface  affected  by  the  phase  transitions.  By  relating  a  to  a  BCS-type  gap  and  modifying  the  calculations  of 
Fedders  and  Martin  on  the  spin-density-wave  state  in  Cr  we  have  obtained  a  theoretical  curve  for  a  which 
agrees  well  with  the  experimental  data.  This  good  agreement  provides  new  indirect  evidence  for  the  charge- 
density-wave  (CDW)  state  in  NbSej.  The  derivative  of  the  low-field  resistivity  versus  temperature  has  also 
been  measured  and  is  shown  to  be  divergent  at  T ,  and  7Y  The  index  of  this  divergence  is  presented  and  its 
significance  to  the  dimensionality  of  the  system  is  discussed.  Finally,  we  briefly  discuss  the  possibility  of 
explaining  the  non-Ohmicity  below  f,  in  terms  of  nonlinear  excitations  in  a  CDW  condensate. 


I.  INTRODUCTION 

The  transition-metal  trichalcogenide1-3  NbSe3 
undergoes  two  phase  transitions  at  the  temper¬ 
atures  Tt  (142  K)  and  T2  (58  K).  The  electronic 
conductivity  is  greatly  perturbed  by  processes 
which  set  in  at  T,  and  T2,  as  evidenced  by  the  ap¬ 
pearance  of  two  giant  anomalies  in  the  resistivity. 
The  conductivity  within  these  anomalies  has  been 
shown3" *  to  be  highly  non-Ohmic,  with  the  conduc¬ 
tivity  increasing  rapidly  as  the  applied  electric 
field  is  increased.  The  frequency  dependence  of 
the  conductivity  is  also  anomalous.3  At  a  frequency 
of  9  GHz,  no  evidence  of  the  Tz  transition  appears 
in  the  microwave  resistivity  while  the  7\  anomaly 
is  sharply  reduced  in  size  compared  to  the  dc  case. 
In  conjunction  with  other  indirect  experimental 
evidence  which  will  be  discussed  in  Sec.  IV A,  the 
two  transitions  have  been  tentatively  identified  with 
the  formation  of  charge-density  waves  (CDW’s). 

If  this  hypothesis  is  valid,  the  non-Ohmic  behavior 
and  the  anomalous  frequency  dependence  raise  the 
intriguing  possibility  that  nonlinear  excitations  are 
being  observed  in  a  CDW  condensate.  Such  exci¬ 
tations  have  been  the  subject  of  much  theoretical 
discussion.  In  particular,  Rice,  Bishop,  Krum- 
hansl  and  Trullinger5  have  described  the  excitation 
of  solitionlike  <t>  particles  in  a  one-dimensional 
(ID)  CDW  system.  Maki6  has  also  discussed  the 
electric- field  generation  of  solitons  as  a  possible 
source  of  non-Ohmicity.  The  lack  of  direct  experi¬ 
mental  evidence  (from  x-ray  or  neutron  scattering) 
for  superlattice  formation  in  NbSe3  makes  the  CDW 
model  rather  provisional  and  it  is  entirely  possible 
that  the  phase  transitions  are  caused  by  a  different  in¬ 
stability.  (See  note  added  in  proof.)  However,  all  ex¬ 
perimental  evidence  obtained  to  date  is  consistent 
with  the  CDW  model  and  recent  Hall-effect  measure¬ 
ments7  have  provided  direct  evidence  of  a  large  loss 


of  carriers  at  both  Tt  and  Tv  This  is  to  be  expected 
from  a  theory  in  which  sections  of  the  Fermi  surface 
(FS)  are  destroyed  by  gap  formation. 

In  this  paper  the  behavior  near  T,  and  T2  of  the 
various  parameters  associated  with  the  non-Oh¬ 
micity  are  presented  and  the  results  are  analyzed 
using  a  model  in  which  a  fraction  of  the  carriers 
are  removed  by  gap  formation  at  the  FS.  The  tem¬ 
perature  dependence  of  the  parameter  a  (defined 
in  Eq.  (2)]  is  interpreted  in  terms  of  thermal  ex¬ 
citations  of  carriers  across  a  growing  gap.  An 
equation  is  derived  relating  a  to  the  gap  function 
A  and  good  agreement  between  experiment  and  the¬ 
ory  is  obtained,  using  appropriately  scaled  Bar- 
deen-Cooper-Schrieffer  (BCS)B  values  for  the  gap 
A.  To  obtain  the  transition  temperatures,  the  de¬ 
rivative  of  the  low-electric-field  resistivity  with  re¬ 
spect  to  temperature  has  been  measured  and  shown 
to  diverge  at  the  two  transitions.  The  critical  in¬ 
dex  of  this  divergence  provides  an  indication  of  the 
conductivity  anisotropy  but  the  conclusions  are  not 
firmly  established  because  of  the  rounding  off  of 
the  divergence  very  near  Tv  and  T2  in  the  present 
data. 

Below  T,  and  T2,  the  non-Ohmicity  of  the  conduc¬ 
tivity  obeys  the  equation3 

o(E,T)  =  a„(T)  +  o*,(T)e-W  =1,2,  (1) 

where  E  is  the  applied  electric  field,  T  is  the  tem¬ 
perature,  and  the  subscripts  1  and  2  identify  the 
transitions.  The  activation  field  E0„  is  strongly 
temperature  dependent.  Its  temperature  depen¬ 
dence  as  well  as  that  of  va  and  have  been  de¬ 
scribed  in  Ref.  3,  where  the  parameter  <*„  defined 
by 

«»  =  **./( +  8  =  1,2  (2) 

was  introduced.  Arguing  from  the  fact  that  the 
conductivity  is  proportional  to  the  FS  area  we  in- 
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terpreted3  a„  as  measuring  the  fraction  of  FS  de¬ 
stroyed  by  gap  formation.  orn  was  shown  to  grow 
from  zero  at  Tn  and  saturate  to  a  constant  value  at 
temperatures  below  that  at  which  the  low-field  re¬ 
sistivity  peaks.  However,  insofar  as  the  above  in¬ 
terpretation  of  an  is  not  unique,  we  shall  discuss 
a  slightly  different  interpretation  of  an  here.  The 
new  interpretation  considers  the  a„  (suitably  nor¬ 
malized)  to  measure  the  number  of  carriers  ther¬ 
mally  excited  across  a  growing  gap.  As  opposed 
to  the  interpretation  in  Ref.  3,  the  fraction  of  FS 
destroyed  by  gap  formation  is  considered  to  be 
temperature  independent.  The  loss  of  carriers  is 
still  a  smooth  function  of  temperatures  because  the 
size  of  the  gap  increases  smoothly  from  zero,  and 
thermal  excitation  of  carriers  across  this  gap  de¬ 
creases  in  relation  to  the  size  of  the  gap. 

II.  EXPERIMENTAL 

All  the  measurements  were  done  in  the  four- 
probe  configuration  with  the  leads  attached  to  the 
sample  by  silver  paint.  The  non-Ohmic  data  were 
taken  using  pulsed  currents  with  a  pulse  width  of 
2  psec  and  a  duty  cycle  of  10"4.  Typically  the 
maximum  power  dissipation  across  a  sample  of 
30  fi  was  5  x  10_7W,  This  was  inferred  to  be  well 
within  the  thermal  capacity  of  the  cryogenic  sys¬ 
tem  because  a  dc  current  dissipating  the  same  po¬ 
wer  in  the  sample  revealed  no  signs  of  temperature 
increase  in  the  sample.  Pulse  heights  were  mea¬ 
sured  by  means  of  a  bucking  variable-voltage  refer¬ 
ence  as  described  in  Ref.  3.  The  uncertainty  and 
reproducibility  of  the  measurements  were  1  in  103, 
decreasing  to  1  in  102  for  the  smallest  pulses.  To 
ensure  that  there  was  no  hysteresis  in  the  mea¬ 
sured  conductivity  (as  the  E  field  was  varied  at  a 
fixed  temperature)  dc  measurements  using  a  10-pA 
current  were  performed  before  and  after  each 
pulse  measurement.  In  samples  which  were  inad¬ 
vertently  damaged  by  excessively  heavy  pulses 
(typically  25  mA  into  a  30- O  sample)  the  measured 
conductivity  showed  large  hysteresis  and  was  not 
reproducible.  Microscopic  examination  of  these 
damaged  samples  invariably  showed  the  presence 
of  broken  strands  which  separate  from  the  main 
body  of  the  sample.  For  undamaged  samples  no 
hysteresis  under  field  variation  or  thermal  cycling 
was  observed  to  the  accuracy  of  the  measurements 
and  results  were  reproducible  after  a  warm-up  to 
room  temperature  and  subsequent  cool  down.  The 
temperature  was  measured  by  a  calibrated  silicon 
diode  sensor  and  the  stability  was  maintained  to 
±20  mK  for  the  20  min  required  to  complete  the 
non-Ohmic  measurement  at  each  temperature 
point. 

In  the  measurements  of  the  low- field-resistivity 


derivative  the  temperature  was  allowed  to  drift  at 
a  rate  of  0.1  K/min  pass  the  transitions,  and  the 
voltage  across  the  sample  was  measured  in  70-mK 
steps  with  an  HP  standard  differential  voltmeter. 
The  derivative  was  computed  by  differentiating  the 
best  quadratic  fit  to  the  resistance  at  five  adjacent 
temperature  points.  There  was  a  hysteresis  of 
approximately  0.1  K  between  cooling  and  warming 
runs,  but  the  two  transitions  were  inferred  to  be 
second  order  because  the  displacement  between 
cooling  and  warming  curves  was  temperature  inde¬ 
pendent  to  ±20  mK. 

Ill  RESULTS 

From  Eq.  (1)  we  have  the  conductivity  in  the  two 
limits  of  weak  and  strong  fields 

a(E~  0)  =  <t„,  (3) 

o(E-°°)  =<Jc„  +  cM.  (4) 

In  Eq.  (3),  cran  is  just  the  low-field  conductivity.  If 
Eq.  (1)  is  valid  then  by  plotting  ln(o-aa)  vs  1/E, 
we  should  obtain  a  straight  line.  From  the  slope 
and  intercept  of  this  straight  line  the  values  of  E0 
and  ctb„  may  be  obtained.  The  parameter  a,  is  then 
computed  from  Eq.  (2).  In  Ref.  3  preliminary  mea¬ 
surements  showed  a„  to  increase  from  zero  at  Tn 
and  then  saturate  to  a  temperature- independent 
value  at  lower  temperatures.  E0„  appears  to  attain 
a  large  value  just  below  Tn,  decrease  to  a  mini¬ 
mum,  and  then  increase  rapidly  as  the  temperature 
is  lowered.  In  this  paper  more  extensive  data  on 
a„  and  E0„  are  reported,  in  particular  near  T,  and 
T2.  Figure  1  shows  a2  and  E02  at  temperatures 
near  T2.  a2  (solid  circles)  grows  smoothly  from 
zero  as  (7\,  -  T)1'2  close  to  T2.  The  solid  line 
obeys  the  equation  2.05  (1  -  T/T2),/2.  At  lower 
temperatures,  the  experimental  points  deviate 
from  the  solid  line.  This  behavior,  reminiscent 
of  the  behavior  of  the  gap  function  A  in  the  BCS 
theory  of  superconductivity ,  has  prompted  us  to 
relate  a,  to  A.  This  will  be  discussed  in  Sec.  IV  B. 
The  open  circles  represent  the  measured  activation 
field  Em.  It  may  be  seen  to  diverge  as  T  ap¬ 
proaches  T*.  The  source  of  this  divergence  is  still 
unclear.  The  broken  line  indicates  the  low- field 
resistivity  p  (zero  suppressed).  It  may  be  seen 
that  deviation  from  the  pretransition  linear  behav¬ 
ior  occurs  2  K  above  T2.  This  is  ascribed  to  scat¬ 
tering  from  fluctuations  in  the  gap  in  analogy  to  the 
case  in  itinerant  antiferromagnets  (see  Sec.  IV D). 
Figure  2  displays  the  same  quantities  a,  and  E01 
close  to  the  T ,  transition.  The  larger  error  bars 
reflect  the  greater  uncertainty  in  the  measurements 
caused  by  the  smaller  change  in  a  induced  by  a 
given  field.  (The  activation  field  near  r,  is  about 
an  order  of  magnitude  larger  than  that  near  T2 
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FIG.  1.  Behavior  of  the  parameter  a2  (solid  circles) 
and  activation  field  E02  (open  circles)  near  T2.  The 
solid  line  is  the  mean-field  prediction  of  the  tempera¬ 
ture  dependence  (Tc—  T)1'*.  The  broken  line  is  the  low- 
field  resistivity  (zero  suppressed). 


while  the  maximum  relative  increase  in  a  is  three 
times  smaller  near  Tt.)  Nevertheless,  behavior 
similar  to  that  of  a2  and  E02  is  obtained  for  al  and 
£0l.  The  solid  line  has  the  equation  0.849(1  -  T/ 
T,),/2. 

In  Figs.  3  and  4  the  critical  behavior  of  a„  and 


FIG.  2.  Behavior  of  the  parameter  0|  (solid  circles) 
and  activation  field  E0I  (open  circles)  near  T , .  The 
solid  line  is  mean-field  behavior.  The  broken  line  is 
the  low-field  resistivity  with  zero  suppressed. 
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FIG.  3.  Log-lot;  plot  of  a,  (circles)  and  F02  (triangles) 
versus  the  reduced  temperature  t  =  1  —  T/7%.  The  open 
(solid)  circles  and  triangles  are  for  values  of  (  with  T, 

'  57.94  K  (57.96  K).  From  the  slopeof  the  straight  lines 
the  exponent  of  a2(E02)  is  +0.50  i  0.03  (—0.45  *0.03). 


Etm  is  shown.  The  data  for  both  o„  and  £0lt  at  the 
two  transitions  may  be  fitted  to  straight  lines  in 
the  log- log  plot  against  the  reduced  temperature 
t  =  (T/Tn-  1).  The  slope  of  the  straight  lines  pass¬ 
ing  through  a,  and  a2  (circles)  is 

3^  0.50  *  0.07.  (5a) 

32  =  0.50  *  0.03,  (5b) 

where  ano:|/lfl».  In  anticipation  of  the  discussion 
which  relates  a„  to  the  order  parameter  the  crit¬ 
ical  index  of  a„  has  been  identified9  with  3.  The 
slope  of  the  Eon  lines  is  -0.45*0.05  (for  £01)  and 
-  0.45*  0.03  (for  £02).  Due  to  the  larger  error  bars 
in  the  £0n  data  the  case  for  a  power- law  behavior 
near  Tn  is  less  firmly  established  than  for  a„.  As 
is  well  known  in  the  interpretation  of  such  data  the 
slope  of  the  “best  fit”  is  affected  by  the  choice  of 
the  critical  temperature  Tn.  In  Fig.  3  the  solid 
circles  correspond  to  choosing  Tz  as  57.96  K  while 
the  open  circles  correspond  to  T2-  57.94  K.  How- 


FIG.  4.  Log-log  plot  of  cv j  (circles)  and  Eol  (triangles) 
versus  the  reduced  temperature  /  =  1  -  T/Tt.  From  the 
two  lines  the  exponents  for  cV|  and  E01  arc  +0.50  *0.07 
and —0.45  +  0.05,  respectively. 
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FIG.  5.  Resistivity  (solid  circles)  and  its  derivative 
with  respect  to  temperature  (open  circles)  near  T,. 

T;  is  taken  to  be  where  rfp/rfT  peaks. 

ever,  precision  resistivity  measurements10  on  the 
ferromagnetics  near  the  Neel  point  have  shown 
that  the  derivative  of  the  resistivity  is  divergent 
at  the  Neel  point.  Recently,  similar  measure¬ 
ments11’12  on  the  one-dimensional  CDW  metal  tet- 
rathiafulvalene-tetracyanoquinodimethane  (TTF- 
TCNQ)  and  its  selenium  analog  have  shown  the  ex¬ 
istence  of  a  divergence  in  the  derivative  of  the  re¬ 
sistivity  with  respect  to  temperature.  The  tran¬ 
sition  temperature  Tc  in  these  systems  is  taken  to 
be  where  the  divergence  occurs  and  it  is  found  ex¬ 
perimentally  that  the  derivative  of  the  resistivity 
behaves  as  a  power  law  in  the  reduced  temper¬ 
ature  on  either  side  of  Tc.  The  critical  index  is 
different  on  both  sides  of  the  transition.  As  an  aid 
to  locating  T,  and  T2  in  NbSe3  the  low-(electric)- 
field  resistivity  has  been  measured  at  70  mK  in¬ 
tervals  near  the  two  transitions.  Figures  5  and  6 
show  for  the  T2  and  Tl  transitions  the  derivative 
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FIG.  6.  Resistivity  (solid  circles)  and  its  derivative 
with  respect  to  temperature  (open  circles)  near  T ,. 
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FIG.  7.  Log-log  plot  of  the  derivative  of  ihc  low-field 
resistivity  (with  background  subtracted)  vs  the  re¬ 
duced  temperature  t  =  (T-  The  solid  (open)  cir¬ 

cles  correspond  to  choosing  7"2  equal  to  58.05  K  (77.95 
K).  The  slopes  of  the  two  lines  are  -1,0  and  -1.5. 

(open  circles)  after  the  subtraction  of  the  flat  pre¬ 
transition  background  B.  As  may  be  seen,  the  de¬ 
rivative  peaks  at  a  temperature  below  where  the 
derivative  changes  sign.  In  both  figures  the  solid 
circles  represent  the  low-field  resistivity  (with 
zero  suppressed).  Unfortunately,  substantial 
rounding  of  the  divergence  occurs,  especially  at 
Tu  and  it  is  not  possible  to  ascertain  the  transition 
points  to  better  than  0.1  K.  For  the  T2  transition 
the  divergence  is  at  58.00  ±  0.05  K.  Although  this 
is  consistent  with  the  value  of  T2  (57.96  K)  used  in 
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FIG.  8.  Log-log  plot  of  the  derivative  of  the  resis¬ 
tivity  (with  background  subtracted)  versus  the  reduced 
temperature  near  the  T,  transition.  As  discussed  in 
the  text  the  rounding  or  the  divergence  prevents  an 
unambiguous  measurement  of  the  critical  exponent 
near  T,.  But  at  temperatures  several  K  above  T,,  the 
exponent  is  -2.0.  This  extends  to  about  15  K  above  Tv 
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studying  the  critical  behavior  of  a2  the  uncertainty 
of  100  niK  prevents  an  unambiguous  determination 
of  the  critical  index  of  dp/dT.  In  Fig.  7  the  deriv¬ 
ative  (after  background  subtraction)  is  plotted 
against  the  reduced  temperature  l  in  log- log  scale. 
(The  open  circles  are  for  T2  =  57.95  K  and  the  solid 
circles  are  for  T.2  =  58.05  K.)  The  critical  index 
lies  between  -1.0  and  -1.5  depending  on  the  choice 
of  T2.  The  flattening  of  the  curves  for  (<3x  io*3 
reflects  the  rounding  off  of  the  divergence.  A  sim¬ 
ilar  situation  obtains  for  the  7',  transition.  The 
peak  determined  by  the  derivative  data  occurs  at 
141.0*0.1  K.  This  is  substantially  lower  than  the 
value  of  142  K  used  in  the  analysis  of  the  critical 
behavior  of  ax.  This  discrepancy  is  larger  than 
the  uncertainty  in  both  measurements  and  is  not 
understood.  Figure  8  shows  the  log-log  plot  of  the 
derivative  versus  t.  As  in  Fig.  7  substantial 
rounding  occurs  near  T2  and  it  is  not  possible  to 
extract  a  critical  index  in  this  region.  However, 
for  f  >2  x  10"  the  derivative  obeys  dp/dT ix/"2*0 
over  an  interval  of  15  K  above  Tx.  The  index  is  in¬ 
sensitive  to  any  reasonable  choice  of  Tv  This 
fluctuation  contribution  to  the  resistivity  over  such 
a  wide  range  of  temperature  makes  the  7\  tran¬ 
sition  appear  much  “broader”  than  the  T2  tran¬ 
sition  in  the  resistivity-temperature  profile. 

IV.  THEORETICAL  AND  DISCUSSION 
A.  Charge-density-wave  model 

The  arguments  in  favor  of  the  CDW  model  in 
NbSe3  have  been  reviewed  in  some  detail  in  Ref. 

3.  Briefly,  these  are  based  on  data  from  the  pres¬ 
sure  dependence2  of  7\  and  T2,  the  detection  of  a 
heat-capacity  anomaly2  at  T2,  the  rapid  rise  of  the 
superconducting  transition  temperature13  Tc  under 
pressure,  and  the  extraction  of  the  parameters3 
a„  which  are  consistent  with  gap  formation  at  the 
FS.  Recently,  the  Hall  effect1  and  Young’s  mod¬ 
ulus14  have  been  measured.  The  Hall  constant  RH 
shows  an  abrupt  rise  at  both  7\  and  T2,  indicating 
a  decrease  in  carrier  concentration  at  these  tran¬ 
sitions.  The  elastic  modulus  along  the  chain  direc¬ 
tion  shows  a  small  anomaly  at  T,  (but  none  was  de¬ 
tected  at  T2).  Both  these  experiments  are  consis¬ 
tent  with  the  CDW  model,  since  the  Young- mod¬ 
ulus  measurement  indicates  a  structural  change 
(at  least  at  7\)  while  the  Hall  measurement  pro¬ 
vides  direct  evidence  for  a  loss  of  carriers  assoc¬ 
iated  with  these  structural  transitions.  In  the  face 
of  the  accumulated  indirect  evidence  it  appears 
worthwhile  to  pursue  the  CDW  hypothesis  further 
even  though  more  direct  evidence  (x-ray  or  neu¬ 
tron-scattering  data)  is  not  available. 

We  consider  the  Frohlich  electron-phonon  (e-p) 
Hamiltonian 


+  XI  tiu>~blb; 


SN 


? gicl .JcJ^  +  hV  ■ 
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(7) 


where  and  hj  are  the  creation  operators  for 
electrons  and  phonons,  respectively,  N  is  the 
number  of  ions,  and  w;  is  the  unrenormalized  pho¬ 
non  frequencies.  Within  the  random-phase  approxi¬ 
mation  (RPA)  the  phonon  frequencies  will  be  renor¬ 
malized  15  to 


Cii  =  cu5  (l  + 
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(8) 


In  Eq.  (8)  the  second  term  is  proportional  to  the 
electronic  susceptibility  and  f(e()  is  the  Fermi- 
Dirac  distribution.  Equation  (8)  may  be  derived 
by15  computing  the  lowest-order  contribution  to  the 
phonon  self-energy  and  examining  the  pole  of  the  pho¬ 
non  propagator  computed  within  RPA.  In  most  me¬ 
tals  the  correction  is  small  except  at  (5  =  2icf  where 
the  derivative  of  the  correction  with  respect  to  q 
becomes  infinite  and  Eq.  (8)  gives  rise  to  the  Kohn 
anomaly.16  A  much  larger  effect  may  arise  in  some 
systems  if  band  pathologies  allow  the  susceptibility 
to  diverge  for  certain  values  of  q.  As  first  pointed 
out  by  Lomer,17  the  existence  of  large  portions  of 
FS  that  are  parallel  (nest)  would  cause  the  denom¬ 
inator  of  the  summand  in  Eq.  (8)  to  vanish  over  a 
large  fraction  of  phase  space,  for  the  value  of  q 
(called  the  spanning  vector  Q)  that  brings  about 
this  nesting  condition.  This  has  been  established 
as  the  mechanism  which  stabilizes  spin-density 
wave18  (SDW)  state  in  Cr.  and  is  also  presumably 
the  case  in  the  layered  compounds  where  CDW’s 
occur. 

The  stability  criterion  may  be  taken  as  the  con¬ 
dition  Rg  =  0.  From  Eq.  (8)  this  is  equivalent  to 
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=  0. 


(9) 


Extensive  calculations  on  the  SDW  state  exist  in  the 
literature.  With  a  slight  modification  the  results 
may  be  applied  to  CDW  systems  as  well.  From  the 
paper  by  Fedders  and  Martin19  (FM)  the  stability 
criterion  for  the  SDW  state  is 


R  r  £>)  -F9(kTQj  ’ 


(10) 


V(0)  =  Aite?  /  iSpT ,  (11) 

where  R  is  the  volume,  V(0)  is  the  long-wavelength 
screened  Coulomb  interaction.  ■>  is  a  matrix  ele¬ 
ment,  c  is  the  electronic  charge,  and  6'rT  is  the 
Fermi- Thomas  screening  length.  The  subscripts 
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a ,  h  refer  to  the  electron  and  hole  bands,  respec¬ 
tively.  The  occurrence  of  two  (or  more)  types  of 
carriers  is  indicated  in  NbSe.,  by  the  Hall  effect7 
which  reverses  sign  at  15  K.  It  is  possible  that 
nesting  may  occur  between  hole  and  electron  poc¬ 
kets  as  inCr.  In  any  case,  the  results  we  shall  use 
from  the  FM  calculation  will  not  be  affected  by 
this  assumption.  The  energies  measured  from  the 
FS  are  given  by 

^  =«•„(*-  kc),  (12) 

«*(«-£*(*  + 9) -*>«rk(*-*«),  (13) 

where  i\,vb  are  velocities  at  the  FS  in  the  two 
bands  and  kc  is  an  average  Fermi  momentum.  The 
FM  results  are  analogous  to  those  of  the  BCS  the¬ 
ory9  of  superconductivity.  They  obtain  for  the 


transition  temperature  7e, 

ka  Te  =  (2 y,rk<//ir)e'>nso/l 

(14) 

where 

lqf  -(K2 \  2Kkc -  6  k2)/ 4£2 , 

(15) 

XSDw  »  (r2  V(0)*2/2ir2r],  lny,  =  0.577  ... , 

(16) 

kl  =  kc(E  -  kc),  V=(r„vb),n,  i'  =  ^(t>a  +  nt). 

(17) 

The  assumptions  made  in  evaluating  Eq.  (10)  are 
that  the  surfaces  of  constant  energy  in  both  bands 
are  spherical  and  congruent.  The  first  Brillouin 
zone  has  also  been  approximated  by  a  sphere  of 
radius  H.  FM  also  derive  an  equation  for  the  gap 
function  A.  In  the  two  limits  T  -  0,  T~Tc,&is 
given  by 

MT-  0)  =  wvkB  Tj yjv  =  l.TG(v/v)k B  Tc,  (18) 

4(7- -  r;)  =  3.0G(v/v)(kB  Tc)(l  -  r/TcYn.  (19) 

The  results  in  Eqs.  (14)- (20)  may  be  applied  to 
a  CDW  system  provided  we  change  the  coupling 
constant  \SDW.  Comparing  Eqs.  (9)  and  (10)  we  may 
define  the  equivalent  coupling  constant  given  by 

\cw=g%nk2c/nwQN2«\’.  (20) 

In  Eq.  (19)  the  gap  increases  as  (1  -  r/Te)W2 
near  Tc.  This  may  also  be  derived  from  the  Lan¬ 
dau  expansion20  of  the  free  energy 

tf**,  n  =  a'{T  -  Tc). 

(21) 

Minimizing  (ECDW  -  Ey)  with  respect  to  ip  we  have 

Z*<l,  =  {a'(Tc-T)/h)''2.  (22) 

The  specific- heat  discontinuity  is  given  by 

Cv/n=7efl'2/6  =-2Tc(FCDW-FAf)/(Te-  T)2.  (23) 

Using  FM’s  result19  for  the  decrease  in  the  free 
energy  near  the  transition 


f,™  -  t's  - -{k'iv2/ 4itW8  Ti)[lz(y]*\  (24) 

we  have  for  the  specific-heat  discontinuity 

C„/f2  =  (4/7£  (3)  }(k2cv /v2)(k~B  Tc).  (25) 

B.  a  parameter  and  the  energy  gap 

In  NbSe3  the  non-Ohmic  results  reported  here  in¬ 
dicate  a  strong  correlation  between  the  behavior  of 
a,  in  Figs.  1  and  2  and  the  order  parameter  de¬ 
scribed  by  Eq.  (22).  In  the  following,  arguments 
will  be  developed  to  establish  the  relationship  of 
a  with  the  gap  A.  Although  the  non-Ohmic  behavior 
is  still  not  understood  it  is  possible  to  extract  in¬ 
formation  on  the  FS  by  examining  the  E  =  0  and  E 
= «»  limits  of  the  conductivity.  We  make  the  as¬ 
sumption  that  the  drop  in  conductivity  at  Tx  and  T2 
is  due  to  the  loss  of  FS  area  which  results  from 
gap  formation  over  the  nesting  fraction.  Since  the 
gap  is  very  small  for  temperatures  near  Tc  there 
will  be  substantial  thermal  excitation  of  carriers 
across  the  gap  at  these  temperatures.  As  the  gap 
increases,  these  excitations  will  drop  rapidly  in 
number.  In  the  low  E-field  limit  we  may  express 
the  conductivity  as21  (suppressing  the  n  subscript) 

<T(£-0)=<r*  +  <rc3f(3A),  (26) 

/(3A)  =  (csa  +  l)‘1.  (27) 

In  Eq.  (26),  3=  (kBT)~‘,  A  is  the  gap  size,  and 
(<jc)  represents  the  conductivity  associated  with 
the  FS  area  not  affected  by  the  gap  (destroyed  by 
the  gap).  Comparing  Eq.  (26)  with  Eq.  (3)  we  have 

oc=os  +  2foc.  (28) 

In  the  strong-field  limit  excitations  induced  by  the 
infinite  field  saturate  the  non-Ohmic  mechanism 
and  we  recover  the  pretransition  behavior  [obtained 
by  setting  A  =  0  in  Eq.  (26)].  From  Eq.  (4)  we  have 

cr0  4-cr*  =  cr^,  +crc.  (29) 

Using  Eqs.  (28)  and  (29)  in  Eq.  (2),  a  may  be  writ¬ 
ten 

a  =  l<7c/(crAr  +  crc)Jtanh(5  3A).  (30) 

Insofar  as  ct„  and  ac  have  the  same  temperature 
dependence,  the  coefficient  of  tanh(|  3A)  in  Eq.  (30) 
will  be  only  weakly  temperature  dependent.  Near 
Tc  we  may  expand  [using  Eq.  (19)] 

ff(D/a(0) ^3CA*  1-53(1  -  T/Te)l,\v/v).  (31) 

This  accounts  for  the  behavior  of  a(T)  shown  in 
Figs.  1-4.  In  Eq.  (31),  o(0)  is  the  value  of  a(T)a.t 
T  =  0  and  3C=  (kTc)~'.  Taking  the  experimental  val¬ 
ues  appropriate  to  the  T2  transition 

a2(0)  =  0.63  ±0.01 


(32) 
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FIG.  9.  Temperature  dependence  of  <>2(  T)/o2(0>.  The 
line  is  the  theoretical  expression  tanh  ('/?  \ )  (see  Eq  . 

(30) 1  with  A=2.12  \ H c g .  The  solid  circles  are  the  data 
from  Fig.  1  while  the  open  circles,  which  have  more 
scatter,  are  from  previously  published  measurements 
(Ref.  3)  on  another  sample. 

and 

a2(r)  =  2.05(l  -  T/T2)l/2(  7-T*.  (33) 

we  have  for  {v/v)2  the  value  [using  Eq.  (33)  in  Eq. 

(31) ] 

(f/tJ)2  =  2.12±0.05.  (34) 

Equation  (34)  expresses  the  ratio  of  the  gap  mea¬ 
sured  in  NbSe3  (at  T2)  and  the  BCS  gap.  From  Eqs. 
(18)  and  (34)  we  derive  a  theoretical  value  for  the 
zero-temperature  gap, 

a2(0)  =  3.13k  BT2  =  19  meV.  (35) 

At  the  7\  transition  the  corresponding  numbers  are 

(r/p),  =  2.8,  A,(0)  =  60  meV.  (36) 

In  Fig.  9  we  compare  the  experimental  data  on  a2 
with  the  BCS-FM  model.  The  measured  values  of 
a2(T)/ a2(0)  for  two  samples  are  compared  with  Eq. 
(30).  The  line  is  the  theoretical  curve  tanh(|  3A). 
Values  of  A  were  obtained  by  multiplying  the  values 
of  provided  by  Miihlschlegel22  by  the  scale  fac¬ 
tor  ( v/v)2  given  in  Eq.  (34).  The  solid  circles  are 
the  data  in  Fig.  1  while  the  open  circles  which 
show  more  scatter  are  from  earlier  published  re¬ 
sults3  on  another  sample. 

The  Landau-type  mean- field  theory  predicts  that 
the  electronic  specific  heat  shows  a  discontinuity 
at  T2  with  a  contribution  given  by  Eq.  (25).  Al¬ 
though  thermal  measurements  on  NbSe3  are  sparse 
we  may  venture  a  comparison  with  the  early  results 
of  Chaussy  et  al.2  The  specific  heat  shows  a  large 
anomaly  (6  x  105  erg/cm3  K)  at  52  K.  This  is  to  be 


compared  with  the  1.7  x  105  erg/cm’  K  anomaly  at 
the  CDW  transition23  in  2//-NbSe3.  Using  Eq.  (34) 
in  Eq.  (35)  we  compute  a  discontinuity  in  C  J 9.  of 
9  x  io3  (3  x  104)  erg/cm*  K  assuming  a  carrier  den¬ 
sity  of  102°  cm'3  (1022cm'3),  and  an  effective  mass 
of  0.3  times  the  free-eleetron  mass.2'1  The  elec¬ 
tronic  contribution  is  clearly  too  small  to  account 
for  the  large  observed  anomaly,  and  one  has  to  in¬ 
clude  the  lattice  contribution.  However,  a  large 
lattice  contribution  would  appear  to  be  inconsistent 
with  the  Young’s  modulus  («)  measurements'4  which 
detected  no  change  in  e  to  3  parts  in  104  at  T2. 
More-accurate  measurements  of  the  heat  capacity 
at  both  transitions  are  necessary  before  a  realistic 
estimate  of  the  electronic  contribution  can  be 
made. 

C.  Activation  field  and  nonlinear  excitations 

In  many  ways  the  non-Ohmieity  of  the  conduc¬ 
tivity  below  /’,  is  the  most  interesting  and  puzzling 
property  of  NbSe3.  Efforts  to  interpret  the  ob¬ 
served  behavior  of  a  given  by  Eq.  (1)  have  not  been 
successful.  The  activated  form  of  a  strongly  sug¬ 
gests  an  electric-field  induced  quantum  tunneling 
process.  However,  the  very  small  E-field  values 
(0.1-  1  V/cm)  needed  to  induce  non-Ohmicity  have 
posed  great  difficulties  in  comparing  theoretical 
expressions  with  experimental  data.  The  simplest 
theory  which  reproduces  Eq.  (1)  is  that  of  Zener 
tunneling4  across  a  gap  at  the  FS.  However,  the 
magnitude  of  the  gap  arrived  at  using  measured 
values  of  E0  is  two  orders  of  magnitude  smaller 
than  kBT.  The  E- field  creation  of  soliton- antiso- 
liton  pairs  in  a  one-dimensional  CDW  condensate 
has  been  discussed  by  Maki.6  In  this  model  the 
electric  field  induces  a  tunneling  from  a  region  in 
< p  space  (where  <p  is  the  CDW  phase)  to  another  re¬ 
gion  which  corresponds  to  the  appearance  of  a  so¬ 
liton- antisoliton  pair.  Interpreting  Eq.  (1)  in  terms 
of  Maki’s  mechanism  again  leads  to  values  of  the 
soliton  rest  energy  two  orders  of  magnitude  smal¬ 
ler  than  kBT.  Thus,  the  thermal  creation  of  these 
solitons  would  appear  to  overwhelm  the  quantum 
creation  at  temperatures  near  r,  and  7  2.  Recently, 
Larkin  and  Lee25  (LL)  have  considered  the  problem 
of  a  one- dimensional  CDW  pinned  by  impurities, 
and  the  nonlinear  conductivity  arising  from  the 
quantum  tunneling  between  impurity  sites.  In  the 
strong  pinning  case  they  obtain  for  the  non-Ohmic 
resistivity 

p  «  exp{|2irrF(»i/*/»/)1/2  ln«  ]/eEl2},  «  =  W/»'f, 

(37) 

where  rF  is  the  characteristic  velocity  related  to 
the  elastic  energy  for  deformation  of  0,  E  is  the 
electric  field,  V  is  the  impurity  strength,  and  l  the 
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average  distance  between  impurities.  Applying 
their  calculations  to  NbSe3  LL  conclude  that  the  ob¬ 
served  value  of  £'„  would  require  l  to  be  extremely 
large.  Furthermore,  as  in  the  previous  two  cases, 
thermal  energy  would  dominate  the  energy  supplied 
bv  the  field.  LL  conclude  that  their  impurity-pin¬ 
ning  model  is  inapplicable  to  NbSe3. 

D.  Critical  divergence  of  dp/dT 

The  divergence  of  the  derivative  of  the  resistivity 
in  pure  metals  and  alloys  at  a  phase  transition  ap¬ 
pears  to  be  very  widespread.  In  ferromagnetics10 
such  as  Ni  the  derivative  of  p  diverges  logarithmi¬ 
cally  at  the  critical  point  T„.  In  antiferromagnetic 
materials  the  divergence  has  a  higher  crit¬ 
ical  index).  Fisher  and  Langer26  (FL)  have  pro¬ 
posed  that  the  divergence  is  due  to  critical  scat¬ 
tering  by  the  short-range  spin  fluctuations  and  have 
proposed  a  divergence  |f  |'“  (where  a  is  the  spec¬ 
ific-heat  index)  above  the  Neel  point  T  N.  Below  Ts 
the  scattering  is  dominated  by  the  growing  order 
parameter  and  dp/dt^f26''  (where  .3  is  the  order- 
parameter  critical  index).  This  agrees  with  what 
is  experimentally  observed  in  Ni.  However,  in  an- 
tiferromagnetics  the  divergence  above  7’v  has  a 
higher  critical  index.  Suezaki  and  Mori27  (SM)  have 
pointed  out  that  the  long-range  spin  fluctuations 
are  dominant  (above  TN)  in  the  case  of  antiferro¬ 
magnetics,  where  fluctuations  occur  around  the 
wave  vector  Q,  leading  to  large  angle  scattering 
of  conduction  electrons.  This  is  particularly  true 
in  SDW  metals  such  as  Cr  where  Q=  2kF.  The  mi¬ 
croscopic  calculations  of  Takada28  point  to  the 
same  conclusion,  namely,  long-range  fluctuations 
are  dominant  in  antiferromagnetics  while  short- 
range  fluctuations  are  dominant  in  ferromagnetics. 
Recently  Horn  and  Guidotti12  (HG)  have  studied  the 
divergence  of  dp/dT  in  the  pseudo-one-dimensional 
organic  metals  [TTF-TCNQ  and  tetraselenafulval- 
ene-tetracyanoquinodimethane  (TSeF-TCNQ)]  at  the 
Peierls  transition  Tp.  They  found  that  (above  T p) 
the  critical  index  is  -1.0  along  both  axes  in  TTF- 
TCNQ  and  -1.5  along  the  6  axis  in  TSeF-TCNQ. 
These  indices  are  substantially  larger  than  those 
found  in  the  magnetic  materials  (0.0  to  -0.3).  To 
explain  their  results  they  considered  the  effect  on 
dp/dT  of  restricting  the  allowed  final  states  for 
scattering  due  to  the  reduced  dimensionality  d  of 
the  metal.  Using  mean- field  values  for  the  critical 
exponents  HG  obtain  the  results27 

w  =  i>, 

w  =  2),  (38) 

(r0-5  {d  =  3). 

The  derivative  of  the  resistivity  arising  from  the 
long-range  fluctuations  may  be  written  as  (cf.  Eq. 


(3.9)  of  Ref.  27] 

%rT  <*  QUQ, A'es  £/;(i  -/a.j) 

x  Re/  dl  exp(-iu)gt^;f  -  17) 

O 

(39) 

k 

where  Q  is  the  spanning  vector  of  the  CDW,  g$  is 
the  electron-phonon  coupling  constant,  /•  is  the 
Fermi-Dirac  distribution,  T'1  is  the  lifetime  of  the 
electrons,  and 

—  (40) 

In  Eq.  (40),  xr  is  the  density- density  correlation 
function  which  diverges  at  7  =  Tc  for  K  =  Q.  To  ev¬ 
aluate  Eq.  (39)  the  assumption  of  critical  slowing 
down  is  made  so  that  we  may  neglect21  the  time  de¬ 
pendence  of  I'j  in  Eq.  (40)  and  it  is  valid  to  use  the 
static  scaling  laws,9 

X5.k(0  )  =  PT(*/k).  (41) 

In  Eq.  (41),  y  is  the  susceptibility  critical  index, 

F  is  an  unknown  function,  and  k  is  the  inverse  cor¬ 
relation  length  which  behaves  as  (for  T  -  Te) 

«ocr.  (42) 

The  dimensionality  d  of  the  electron  gas  enters 
when  we  evaluate  £j*g.£.  Using  Eqs.  (41)  in  Eq. 
(40)  and  performing  the  sum  over  k  we  have 

E*«.|(0)«  f  dxx*G(x).  (43) 

l  O 

Substituting  Eq.  (43)  in  Eq.  (39)  we  finally  have  [us¬ 
ing  Eq.  (42)] 

*  OW'"*  ■  (44) 

Using  mean- field  values  (y  =  1,  v  =  j)  we  have  HG’s 
results  [Eq.  (38)]. 

In  NbSe3  the  critical  exponent  lies  between  -1  and 
-1.5  at  the  T2  transition.  Using  the  above  theory 
this  implies  that  the  dimensionality  of  the  system 
is  2  or  1,  the  data  being  insufficiently  unambig¬ 
uous  to  decide  between  the  two.  However,  a  value 
of  -0.5  or  smaller  for  the  critical  exponent  can  be 
excluded  rather  definitely.  At  the  Tt  transition  the 
rounding  of  the  divergence  very  near  Tl  prevents 
an  extraction  of  the  index.  For  /  >  2  x  10"2,  how¬ 
ever  ,  an  exponent  of  2  is  obtained.  This  is  larger 
than  the  result  in  Eq.  (44)  using  reasonable  values 
of  y,  v,  and  d.  It  is  possible  that  the  assumptions 
underlying  the  derivation  of  Eq.  (44)  no  longer  hold 
at  temperatures  well  above  7\.  In  particular  the 
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critical  slowing-down  assumption  is  expected  to  be 
invalid  far  above  7",  and  one  needs  to  examine  the 
dynamics  of  the  fluctuations  without  this  simpli¬ 
fying  assumption. 

V.  SUMMARY  AND  CONCLUSION 

The  non-Ohmic  conductivity  of  NbSe3  has  been 
measured  at  temperatures  in  the  vicinity  of  the  two 
phase  transitions  at  T,  and  T,.  At  each  temper¬ 
ature  the  conductivity  may  be  decomposed  into  a 
field- independent  part  oa  and  a  field-dependent  part 
obe'E°/E  (suppressing  the  subscript  n,  which  iden¬ 
tifies  the  transition).  From  cra  and  a  parameter 
a  is  obtained  which  (at  T  =  0  K)  gives  the  fraction 
of  FS  affected  by  the  gap.  (In  a  multiband  model 
the  pockets  of  FS  which  do  not  participate  in  the 
CDW  formation  are  included  in  cr4.)  At  elevated 
temperatures  (particularly  a  few  K  below  Tx  and  T2) 
thermal  excitations  across  the  growing  gap  reduce 
the  value  of  a  from  its  T  =  0  K  value.  Therefore 
the  temperature  dependence  of  a  provides  infor¬ 
mation  on  the  temperature  dependence  of  the  gap. 
By  applying  Fedders  and  Martin’s  theory  to  a  CDW 
system  we  have  theoretically  computed  the  temper¬ 
ature  variation  of  a  using  a  BCS-type  gap.  Good 
agreement  between  the  calculations  and  the  data  is 
obtained.  This  agreement  may  be  interpreted  as 
new  evidence  for  gap  formation  at  the  FS  at  both 
transitions.  The  loss  of  carriers  due  to  gap  for¬ 
mation  has  also  been  observed  directly  in  recent 
Hall  measurements.  From  the  fit  of  a  to  the  BCS 
gap  theoretical  values  for  the  zero  temperature 
gaps  have  been  calculated.  The  magnitude  is  60 
meV  (19  meV)  for  the  gap  corresponding  to  the  7\ 
(Ta)  transition.  The  transport  measurements  re¬ 
ported  here  also  provide  some  information  on  the 
dimensionality  of  the  electronic  band  structure. 

By  studying  the  divergence  of  the  resistivity  de¬ 


rivative  at  T.,  we  conclude  that  the  dimensionality 
is  less  than  three.  In  conjunction  with  the  large 
Hall  effect  and  strong  galvanomagnetic  response 
at  low  temperatures,  and  the  absence  of  a  metal- 
to- insulator  (Peierls)  transition  we  believe  that  the 
dimensionality  is  closer  to  two  than  to  one. 

Although  some  progress  has  been  made  in  inter¬ 
preting  the  non-Ohmic  data  the  actual  mechanism 
for  the  striking  electric  breakdown  of  the  anom¬ 
alies  remains  a  puzzle.  A  number  of  theories 
based  on  quantum  tunneling  are  successful  in  ex¬ 
plaining  the  electric-field  dependence  of  the  con¬ 
ductivity,  but  numerical  comparison  with  the  ex¬ 
perimental  data  seems  to  come  up  against  the 
paradoxical  situation  that  the  energy  supplied  by 
the  electric  field  is  much  smaller  than  the  thermal 
fluctuation  energy  at  T ,  and  Ta. 

Note  added  in  proof.  Electron  diffraction  evi¬ 
dence  of  superlattice  formation  in  NbSe.,  at  T, 
has  been  reported  by  Tsutsumi,  Takagaki,  Yama¬ 
moto,  Shiozaki,  ldo,  Sambongi,  Yamaya,  and 
Abe,  Phys.  Rev.  Lett.  39,  1675  (1977).  X-ray 
diffraction  evidence  for  a  new  superlattice  at  T2 
has  also  been  obtained  by  Fleming  and  co-workers 
(P.  A.  Lee  (private  communication)].  Ayrolles  and 
Roucan  have  also  observed  the  superlattice  at  7', 
by  electron  diffraction  .  I P.  Monceau  (private 
communication)]. 
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We  present  measurements  of  the  low-field  Hall  constant  R  as  a  function 
of  field  and  temeprature  in  the  linear-chain  metal  NbSe_,  The  data  range 
from  2  K  to  200  K  in  temperature  and  up  to  15  kG  in  field.  At  the  two 
phase  transitions  T^  (142  K)  and  T2  (58  K)  RH  shows  an  abrupt  increase 
in  magnitude  but  no  change  in  sign.  Below  30  K  R  becomes  strongly 
field  dependent.  The  zero-field-limit  changesMfrom  n-type  to  p-type 
at  15  K.  It  is  argued  that  the  results  are  consistent  with  a  two-band 
model  in  which  the  difference  in  population  p  -  n  is  equal  to  3  x  10*® 
cm~3  below  T2*  The  Hall  data  are  also  consistent  with  a  model  in  which 
both  a  loss  of  carriers  and  a  drastic  change  in  the  conductivity  aniso¬ 
tropy  occur  at  T^  and  T2  and  lend  support  to  the  charge-density-wave 
hypothesis  which  has  been  used  to  characterize  the  two  transitions. 


The  linear  chain  metal  niobium  triselenide 
has  rather  unusual  transport  properties.  The 
resistivity-temperature  profile^  shows  two  giant 
anomalies  associated  with  phase  transitions  at 
T^  (142  K)  and  T2  (58  K).  ^Tjie  electronic  conduc¬ 
tivity  is  highly  non-Ohmic  *  at  temperatures 
below  T. .  The  non-Ohmicity  manifest  itself  as 
a  breakdown  of  the  anomalies  as  the  applied 
E-field  is  Increased.  Furthermore  the  resistiv¬ 
ity  is  a  strong  function  of  frequency  at  microwave 
frequencies^,  No  trace  of  the  T2”transition  is 
detected  in  the  resistivity  at  9  GH^,  while  the  T^ 
anomaly  is  sharply  reduced  in  size. 

The  two  phase  transitions  have  been  dis- 
cussed^*-*  in  terms  of  charge-density-wave  (CDW) 
formation^.  In  this  model  Fermi  surface  (FS) 
pathologies  cause  the  electron-phonon  interaction 
to  become  divergent  for  a  particular  phonon  mode 
^  which  softens.  A  superlattice  of  wave  vector 
Q  appears  and  Induces  gaps  at  joints  on  the  FS 
serparated  by  the  wave  vector  Q.  The  gain  in 
electronic  energy  due  to  gap  formation  stabilizes 
the  new  phase,  gap  formation  also  leads  to  a 
loss  of  carriers  and  accounts  for  the  rise  in 
resistivity  at  the  transitions. 

Recently^  a  model  based  on  these  ideas  has 
been  used  to  interpret  the  non-Ohraic  data  near 
the  phase  transitions.  Using  a  modification  of 
Fedders  and  Martins*'  calculations  on  the  spin 
density  wave  transition  we  have  shown  that  the 
temperature  dependence  of  one  of  the  non-Ohmic 
parameters,  a,  can  be  described  by  a  BCS-type 
gap  equation.  The  parameter  ct  also  provides 
a  number  for  the  fraction  of  FS  affected  by 
the  transit tone  at  Tj  and  T2»  The  CDW  formation 
destroys  20%  of  the  FS  at  Ti  and  60%  of  the 
remaining  at  T2* 


The  most  direct  way  to  observe  a  change  in 
carrier  concentration  in  a  conductor  is  to  mea¬ 
sure  the  Hall  effect  as  a  function  of  temper¬ 
ature.  Recently  we  have  obtained  samples  with 
sufficient  width  to  observe  the  Hall  voltage. 

NbSe^  has  a  rather  large  Hall  signal  at  low 
temperatures.  At  2.0  K  the  Hall  constant  R^. 
is  4.1  x  10“ 7  m^/C  (in  the  limit  of  zero  field). 

As  the  field  H  increase  rapidly  saturates  to 
the  value  of  2.3  x  10“®  mVc.  The  sign  of 
changes  from  p-type  to  n-type  at  15  K  aS  T  is 
increased.  At  both  T.  and  T ^  ^  shows  an  abrupt 
change  in  magnitude  (out  not  in  sign).  remains 
n-type  to  room  temperature  where  it  is  of  the 
order  of  -1.0  x  10“®  m^/C. 

The  samples  were  synthesized  by  heating 
stoichiometric  proportions  of  Nb  and  Se  in  a 
quartz  tube  for  4  weeks  at  725°C  as  discussed 
in  Ref.  (1).  Four  large  crystals  were  found 
suspended  in  the  middle  of  the  tube  away  from 
the  walls.  Microscopic  examination  of  these 
samples  showed  them  to  be  single  crystals 
typically  8  x  0.40  x  0.008  mm*  in  size.  Two  of 
these  samples  were  mounted  on  sapphire  substrates 
and  gold  leads  were  attached  to  them  in  the  stan¬ 
dard  Hall  configuration  by  means  of  silver  paint. 
Temperature  was  monitored  and  automatically 
controlled  by  a  capacitance  sensor  to  +  20  mK 
from  2  K  to  200  K.  A  calibrated  silicon  diode 
sensor  was  used  to  measure  the  absolute  temper¬ 
ature  when  the  magnetic  field  was  off.  The 
longitudinal  dc  current  was  maintained  at  300uA 
and  the  signal  from  the  Hall  leads  was  bucked  by 
a  constant  voltage,  amplified  and  displayed  on 
an  x-y  recorder  as  a  function  of  field.  At  all 
temperatures  the  Hall  signal  for  normal  and 
reversed  fields  was  measured  for  one  direction 
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rig.  1.  The  Hall  voltage  of  sample  1  at  selected 
temperatures  above  Tj.  The  two  branches 
at  each  temperature  correspond  to  the 
two  directions  of  the  applied  magnetic 
field.  The  Hall  constant  RH  is  field 
independent  up  to  15  kG  at  temperatures 
above  AO  K. 


of  the  current.  At  some  selected  temperatures 
the  current  direction  was  also  reversed  to  correct 
for  the  IR  drop  due  to  misalignment  of  the  Hall 
leads.  This  correction  was  necessary  for  only 
one  of  the  two  samples  measured.  Typical  sweep 
rates  were  4  kG/min.  At  high  temperatures  the 
large  amplification  used  required  some  care  in 
eliminating  spurious  effects.  The  sweep  rate 
was  varied  to  detect  time  dependent  emf?s  due 
to  varying  flux  threading  the  Hall  leads.  In 
general  only  those  runs  were  retained,  in  which 
both  the  voltage  across  the  Hall  leads  and  the 
voltage  from  the  Si  cryogenic  sensor  returned 
to  their  initial  values  after  the  field  was 
cycled. 

In  both  samples  the  current  I  was  parallel 
to  the  b  axis  and  the  Hall  voltage  V  was  measured 
along  the  a  axis  with  the  magnetic  field  either 
parallel  or  anti-parallel  to  b  x  'a.  Figure  1 
shows  typical  runs  at  selected  temperatures 


2.  The  hall  voltage  of  sample  2  at  tempera¬ 
tures  below  30  K,  showing  the  strong 
field  dependence  u)  the  Hall  constant 
and  the  reversal  of  sign  that  occurs 
below  20  k.  Similar  behavior  is  obtained 
in  sample  1 . 


for  sample  1.  The  Hall  constant 

Rh  =  (VHt)/(IH)  U) 

(where  t  is  the  sample  thickness  and  H  the 
magnetic  field)  is  field  independent  up  to  the 
highest  field  used  (15  kG)  for  temperatures 
above  30  K.  Below  30  K  R^  is  strongly  field 
dependent.  As  shown  in  Figure  2  (data  on  sample 

2)  V  is  strikingly  non-linear  in  the  vicinity 
_  H 

ot  15  K  where  changes  sign. 

In  Figure  3  we  show  the  temperature  depen¬ 
dence  of  (0)  at  both  T]  and  T?.  Just  ’^elow 
T,  R  (0)  increases  ut  an  inital  rate  of 

»  H 


1 


dR^tOl 


, -  I 


M0) 


dT 


-  0.14  K  (T->T|  ) 


(2) 


As  T  decreases  continues  to  rise  with  a 
decreasing  slope  until  T  crosses  the  T2  tran¬ 
sition  where  R  (0)  again  undergoes  an  abrupt 
increase.  Just  below  the  T^  transition  we  have 


J_ 

R  *(0) 


dRH(0) 


-0.58  K 


-1 


(t+t2) 


(3) 


which  is  four  times  larger  than  the  fractional 
increase  at  Tj.  Below  30  K  the  f ield- dependence 
of  R.j  becomes  manifest  at  the  fields  we  employed. 
The  values  dispiaveu  in  Figure  3  are  the  values 


Vo l .  26,  No.  8 


HALL  EFFECT  OF  A  LINEAR-CHAIN  METAL:  NbSe3 


489 


TOO 


Fig. 3.  The  temperature  dependence  of  Rfj(O)  (the  Hall  constant  at  zero  field). 
The  open  (solid)  circles  are  data  from  sample  1  (2).  The  data  above 
50  K  is  also  plotted  with  a  factor  of  ten  magnification.  The  data 
from  the  two  samples  were  normalized  to  agree  at  2  K  and  14  kG.  At 
both  transitions  T  and  T  Rfj(O)  undergoes  an  abrupt  change  in  magni¬ 
tude.  It  changes  sign  at  15  K,  becoming  p-type  at  low  temperatures. 

The  jump  in  magnitude  at  T  and  T  is  consistent  with  a  change  in  both 
carrier  concentration  and  mobility  ratios  at  the  transitions.  The 
solid  lines  are  drawn  to  guide  the  eye  and  have  no  further  significance. 
The  broken  line  is  the  longitudinal  resistivity. 


of  obtained  by  extrapolating  to  zero  field 
in  a  plot  of  vs.  H.  R^O)  changes  sign  at 
15  K  and  the  Hall  signal  is  p-type  at  low  temper¬ 
atures.  The  strong  field-dependence  of  R^  at 
these  temperatures  is  shown  in  Figure  4  for 
several  temperatures.  The  qualitative  features 
of  R^(0)  “gy  *>e  understood  by  examining  the 
expression0  for  appropriate  to  a  two-carrier 
model.  We  assume  that  the  conductivity  matrix 
is  given  by 

_ _ _  ^  (4) 

CT  =  a  +  a 

n  p 


a 

n 


_ ne _ 

(1  +  U^H2) 


vi  -Uiu2h 

b1U2H  u2 


(5) 


where  subscript  n  (p)  refers  to  electrons 
(holes),  y  (y2)  is  the  mobility  of  electrons 
in  the  x(yj  direction,  and  n  is  the  electron 
concentration.  For  clarity  e  and  y^  are  con¬ 
sidered  to  be  positive  quantities.  The  sign 
of  the  carriers  will  be  written  out  explicitly. 
The  matrix  for  o  is  obtained  from  Eq.  (5)  by 
replacing  n  by  pP(the  hole  concentration)  and  ^ 
y,  by  Vj  (the  hole  mobility).  It  may  be  shown 
that  the  Hall  constant  is  given  by 


RJ)(H)  -  -N/fA^  +  H2  N2) 


(6) 


where 
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In  the  two  limits  of  zero  field  and  infinite 
field  Eq.  (6)  reduces  to 


Rh  (H+o) 


(nyjy2  -  PVjV2) 

efny^pVj)  (ny2+pv2) 
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+  0(H 


3 


(9) 


•s.  (H~>  -  -  ror^r  £  -  0(&\ 


(10) 


For  small  values  of  H,  is  quadratic  in 
H  while  in  the  limit  of  infinite  field  IL 
approaches  the  saturation  value  RH(“>  )  ” 
l/(e(n  -  p))  as  1/H2.  These  features  are 
evident  in  all  the  curves  shown  in  Figure  4. 
Another  prediction  of  Eq.  (6)  is  that  as  the 
mobilities  Increase  (with  decreasing  T)  the 
characteristic  field  at  which  the  behavior  of 
R^  switches  from  H2  to  (1  -  (H-2))  should 
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Fig.  4.  The  field  dependence  of  R^  at  low  temp¬ 
eratures  for  fields  up  to  14  kG.  The 
general  features  at  each  temperature 
are  the  quadrati c  increase  from  its 
zero  field  value  at  low  fields,  and 
the  asymptotic  approach  to  the  infinite- 
field  vaLue  R^(«)  at  high  fields.  As 
discussed  in  the  text,  these  features 
are  consistent  with  a  two  band  model 
in  which  the  mobility  increases  with 
decreasing  temperature.  The  data  are 
from  sample  2. 


move  to  successively  lower  fields.  This  is 
because  H  scales  as  (p^)  or  in 

Eq.  (7).  The  low  temperature  curves  in  Figure 
4  also  show  this  behavior.  For  example  at  the 
lowest  temperature  R  saturates  at  the  rela¬ 
tively  modest  field  of  12  kG.  Although  Eq .  (6) 
can  reproduce  Figure  4  with  a  judicious  choice 


of  the  six  unknowns  n,  p,  a  meaningful 

test  of  the  simple  model  used  in  Eqs.  (4)  and 
(5)  will  be  possible  only  if  the  analysis  also 
agrees  with  ether  transport  measurements  such 
as  conductivity  anisotropy  and  magnetoresis¬ 
tance.  A  detailed  report  will  be  published 
elsewhere. 

We  conclude  with  a  few  remarks.  The 
saturation  value  of  at  high  fields  provides 
unambigously  the  value  of  2.7  x  10*  cm  for 
the  quantity  (p  -  n).  Equation  (10)  implies 
that  at  sufficiently  high  fields  R  attains 
the  same  saturation  value^,  R  (  00  7 ,  at 
all  temperatures  below  T  .  This  behavior 
is  consistent  with  the  low  temperature 

behavior  shown  in  Figure  4.  Above 
T2  changes  in  the  FS  area  alter  (p  -  n)  and 
we  obtain  a  different  R^00  )•  Assuming  a 
characteristic  field  of  M  kG  at  T  *  2  K  we 
deduce  that  the  characteristic  mobility  of  the 
faster  carrier  is  of  the  order  of  10  cm^/Vs 
at  this  temperature.  Below  the  T2  transition 
the  Hall  data  implies  that  p  is  greater  than  n. 
However  the  large  value  of  the  ratio 
keeps  Rh(0)  negative  until  T  decreases  Below 
15  K.  The  sharp  increase  of  R^(0)  at  both  tran¬ 
sitions  is  too  large  to  be  explained  by  a  change 
in  carrier  concentration  alone.  It  is  possible 
that  a  change  in  the  mobility  anisotropies  is 
occuring  simultaneously.  This  will  be  consis¬ 
tent  with  the  appearance  of  a  superlattice 
Bragg  plane  which  preferentially  scatters  the 
carriers  moving  parallel  to  the  spanning 
vector.  In  fact  a  comparison  of  Figure  3 
with  the  prediction  of  Eq.  9  implies  that  the 
conductivity  anisotropy  (nl^  +  pv  )/(nL0  +  pv^) 
should  be  a  strong  function  of  temperature 
below  T^  and  T2.  The  reason  is  that  if  n  and 
p  are  temperature  independent,  a  strongly 
varying  R^O)  implies  that  (Lb)  and  (v  v  ) 
cannot  have  the  same  temperature  dependence 
Preliminary  conductivty  anisotropy  measurements 
using  the  Montgomery  technique  have  confirmed 
this  inference. 

After  this  work  was  completed,  evidence 
for  superlattice  formation  at  T..  was  reported 
by  Tsutsumi  et  al. 
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Conductivity  anisotropy  and  transverse  magnetoresistance  of  NbSe3 
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The  first-transverse  dc  conductivity  measurements  on  a  member  of  the  transiton-melal  tnchalcogemdes 
(NbSej)  are  presented.  The  conductivity  anisotropy  crb/ar ,  (where  b  is  the  chain  axis)  was  obtained  using  the 
Montgomery  technique,  and  varies  from  10  to  20  as  a  function  of  temperature.  The  transverse  resistivity  pr 
also  shows  anomalies  associated  with  the  charge-density-wave  transitions  at  /  (142  K)  and  7;  (58  K) 
However,  these  anomalies  are  smaller  than  the  corresponding  ones  in  pb.  The  quantity  RH(0)/(pb  pr),  where 
Rm  is  the  Hall  constant,  demonstrates  explicitly  the  freeze-out  of  thermally  excited  quasiparticles  below  each 
transition.  We  show  that  at  both  transition  gaps  appears  on  the  hole  surface,  and  the  hole  concentration 
decreases  rapidly.  The  transverse  magnetoresistance  at  low  fields  has  also  been  measured.  Its  monotomc 
power-law  decrease  with  increasing  temperature  shows  explicitly  that  the  resistivity  anomalies  are  due  to 
changes  in  the  carrier  concentration,  and  not  in  the  lifetimes 


I.  INTRODUCTION 

The  transition-metal  chalcogenides  have  proved 
to  be  a  remarkably  fertile  series  of  compounds 
for  the  study  of  instabilities  of  the  electron  gas  of 
reduced  dimensionality.  The  MX,  (where  .1/  is  a 
transition  metal  element  and  -V  is  one  of  the  chal- 
cogens)  compounds  form  layered  structures1  and 
the  electron  gas  has  a  quasi-two-dimensional 
character.  A  host  of  phase  transitions  into  the 
incommensurate  charge-density-wave  (ICDW)  and 
commensurate  charge-density-wave  state  (CCDW) 
occurs  in  the  various  poly -types  of  these  com¬ 
pounds.  Much  theoretical  effort  has  been  expended 
to  study  the  nature  of  the  incommensurate  phase 
and  how  an  ICDW  "locks -in”  to  the  underlying 
host  lattice.  For  instance,  McMillan2  has  pro¬ 
posed  that  the  ICDW  state  in  2H-TaSe2  is  in  fact 
one  in  which  the  charge-density-wave  (CDW) 
achieves  commensurability  over  microdomains 
which  are  separated  by  domain  walls  (discom- 
mensurations)  where  the  CDW  phase  undergoes 
an  abrupt  change.  It  remains  unclear  whether 
such  discommensurations  or  similar  defects 
occur  also  in  other  ICDW  systems  such  as  2H-NbSe2, 
or  whether  the  canonical  CDW  theory  of  Chan  and 
Heine3  is  adequate. 

In  the  last  two  years  compounds  of  the  type 
MX,  have  been  shown  to  exhibit  the  same  kinds  of 
electronically  driven  instabilities.  Specifically, 
the  compounds  NbSe3,4,s  NbS3,6  TaSe,,5’7  and  8TaS3 
have  been  synthesized  and  studied.  The  occurrence 
of  CDW’s  in  the  MX,  family  is  significant  because 
the  members  generally  have  a  linear-chain  struc¬ 
ture  with  high  conductivity  along  the  chain  direc¬ 
tion.  In  contrast  to  the  layered  structure  of  the 
MX,  family  where  the  high  conductivity  occurs 


in  the  basal  plane,  the  MX,  metals  could  possibly 
show  quasi-one-dimensional  behavior.  Indeed, 
Sambongi  et  al.,1  have  shown  that  TaS3  undergoes 
a  transition  into  the  insulating  CCDW  state,  sim¬ 
ilar  to  the  well-studied  quasi-one-dimensional 
conductors,  tetrathiafulvalene  tetracy  anoquino- 
dimethane9  (TTF-TCNQ)  and  the  mixed-valence 
compounds9  KCP  and  their  numerous  derivatives. 
Another  member  TaSe3,  shows  metallic  behavior9 
at  all  temperatures  down  to  2.2  K  where  it  be¬ 
comes  superconducting.8  No  superlattice  forma¬ 
tion  has  been  reported  for  TaSe3. 

In  many  ways  the  most  interesting  member  of 
the  MX 3  family  is  NbSe3  which  was  first  struc¬ 
turally  characterized  by  Meerschaut  and  Rouxel.4 
Monceau  and  co-workers,5  first  showed  the  occur¬ 
ence  of  two  giant  r»  sitivity  anomalies  associated 
with  phase  transitions  at  Tt  (-142  K)  and  T,  (58  K) 
which  they  inferred  to  be  CDW  transitions,  from 
the  pressure  dependence10  of  Tt  and  T,.  Support 
for  this  hypothesis  was  provided  by  measurements 
on  the  Hall  effect11  and  Young’s  Modulus12  which 
showed  respectively  a  loss  of  carriers  at  both  T, 
and  r2,  and  an  elastic  anomaly  at  7t.  Further¬ 
more,  the  growth  of  a  gap  at  the  Fermi  surface 
(FS)  at  both  transitions  was  inferred13  by  applying 
the  BCS  gap  equation  to  data  obtained  from  the 
non-Ohmic  measurements  of  the  longitudinal  con¬ 
ductivity.  Since  then  other  groups14"16  have  ob¬ 
tained  direct  supperlattice  evidence  of  ICDW’s 
in  both  tow-temperature  phases  by  electron  dif¬ 
fraction  and  x-ray  studies. 

The  confirmation  of  the  CDW  state  in  NbSe,  is 
significant  because  of  the  anomalous-electronic 
transport  properties  which  occur  below  Tr  The 
giant  resistivity  anomalies  show  non-Ohmic 
breakdown  effects17  under  weak  dc  electric  fields. 
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The  conductivity  also  shows"  strong  frequency 
dependence  in  the  gigahertz  region  even  when  the 
ae  field  is  much  weaker  than  the  value  required 
to  observe  non-Ohmic  behavior  at  zero  frequency. 
The  great  difficulties  m  accomodating  these  re¬ 
sults  within  the  framework  of  a  smgie-particle 
picture  have  been  discussed  in  Kefs.  13  and  lb. 

In  these  papers  it  was  also  suggested  that  the 
significant  enhancement  ol  the  conductivity  may 
arise  from  the  depinning  of  the  rigid  CDW  con¬ 
densate  which  is  then  accelerated  by  the  K  field. 
This  possibility  has  been  discussed  by  many 
workers"-’0  in  the  context  of  one-dimensional 
“Frohlich  superconductivity.”  Bardeen-1  re¬ 
cently  applied  the  two-iluid  model  to  account  for 
the  temperature  dependence  of  the  activation  field 
in  NbSej,  and  finds  that  (lie  depinning  field  should 
first  decrease,  then  increase  as  the  temperature 
decreases,  in  qualitative  agreement  with  the  ex¬ 
periment.  However,  little  else  is  understood, 
especially  quantitatively,  regarding  the  anomalous 
transport  behavior  aside  from  the  strong  suspi¬ 
cion  that  a  collective  mode  of  the  CDW  condensate 
is  being  excited. 

In  this  paper  we  present  results  on  the  galvano- 
magnetic  and  other  transport  properties  of  NbSe3 
in  the  zero-frequency -Ohmic  regime,  restricting 
current  values  to  below  the  Ohmic  breakdown 
values.  Although  the  nature  of  the  current-carry¬ 
ing  excitations  in  the  non-Ohmic  and/or  high- 
frequency  regime  is  of  paramount  interest,  the 
understanding  derived  from  the  analysis  of  the 
galvanomagnetic  and  conductivity  anisotropy  data 
in  the  Ohmic  region  is  a  necessary  first  step 
towards  an  understanding  of  the  anomalous  be¬ 
havior.  It  will  be  seen  in  the  following  paper22 
(II)  that  the  various  transport  measurements 
taken  together  offer  a  rather  complete  view  of 
the  normal  uncondensed  electrons  and  holes  in 
this  metal.  The  galvanomagnetic  measurements 
in  the  high-frequency  a..J  high  electric  field  re¬ 
gimes  are  currently  in  progress  and  will  be  re¬ 
ported  in  a  later  publication. 

Hall  measurements  on  NbSe3  from  2  to  300  K, 
with  the  B  field  along  a*  and  the  current  parallel 
to  b  (the  chain  axis)  have  been  reported  by  Ong 
and  Monceau11  (OM).  Monceau23  has  measured 
the  Shubnikov-de  Haas  oscillations  in  the  mag¬ 
netoresistance.  Fleming,  Polo,  and  Coleman2'* 
(FPC)  have  also  reported  oscillatory  magnetore¬ 
sistance  measurements  in  Nbi>e3  and  TaSe3,  as 
well  as  high-magnetic-field  Hall  measurements 
in  NbSe3.  A  more  complete  mapping  of  the  Shub- 
nikov  frequencies  has  recently  been  obtained  by 
Monceau  and  Briggs25  (MB).  The  preceding  mea¬ 
surements  are  all  in  substantial  quantitative 
agreement  except  for  the  disagreement  on  Hall 


polarity  reported  by  OM  and  FPC.  Together  with 
tin'  conductivity  anisotropy  in  the  b-c  plane  and 
the  transverse  low-field  niagnetoresistance  mea¬ 
surements  to  be  reported  here  one  should  have 
sufficient  data  to  unravel  the  contributions  from 
the  electron  and  hole  pockets.  In  particular  one 
would  hope  for  a  direct  observation  of  the  loss 
of  carriers  near  the  CDW  transitions,  as  thermally 
excited  quasiparticles  across  the  developing  gaps 
tire  frozen  out  with  decreasing  temperature.  Evi¬ 
dence  for  this  behavior  is  obtained  by  combining 
the  Hall  and  anisotropy  measurements.  A  full 
analysis  of  the  galvanomagnetic  data  is  presented22 
in  II. 

II  EXPERIMENTAL  DETAILS 

Conductivity  anisotropy  measurements  were 
done  using  the  Montgomery26  technique.  The 
fibrous  morphology  of  the  crystals  and  the  minute 
transverse  dimensions  presented  the  most  serious 
difficulties  in  measuring  the  dc  transverse  con¬ 
ductivity.  Samples  were  mounted  on  a  sapphire 
substrate  and  glued  to  the  substrate  by  four 
silver  paint  (du  Pont  4922)  droplets,  which  also 
served  as  contacts  for  the  gold  leads.  This 
mounting  configuration  is  less  satisfactory  in  one 
respect  than  some  that  have  been  reported27  in 
the  literature.  The  differential  contraction  of  the 
sample  and  substrate  sets  up  stresses  in  the  sam¬ 
ple  which  lead  to  cracks  parallel  to  the  chain  axis. 
These  stress-induced  cracks  show  up  as  irrepro- 
ducible  jumps  in  the  measured  transverse  con¬ 
ductivity.  Nonetheless,  this  mounting  configura¬ 
tion  was  used  because  one  could  achieve  better 
control  over  the  size  of  the  contacts  to  the  sample. 
Attempts  to  use  techniques  such  as  that  described 
by  Coleman27  resulted  in  the  droplets  completely 
enveloping  the  thin  sample.  The  Montgomery  con¬ 
figuration  requires  the  contacts  to  be  at  the  four 
corners  of  a  rectangular  sample.  Because  the 
ends  of  single  crystals  tend  to  be  frayed  the  con¬ 
tacts  were  placed  approximately  3  mm  away  from 
the  ends.  Since  the  transverse  dimension  in  the 
<•  direction  was  typically  300  pm  and  the  diameter 
of  the  paint  droplets  was  at  least  100  pm  great 
care  was  required  to  avoid  shorting  out  contacts 
on  opposite  sides  of  the  sample.  The  technique 
which  proved  most  successful  was  to  place  the 
sample  on  the  substrate  and  to  wash  it  with  a  drop 
of  butyl  acetate.  This  improved  adhesion  of  the 
sample  to  the  substrate  as  well  as  the  contact  be¬ 
tween  the  sample  and  paint.  When  the  butyl  ace¬ 
tate  had  evaporated  a  small  droplet  of  paint  was 
placed  on  the  substrate  close  to  the  sample  and 
edged  in  with  a  3-mil  wire  until  contact  was 
achieved.  The  optimum  viscosity  was  maintained 
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FIG.  1.  Conductivity  anisotropy  (in  the  b-c  plane), 
the  transverse  resistivity,  and  the  longitudinal  resis¬ 
tivity  of  NbSe3  (in  descending  order).  The  transverse 
resistivity  has  been  reduced  by  a  factor  of  10  for  con¬ 
venience.  The  longitudinal  resistivity  profile  is  in 
good  agreement  with  previous  results  obtained  by  the 
conventional  four-probe  technique.  Like  pb  the  trans¬ 
verse  resistivity  pc  is  also  greatly  perturbed  by  the 
CDVV  transitions  at  T i  and  T2.  The  inset  shows  the 
actual  sample  orientation  in  the  Montgomery  technique 
used. 


by  diluting  a  fresh  “pot”  of  paint  with  butyl  ace¬ 
tate  for  each  contact.  To  ensure  that  the  sample 
could  withstand  the  stresses  associated  with 
cooling  only  those  crystals  with  substantial  thick¬ 
ness  a  in  the  a*  direction  (normal  to  substrate) 
were  used.  A  successful  cool  down  was  judged 
to  be  one  in  which  no  measurable  hysterisis  was 
observed  on  subsequent  warm  up.  The  amount  of 
stress  could  also  be  estimated  by  comparing  T, 
and  T2  as  well  as  the  size  of  the  two  anomalies  in 
the  longitudinal  resistivity  with  the  ambient  pres¬ 
sure  results.  Chaussy  have  shown  that  these 

parameters  are  sensitive  to  applied  hydrostatic 
pressure.  In  the  results  reported  here  the  long¬ 
itudinal  resistivity  and  the  transition  tempera¬ 
tures  show  no  measurable  difference  from  the 
ambient  pressure  results.  Thus  we  may  safely 
ignore  the  effect  of  these  stresses.  The  sub¬ 
strate  was  thermally  anchored  to  a  copper  plat¬ 
form  which  was  surrounded  by  He  gas  and  en¬ 
closed  in  a  double-vacuum  jacket.  Thermal  sta¬ 
bility  was  maintained  to  ±10  mK  while  the  two 
measurements  were  made.  The  sample  whose 


conductivity  anisotropy  is  shown  in  Fi<-  1  had 
dimensions  1980x290*  50  Mm.  Since  the  contact 
dimensions  varied  from  3'  to  50  pm  there  is 
some  uncertainty  over  what  the  appropriate  effec¬ 
tive  ratio  of  the  sides  of  the  crystal  b/c  should 
be.  This  is  discussed  in  Sec.  III. 

For  the  transverse  magnetoresistance  measure¬ 
ments  much  thinner  fibers  can  be  used  (typically  20- 
pm  wide).  Leads  were  mounted  in  the  conventional 
four-probe  configuration  with  the  B  field  along  the  a* 
axis  and  the  current  in  the  b  direction.  A  quad¬ 
ratic  fit  to  the  resistivity  at  low  fields  was  made 
at  each  temperature.  This  is  important  below 
-12  K  where  the  resistivity  starts  deviating  from 
quadratic  behavior  at  relatively  low  fields.  The 
low-field  magnetoresistance  ±p/p0B?  where  p0 
is  the  zero-field  resistivity  is  very  large  at  4.2  K 
(corresper.ding  to  an  equivalent  mobility  of 
-40  000  rnr’/V  sec),  but  rapidly  decreases  to  an 
equivalent  mobility  of  200  cmW  sec  at  55  K.  No 
magneto-resistance  could  be  measured  to  our  sen¬ 
sitivity  at  fields  of  70  kG  above  the  T2  transition.  Our 
highest  temperature  data  is  at  55  K.  At  1 5  kG  the 
measured  fractional  change  in  the  resistivity  at 
55  K  is  7xi0~4.  Since  the  temperature-coeffi¬ 
cient  of  the  resistivity  is  very  large  at  this  temp¬ 
erature  o'1  dp/tlT?*  -8xiO~2K"‘  the  fractional  in¬ 
crease  in  p  would  be  suppressed  by  a  temperature 
drift  of  9  mK  during  the  time  of  the  field  sweep. 
This  represents  an  upper  limit  on  the  cryogenic 
stability  of  our  apparatus.  The  temperature  of 
the  sample  chamber  was  controlled  by  a  capaci¬ 
tance  sensor  which,  unfortunately,  was  least 
sensitive  around  60  K.  Transients  caused  by  the 
switching  relays  in  the  magnet  supply  were  sup¬ 
pressed  by  careful  screening  of  the  sample  leads. 
Altogether  data  from  three  samples  were  taken. 
Despite  the  scatter  they  are  in  reasonable  agree¬ 
ment  with  each  other. 


III.  RESULTS  AND  ANALYSIS 

In  Montgomery’s26  technique  for  measuring  the 
conductivity  anisotropy  of  a  rectangular  sample 
of  sides  b  and  c,  the  potential  drop  Vb  across  con¬ 
tacts  3  and  4  (see  Fig.  1)  is  first  measured  with 
the  current  /  going  through  contacts  1  and  2.  Next, 
the  potential  drop  Vc  across  contacts  2  and  3  are 
measured  with  the  same  current  going  through 
contacts  1  and  4.  Specializing  to  the  case  where 
the  third  dimension  a  (along  the  a*  axis)  is  neg¬ 
ligible  compared  to  c,  Montgomery  and  Logan 
et  al.2a  show  that  the  resistivity  in  the  b  and  c 
directions  are  given  by 

pb  =H(x)a(Vc/l)(cx/b) ,  (1) 
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•utd  i  is  the  solution  to  the  equation 
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In  practice  at  each  temperature  the  ratio  r  is 
formed  from  the  measured  ls  and  l^.  Equation 
(4)  is  then  solved  '.or  i ,  and  II  computed  from 
Eq.  (3).  Substitution  into  Eqs.  (1)  and  (2)  gives 
p,  and  pc.  For  reference  we  note  that  the  product 

(5) 

is  independent  of  h/c.  To  solve  Eq.  (4)  for  v  we 
have  found  it  convenient  to  use  the  following 
series  which  gives  an  accuracy  of  tO.OlT  when 
l  <  r  <  30. 


x  A  ■  H  lnr  »c'(lnr)-’  +  £>(lnr)5*  •  •  .  (6) 

A  0.997229.  «  0.160697,  C  -- 1.231  448x  10*2, 

/>  -  3.797  40X10"4. 

Figure  1  st  aws  the  c  imputed  anisotropy  in  the 
b-c  plane  and  'he  separate  resistivities  along  the 
two  axes.  (The  transverse  resistivity  pc  has  been 
red  "red  bv  a  lector  of  10.)  As  mentioned  in  Sec. 

1  the  longitudinal  anomaly  sizes  in  pb  show  good 
agreement  with  conventional  four-probe  measure¬ 
ments  on  single  ribbons.  Furthermore,  the  room 
temperature  value  of  pb  (655  pflcm)  is  in  good 
agreement  with  previously  reported5  measure¬ 
ments.  A  rather  surprising  feature  of  the  trans¬ 
verse  resistivity  pe  is  the  sizeable  anomaly  at  7j. 
This  is  unexpected  from  the  orientation  of  the 
superlattice  Bragg  plane  which  is  normal  to  the 
b  axis.  Equally  unexpected  is  the  fact  that  the 
appearance  of  the  CDW  at  T2  which  has  spanning 
vector15  (0.5,  0.26,  0.5)  seems  to  disrupt  pc  less 
thanp4.  Below  T2  the  transverse  resistivity  pc 
decreases  much  more  rapidly  than  pb  but  at  li¬ 
quid-helium  temperatures  the  anisotropy  has 
returned  to  the  value  it  had  at  59  K.  As  men¬ 
tioned  before,  the  computed  anisotropy  pc/pb 
is  very  sensitive  to  the  choice  of  the  ratio  of  the 
sides  h/c.  Since  the  contacts  have  finite  size  it 
is  not  clear  what  values  for  b/c  is  appropriate. 
The  anisotropy  shown  in  Fig.  1  (~16  at  290  K)  was 
computed  with  a  value  of  b/c  equal  to  5.0.  Room- 
temperature  measurements  on  two  other  samples 
gave  values  of  pc/p4  which  varied  between  10  and 
30,  again  depending  on  the  choice  of  b/c.  Thus 
the  absolute  value  of  the  anisotropy  is  known  only 
to  ±501?.  Nonetheless  it  may  be  seen  from  Eqs. 


(4)  and  (5)  that  the  value  of .»  and  (p6pc)  are  un¬ 
affected  by  this  uncertainty.  It  will  turn  out  that 
in  analyzing  the  anisotropy  data  in  conjunction 
with  the  galvanomugnetic  data  the  important  quan¬ 
tity  is  p„pc  which  does  not  depend  on  h/c .  The 
quantity  that  figures  prominently  in  these  analy¬ 
ses  is  ltH/p,,Pc  which  may  be  shown  to  be  indepen¬ 
dent  of  all  sample  dimensions.  (Both  :iH  and 
pbpc  are  linear  inu.)  Similarly  the  magnetore¬ 
sistance  is  independent  of  sample  dimensions. 

Knowing  the  product  pbpc  (solid  circles  in  Fig. 
2),  enables  us  to  proceed  one  more  step  in  under¬ 
standing  the  Hall  effect  of  NbSe3.  The  puzzling 
feature  of  the  Hall  constant  at  zero  field,  R„( 0), 
is  the  large  increase  in  |/<w(0)|  below  7,  and  T2. 
Specifically11  |/ftf(0)|  undergoes  an  eightfold  in¬ 
crease  as  the  temperature  decreases  from  145  to 
60  K  and  appears  to  be  still  increasing  when  the 
7j  transition  is  reached.  A  similar  large  increase 
(12 -fold)  occurs  below  the  T2  transition  (between 
59  and  36  K).  In  a  one-band  model  R„(0)  is  in- 


T  (K) 

FIG.  2.  Product  of  the  resistivities  pbpc  (solid  circles) 
and  the  absolute  value  of  the  quantityA  -  Rfl(Q)/pbpl. 

(open  circles)  vs  temperature  in  log-log  scale.  The 
"resonance”  structure  in  I A  I  is  due  to  the  zero  crossing 
of  Rh( 0)  at  15  K.  The  straight  line  through  the  open 
circles  reflects  the  dominant  power-law  behavior  T'“ 
of  the  mobilities. 


Itt 


CONDUCTIVITY  ANISOTHOPV  AND  TH  A  N  S  V  URSK.  . . 


5269 


versely  proportional  to  the  carrier  concentration. 
However,  it  is  clear  that  such  large  increases  in 
ltH( 0)  with  no  sign  of  saturation  is  incompatible 
with  this  simple  interpretation,  since  the  resis¬ 
tivity  shows  only  a  20%  and  50'"  change  at  7,  and 
T2,  respectively.  Nevertheless,  one  would  hope 
that  the  Hall  data  would  be  amenable  to  a  reason¬ 
ably  straight-forward  interpretation,  and  provide 
experimental  evidence  of  loss  of  carriers  as 
required  by  the  canonical  CDW  theory.3  The 
change  of  sign  of  RH(0)  at  15  K  indicates  that  a 
two-band  model  is  appropriate  for  NbSe3.  This 
is  also  consistent  with  the  very  large  transverse 
magnetoresistance  observed  at  low  temperatures. 
The  expression  for  the  Hall  constant  at  zero  field 
for  the  two-band  model  (assuming  a  closed  simply- 
connected  FS  for  both  holes  and  electrons)  is 
given  by22 

R„(0)  =  -  (np,p2  -pi/li/2)/e(npl  +pvl)(np2  +pu2) , 

(7) 

where  e  is  the  magnitude  of  the  electronic  charge, 
n(p)  are  the  electron  (hole)  density,  pl(p2)  is  the 
electron  mobility  in  the  x(y)  direction  and  ult  v2 
are  the  analogous  hole  mobilities.  The  form  of 
the  denominator  suggests  that  a  simpler  expres¬ 
sion  is  obtained  by  dividing  Eq.  (7)  by  (pji,). 

This  gives 

A^RH(0)/plp2  =  e(pvlu2~nplp2) ,  (8) 

whose  absolute  value  is  plotted  in  Fig.  2  (open 
circles)  in  log-log  scale.  It  becomes  immediately 
apparent  that  the  temperature  dependence  of  A  is 
dominated  by  that  of  the  mobilities  (lifetimes  of 
carriers)  which  have  a  power-law  dependence. 
Dividing  out  this  power-law  dependence  we  thus 
isolate  the  temperature  variation  predominantly 
due  tow  and  p.  In  Fig.  3  the  quantity  A  (r)(T/T2)3-875 
is  shown  versus  temperature.  The  deviation 
above  40  K  of  A  from  the  straight  line  in  Fig.  2 
now  manifests  itself  as  a  sharp  change  in  the  hole 
concentration  relative  to  the  electron  concentra¬ 
tion.  To  be  consistent  with  the  conventional  con¬ 
cepts  of  gap  formation  at  the  FS  we  postulate 
from  Fig.  3  that  the  gapping  occurs  at  the  hole 
surface.  This  is  because  (pvlu2  -nplp2)(T/T2)3-975 
increases  algebraically  as  T2  is  approached  from 
below. 

(The  reader  may  question  our  treatment 
of  the  Hall  effect  due  to  the  thermally  ex¬ 
cited  quasiparticles.  On  the  hole  FS  the  excita¬ 
tion  of  holes  across  the  CDW  gap  leaves  vacan¬ 
cies  behind  which  are  electronlike.  These  will 
presumably  contribute  a  negative  Hall  signal, 
and  to  a  first  approximation,  cancel  the  Hall  sig¬ 
nal  from  the  excited  holes.  However,  this  argu- 


FIG.  3.  Temperature  dependence  of  A  with  the  power- 
law  dependence  T'a  divided  out.  Below  T 2  (open  cir¬ 
cles)  a  equals  3.875.  Above  T2  (solid  circles)  o  has 
the  value  2.792.  It  argued  in  the  text  that  the  increase 
in  A  as  the  temperature  increases  towards  T2  is  due 
to  the  thermal  excitation  of  quasi-particles  across  the 
CDW  gap. 


ment  is  incorrect.  By  studying  the  topology  of  the 
constant  energy  surfaces  it  may  be  shown  that  the 
vacancies  will  always  have  the  same  galvanomag- 
netic  response  as  the  excited  quasiparticles.  A 
more  detailed  treatment  is  offered  elsewhere.) 
Much  the  same  analysis  can  be  applied  to  the  Hall 
data  above  T2.  Again  we  find  that  A  vs  T  in  a  log- 
log  scale  falls  on  a  straight  line  of  slope  -2.792 
except  near  Tv.  Dividing  out  this  power-law  de¬ 
pendence  shows  the  same  carrier  freeze-out 
phenomenon  below  T,  (see  Fig.  4)  rather  more 
clearly  than  in  Fig.  3  because  RH( 0)  has  no  zero 
crossing  above  T2.  Above  T„  A(r/T,)UT44  is 
temperature  independent  indicating  no  relative 
changes  in  the  electron  and  hole  concentrations. 
Thus  the  slight  decrease  in  RH{ 0)  above  145  K 
that  is  evident  in  Fig.  3  of  Ref.  11  is  entirely  due 
to  the  temperature  dependence  of  p,pc.  Similar 
to  the  7^  transition  we  conclude  that  at  T,  the  loss 
of  carriers  occurs  on  the  hole  FS.  From  the 
foregoing  the  picture  that  emerges  is  that  above 
7),  Rh  is  negative  (nplp2>  pv^J.  At  T,  a  gap 
occurs  on  the  hole  FS  associated  with  a  CDW  with 
spanning  vector14  (0, 0.24,  0).  The  decrease  in  p 
relative  to  n  drives  R„  more  negative.  However, 

A  (with  the  temperature  dependence  of  the  life¬ 
time  divided  out)  decreases  smoothly  before 
saturating  instead  of  dropping  abruptly  to  its  sat- 
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FIG.  4.  Temperature  dependence  of  >l(777,1)a  around 
the  upper  transition.  Below  Tt  (solid  circles)  a  equals 
2.792  while  above  T,  (open  circles)  a  equals  1.774. 
Between  ~  100  and  145  K  the  steep  rise  in  A(T/Tl)“  is 
ascribed  to  the  thermal  generation  of  quasiparticles 
across  the  CDW  gap  on  the  hole  Fermi  surface.  Note 
that /4  (T/T,)  "  is  temperature  independent  above  T, 
and  below  100  K. 


uration  value.  This  is  due  to  the  thermal  excita¬ 
tion  of  quasiparticles  (which  are  holes)  across  the 
developing  gap.  Below  -90  K  these  quasiparticles 
have  been  frozen  out  and  A  ('T/T’,)*2-798  becomes 
temperature  independent  until  T2  where  a  new 
gap — again  associated  with  the  hole  surface— ap¬ 
pears  in  conjunction  with  a  CDW  of  spanning  vec¬ 
tor15  (0.5,  0.26,  0.5).  A  is  again  driven  to  more 
negative  values  and  the  appearance  of  thermally 
excited  quasiparticles  (again  holes)  is  shown  in 
Fig.  3.  At  -40  K  the  freeze  out  of  these  carriers 
is  complete.  No  further  change  in  carrier  popula¬ 
tion  occurs  below  40  K.  The  reversal  of  sign  of 
A  is  presumably  due  to  the  steady  increase  of  the 
hole  mobility  versus  the  electron  mobility  at  low 
temperatures.  This  picture  is  borne  out  by  the 
analysis  presented  in  II. 

Figure  5  presents  the  transverse  magnetore¬ 
sistance  data  at  low  magnetic  fields.  We  have  de¬ 
fined  the  effective  magnetoresistance  mobility  p„ 
as  the  square  root  of  the  coefficient  of  B2  in  the 
expansion  of  the  resistivity  in  a  Taylor’s  series, 

p(B)  =p(0)(l  +  p]/Bz  +•••)•  (9) 

Figure  5  shows  the  temperature  dependence  of 
in  log-log  scale.  The  striking  feature  of  the  data, 
in  view  of  the  model  sketched  above,  is  the  mono¬ 
tonic  power-law  decrease  of  as  T2  is  approached 
from  below.  Significantly,  whereas  and  pb  pc 
attain  a  maximum  before  decreasing  to  their 
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FIG.  5.  Transverse  magnetoresistance  expressed  as 
an  effective  mobility  pv  vs  temperature  in  log-log 
scale  [see  Eq.  (9)|.  Data  from  three  samples  are  shown 
Note  that  the  mobility  p  v  continues  to  decrease  w  ith  in¬ 
creasing  temperature  between  40  and  54  K  whereas  the 
longitudinal  conductivity  is  rapidly  rising  in  this  tem¬ 
perature  interval.  This  implies  a  steep  increase  in 
carrier  concentration  as  discussed  in  the  text. 


value  at  T2,  pM  shows  no  such  structure.  Insofar 
as  p."  measures  the  lifetime  of  the  carriers  (in 
an  averaged  way)  this  immediately  implies  that 
the  steep  rise  in  conductivity  in  both  b  and  r  direc 
tions  as  the  temperature  increases  from  40  K  to 
T2  is  not  due  to  any  increase  in  carrier  lifetimes. 
Taken  together  with  the  change  in  carrier  concen¬ 
tration  shown  by  the  data  on  A,  the  data  in  Fig.  5 
give  strong  support  to  the  conventional  picture  of 
FS  gapping  at  a  CDW  transition. 

IV  CONCLUSION 

The  first  measurements  on  the  transverse  dc 
resistivity  of  a  transition-metal  trichalcogenide 
are  presented.  Although  the  linear-chain  crystal 
structure  of  NbSe,  suggests  the  possibility  of 
quasi-one-dimensional  behavior  the  measured 
anisotropy  (-20)  is  several  orders  of  magnitude 
smaller  than  the  best  one-dimensional  conduc¬ 
tors.  Nonetheless,  the  peculiar  transport  prop¬ 
erties  of  this  metal  make  it  directly  relevant  to 
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the  issues  of  interest  in  the  field  of  highly  aniso¬ 
tropic  conductors,  particularly  the  question  of 
the  depinning  of  the  CDW  condensate  and  the 
sliding  FrShlich  mode.  Some  clarification  of  the 
electronic  instabilities  have  been  obtained  by 
studying  the  conductivity  anisotropy  in  relation 
to  the  Hall  constant  and  transverse  magnetore- 
sistance.  In  particular  it  is  demonstrated  that 
the  gapping  occurs  on  the  hole  FS  at  both  transi¬ 
tions.  Furthermore,  the  giant  resistivity  ano¬ 
malies  in  the  longitudinal  resistivity  are  shown 
to  be  a  consequence  of  loss  of  carriers  at  Ti 
and  T2  due  to  the  slow  freeze -out  of  thermally 
excited  quasiparticles  across  a  developing  gap. 
This  confirms  the  previously  published13  analysis 
of  the  non-Ohmic  data  which  also  showed  that  the 
decrease  in  conductivity  was  consistent  with 
carrier  freeze-out  across  a  gap  which  could  be 
described  by  the  BCS  equation.  In  particular,  by 
combining  the  product  pbpc  with  previously  pub¬ 


lished"  data  on  the  zero-field  Hail  constant 
UH(fi)  which  shows  a  surprisingly  large  change  at 
7',  and  we  obtain  an  expression  which  shows 
more  transparently  the  decrease  in  the  hole  con¬ 
centration  below  7',  and  72.  Finally,  the  trans¬ 
verse  magnetoresistance  provides  complementary 
information  on  the  averaged  lifetime  of  the  car¬ 
riers  and  shows  explicitly  that  the  resistivity 
anomaly  at  7j  is  not  due  to  pathological  changes  in 
carrier  mobilities,  as  has  been  suggested  by 
other  workers.24 
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The  charge-density-wave  (CDW)  linear-chain  metal  NbSe,  shows  striking  non-Ohmic  behavior  when  the 
applied  electric  field  exceeds  —0.1  V/cm.  Hall  effect,  transverse  magnetoresistance,  conductivity  anisotropy, 
and  Shubnikov-de  Haas  measurements  using  sufficiently  low  current  densities  to  avoid  Ohmic  breakdown 
have  been  published.  We  propose  a  simple  two-band  model  to  account  for  the  temperature  dependence  of 
these  quantities  as  well  as  the  (magnetic)  field  dependence  of  the  Hall  constant  in  the  Ohmic  regime  below 
58  K.  The  model  has  six  unknowns  (carrier  concentrations  and  mobilities)  that  are  fixed  by  six  experimental 
numbers  at  each  temperature  T.  The  solution  shows  that  all  the  mobilities  obey  a  power-law  behavior  versus 
T,  whereas  the  carrier  concentrations  are  both  T  independent  up  to  40  K.  Above  40  K  the  hole  population 
rises  sharply,  analogous  to  the  theoretical  predictions  for  an  excitonic  insulator.  This  implies  that  the  CDW 
gap  occurs  on  the  hole  surface.  Using  the  parameters  of  the  model,  we  have  recomputed  the  resistivities, 
Hall  constant,  and  magnetoresistance,  and  they  have  been  shown  to  agree  with  all  the  available  experimental 
data.  Thus  the  conventional  single-particle  picture  with  the  additional  hypothesis  of  a  BCS-type  gap  on  the 
hole  surface  is  adequate  for  understanding  the  transport  properties  of  NbSe,  in  the  zero-frequency-Ohmic 
regime.  We  also  interpret  the  SdH  data  in  terms  of  the  two-band  model. 


I.  INTRODUCTION 

Below  the  transition  temperature  T,  (~142  K) 
the  transport  properties  of  the  charge-density- 
wave  (CDW)  metal1  NbSe,  are  conveniently  divided 
into  three  regimes  for  the  purpose  of  analysis. 
These  are,  (i)  the  zero-frequency  low-current 
(Ohmic)  regime;  (ii)  the  zero-frequency  high-cur- 
rent  (non-Ohmic)  regime;  (iii)  the  high-frequency 
(above  'TOO  MHz)  low-current  regime.  (There 
clearly  exists  a  fourth  regime:  the  high-frequency 
non-Ohmic  regime,  but  no  experimental  data  exist 
in  this  regime.)  Strong  interest  has  been  generated 
by  the  study  of  regimes  (ii)  and  (iii)  because  of  the 
possibility  of  studying  experimentally  the  dynamics 
of  a  sliding  CDW  in  a  real  metal.  Such  studies  will 
shed  light  on  (a)  the  nonlinear  dynamics  associated 
with  a  drifting  CDW  ,J  (b)  the  nature  of  the  pinning 
forces  impeding  its  motion,  and  (c)  the  feasibility 
of  realizing  a  system  in  which  the  sliding  Frbhlich 
mode  leads  to  superconductivity.  With  regard  to 
(b)  the  ease  with  which  the  pure  compound  is  al¬ 
loyed3  or  intercalated1’5  with  Li  makes  it  an  at¬ 
tractive  system  in  which  to  study  the  pinning  ef¬ 
fects  of  the  CDW'  condensate  due  to  impurities. 

Such  impurity  studies  are  presently  underway  and 
show1  that  the  behavior  in  regimes  (ii)  and  (iii)  are 
exceedingly  sensitive  to  impurity  concentrations. 

Most  of  the  galvanomagnetic''6  as  well  as  Fermi- 
ology7'9  studies  done  to  date  have  been  in  regime 
(i)  where  the  applied  electric  fields  are  sufficiently 
low  to  avoid  Ohmic  breakdown.  Although  regime 
(i)  is  of  less  interest  as  compared  to  (ii)  and  (iii) 
it  is  essential  to  obtain  a  picture  of  the  “normal" 
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regime  that  is  as  complete  as  experimentally  fea¬ 
sible  so  that  the  intriguing  additional  effects  (con¬ 
ductivity  enhancement)  in  regimes  (ii)  and  (iii)  can 
be  isolated  and  analyzed.  In  this  paper  we  con¬ 
struct  a  simple  two-band  model  that  will  account 
for  the  galvanomagnetic  and  conductivity  aniso¬ 
tropy  data  [regime  (i) ]  presented  in  Refs.  1,  6, 
and  8.  Single  crystals  of  NbSe,  usually  grow  in 
the  shape  of  long  ribbons  with  the  long  side  along 
the  (chain)  b  axis  and  the  wider  transverse  dimen¬ 
sion  along  the  c  axis.  Low-magnetic-field  Hall 
measurements  with  H  along  «*  have  been  reported 
by  Ong  and  Monceau  (OM).r>  High-field-Hall  mea¬ 
surements  in  the  same  geometry  have  been  pre¬ 
sented  by  Fleming.  Polo,  and  Coleman  (FPC).8 
Transverse  magnetoresistance  (J/||o*)  were  pre¬ 
sented  in  the  preceding  paper1  as  well  as  conduc¬ 
tivity  anisotropy  measurements.  In  addition,  ex¬ 
tensive  Shubnikov-de  Haas  (SdH)  measurements 
have  been  carried  out  by  Monceau  (M).7  Fleming 
d  al.  (FPC)6,  and  Monceau  and  Briggs  (MB).9  It 
will  be  seen  that  a  two-band  model  can  satisfactor¬ 
ily  account  for  much  of  the  data  in  these  studies 
below  the  T  ,  transition.  Out  of  this  analysis  we 
obtain  the  temperature  dependence  of  the  four  mo¬ 
bilities  in  the  6-c  plane  as  well  as  the  carrier 
concentration  in  both  bands.  In  support  of  the  con¬ 
ventional  CDW  model,10’"  our  analysis  shows  that 
the  carrier  concentrations  are  temperature  inde¬ 
pendent  below  ~40  K.  Above  40  K  the  hole  popula¬ 
tion  rises  dramatically  in  analogy  with  the  increase 
in  thermally  excited  quasiparticles  across  the 
Fermi  surface  (FS)  gap  of  an  excitonic  insulator.11 
In  contrast  the  electron  population  remains  con- 
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stant  up  to  58  K.  Persuasive  evidence  identifies 
the  electron  pocket  with  the  ellipsoidal  pocket  seen 
in  the  SdH  oscillations  of  M.  FPC,  and  MB.  The 
evidence  for  a  hole  pocket  is  unclear  at  present. 
This  two-band  analysis  enables  us  to  disentangle 
the  individual  contributions  of  the  four  mobilities 
to  the  conductivity  anisotropy,  the  Hall  signal  and 
the  transverse  magnetoresistance.  It  will  be 
shown  for  instance  that  the  twelve-fold  increase1' 
in  |/tw(0)  |  at  7 2  is  compatible  with  the  twofold  rise 
in  the  longitudinal  resistivity  and  the  still  smaller 
rise  in  the  transverse  resistivity.  The  zero 
crossing  of  RH{ 0)  at  ~15  K  is  due  to  the  faster  in¬ 
crease  of  the  hole  mobilities  relative  to  the  elec¬ 
tron  mobilities  as  the  temperature  drops,  and  not 
to  any  change  in  the  carrier  population.  On  the 
other  hand,  the  giant  anomaly  in  the  resistivities 
is  due  to  the  rapid  change  in  hole  concentration 
between  ~40  K  and  T ..  No  anomalous  structure  is 
reflected  in  any  of  the  mobilities.  This  is  also 
presumably  the  case  at  the  upper  transition,  al¬ 
though  the  lack  of  sufficient  galvanomagnetic  data 
and  theoretical  considerations  preclude  the  exten¬ 
sion  of  our  analysis  beyond  58  K. 

II  TWtFBAND  MODEL 

The  single-particle  picture  with  a  semiclassical 
treatment  of  the  transport  properties'3  has  been 
remarkably  successful  when  applied  to  ordinary 
metals  and  semiconductors.  To  make  the  analysis 
manageable  we  have  assumed  the  simplest  Boltz¬ 
mann  equation- relaxation  time  approach  in  com¬ 
puting  the  galvanomagnetic  response.  Since  the 
details  are  found  in  many  text  books1'*  we  will  pre¬ 
sent  results  only.  The  assumptions  intrinsic  to  the 
model  are  as  follows,  (a)  At  low-temperatures 
there  exist  two  closed  FS  pockets  containing  elec¬ 
trons  and  holes,  respectively,  (b)  These  pockets 
are  ellipsoids  with  their  principal  axes  along  a‘ , 
f>,  and  c.  (c)  The  relaxation  time  is  assumed  iso¬ 
tropic  and  independent  of  magnetic  field  and  com¬ 
plications  such  as  intense  scattering  from  “hot 
spots”  on  the  FS  will  be  disregarded.  Assump¬ 
tions  (a)  and  (b)  are  indicated  by  the  SdH  data 
which  show  two  branches.  The  prominent  branch 
is  consistent  with  an  ellipsoid  of  moderate  aniso¬ 
tropy  with  principal  axes  along  a*.  6.  and  c.  Less 
conclusive  evidence  is  available  for  the  other  poc¬ 
ket  which  we  take  to  be  ellipsoidal  as  well.  The 
third  assumption  will  be  discussed  together  with 
the  computed  magnetoresistance.  No  further 
justification  of  the  assumptions  will  be  discussed 
except  to  appeal  to  the  comparison  of  the  compu¬ 
tations  with  the  experiments. 

In  the  Boltzmann  equation  approach  the  current 
is  given  by1'* 
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where  vf  =  l<''9tf  9k  is  the  velocity  of  the  carrier 
at  k,  f  is  the  Fermi-dirac  distribution,  t  the  iso¬ 
tropic  relaxation  time,  m,  the  effective  mass  a- 
long  the  i  axis  and  the  integral  in  (1)  is  over  the 
FS.  In  general  the  relaxation  time  may  be  aniso¬ 
tropic  as  well,  especially  if  the  system  is  highly 
anisotropic.  '  iowever,  the  galvanomagnetic  mea¬ 
surements  analyzed  here  cannot  distinguish  be¬ 
tween  mass  and  relaxation  time  anisotropy.  In  the 
spirit  of  the  model  we  will  describe  all  the  ob¬ 
served  anisotropy  in  terms  of  the  mobilities.  As¬ 
suming  a  quadratic  dispersion  it  is  straightforward 
to  reduce  Eq.  (1)  to  the  familiar  conductivity  ma¬ 
trix  for  electrons  (with  A  along  the  three  axis) 
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where  the  mobility  p,seT,  »n(  has  been  introduced 
and  n  the  carrier  concentration  is  given  by 
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Calling  the  hole  concentration  p  and  the  hole  mo¬ 
bilities  vt  we  have  for  the  total  conductivity  matrix 
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For  convenience  we  have  absorbed  e  into  n  and  p. 
All  quantities  in  Eq.  (5)  are  assumed  positive. 

The  sign  difference  in  the  off-diagonal  elements 
has  been  shown  explicitly.  In  the  experimental 
situation15’1"  it  is  actually  the  resistivity  matrix 
p  =  &  *'  which  is  measured.  Explicitly,  if  the  con¬ 
stant  current  is  along  the  two-axis  (6  axis)  and 
the  magnetic  field  along  the  three-axis  (a*  axis) 
the  electric  field  along  the  one  and  two  axes  is 
given  by  E=p’J  where  3  =  (0,  J„).  The  Hall  voltage 
is  proportional  to  E ,  =  p,_.t/0  whereas  the  transverse 
magnetoresistance  is  measured  by  E.,  =  p22J0.  In¬ 
verting  the  matrix  in  Eq.  (5)  we  may  show  that  the 
Hall  constant  /?„=£,/(#«/„)  is  given  by 


R„(H)  --  -A  B.  +AzHi)  ,  (6) 


where 
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and 

S,^nn,/il  +  (plp,H)-]+pv,/[l  (8) 

In  the  limits  of  zero  field  and  infinite  field  Eq. 

(6)  reduces  to 

RosRm(°)  +pvi)(np2 

(9) 

and 

R~sRw(«)  =  -l/(»-/>).  (10) 

We  have  adopted  the  convention  that  R„  is  negative 
if  there  are  no  holes  present  (p  =  0).  From  Eqs. 

(9)  and  (10)  it  is  clear  that  whereas  the  sign  of 
R x  is  determined  by  the  more  populous  carrier 
type  the  sign  of  R  „  at  low  fields  is  also  influenced 
by  the  relative  sizes  of  the  mobilities  and  the 
anisotropy  in  each  band.  The  longitudinal  resis¬ 
tance  in  the  presence  of  a  weak  magnetic  field  be¬ 
haves  as  p(H)  -p„(l  +  P*#2)  where  p„  is  the  effec¬ 
tive  magnetoresistance  mobility.  From  the  diag¬ 
onal  element  of  the  resistivity  matrix  we  calculate 

_>  P.,.(tf)-p,..(0)  npp,  i/,(m  ,  4-i/2)j 

"  (np,  +pv,)(np2  +pVi)  ' 

We  note  that  p„  vanishes  if  either  «  or  p  is  zero. 
This  is  the  expected  result  in  a  single  band  where 
the  relaxation  time  is  energy  independent.  The 
Hall  field  prevents  any  deflection  of  the  carriers 
and  the  measured  longitudinal  resistance  is  unaf¬ 
fected  by  a  transverse  magnetic  field.  Further¬ 
more  p„  vanishes  if  p.,  -  -t/2,  i.e.,  if  the  two  bands 
have  the  same  sign  of  charge  and  the  same  longi¬ 
tudinal  mobility.  The  longitudinal  mobility  is  re¬ 
sponsible  for  this  cancellation  because  the  longi¬ 
tudinal  velocity  sets  up  the  Lorentz  force.  As 
mentioned  above  the  mobilities  themselves  are 
assumed  to  be  field  independent.  In  principle  the 
relaxation  time  can  have  a  field  dependence  such 
as  t'1  =  t~'  +a*H~.  However,  since  the  observed 
magnetoresistance  at  low  temperatures  is  very 
large  we  shall  assume  that  the  dominant  effect  is 
due  to  the  deflection  of  carrier  trajectories  arising 
from  mutual  (partial)  cancellation  of  the  Hall  field 
of  the  holes  and  electrons.  The  existance  of  hot 
spots  on  the  FS  where  electron-phonon  interaction 
is  enhanced  can  also  affect  the  magnetoresistance. 
But  in  the  low-field  limit  where  the  carriers  sam¬ 
ple  only  a  limited  portion  of  the  FS  beiore  being 
scattered  we  may  neglect  such  complications. 

In  the  model  described  by  Eq.  (5)  there  are  six 
unknowns,  n.  p,  p,,  iq,  p2,  and  v2.  At  each  tem¬ 
perature  there  are  four  measured  quantities-,  the 
Hall  constant  R„(0),  the  transverse  magnetoresis¬ 
tance,  and  the  zero-field  resistivities 


The  Hall  constant  at  infinite  field  R®  which  can  be 
measured  only  at  very  low  temperatures  provides 
a  constraint  through  Eq.  (10).  Finally,  the  field 
dependence  of  R„  at  each  temperature  provides  the 
sixth  equation  to  determine  the  six  unknowns, 
leaving  no  adjustable  parameters.  Although  the 
analytical  inversion  of  the  six  equations  is  rather 
difficult  the  numerical  solution  may  be  readily  ac¬ 
complished  by  first  forming  the  dimensionless 
constants 

C1  =  Pi,<p1P2),/V(fl»-ftffl),  (13) 

C2£P*(p,P2),/2/rt°t>.  (14) 

(It  turns  out  that  C,  and  C2  may  be  measured  very 
accurately  even  if  the  sample  dimensions  are  un¬ 
known.  This  is  because  (p,p2)l/‘  (determined  by 
Montgomery's”  technique)  is  independent'_of  6  and 
c,  and  is  linear  in  a  (the  dimension  along  H).  This 
is  also  true  of  R„(H).  p„  is  independent  of  sample 
dimensions.  Thus  C,  and  C2  are  the  ratio  of  volt¬ 
meter  readings).  A  closed  subset  of  three  equa¬ 
tions  can  be  formed  which  expresses  C,,  C2,  and 
Rh(H)/Rx  in  terms  of  the  three  dimensionless  pa¬ 
rameters 

r-np,  7/  =  p,/tq,  and4  =  p2/iq.  (15) 

Finally,  these  three  equations  can  be  solved  nu¬ 
merically  for  r,  £,  and  rj.  Details  of  the  solution 
are  presented  in  the  appendix.  We  note  that  inso¬ 
far  as  C,,  C2,  and  R„/Rx  do  not  depend  on  sam¬ 
ple  dimensions  the  values  r,  rj.  and  4  will  not  be 
influenced  by  uncertainties  in  sample  dimensions. 
However,  to  compute  p,,  iq,  p2,  and  v2  separately 
we  need  to  use  the  anisotropy  p,/p2  which,  of 
course,  depends  on  the  ratio  b/c. 

In  the  foregoing  analysis  we  have  assumed  that 
the  carrier  concentrations  are  temperature  de¬ 
pendent.  This  enables  us  to  use  the  value  of  R® 
measured  at  2  K  to  fix  the  value  of  ( p  -  n).  (It  is 
not  possible  to  measure  R®  at  elevated  tempera¬ 
tures  with  currently  available  magnetic  fields.) 
This  assumption  is  justified  in  a  self-consistent 
fashion  by  the  solutions  below  36  K  which  show 
that  r  in  Eq.  (15)  is  temperature  independent. 
However,  this  assumption  is  invalid  at  tempera¬ 
tures  above  ~40  K  where  a  previous  analysis18  of 
the  non-Ohmic  conductivity  data  shows  that  sub¬ 
stantial  thermal  excitation  of  quasiparticles 
across  the  CDW  gap  occurs.  This  is  also  clearly 
demonstrated  by  Hall  data  in  Figs.  3  and  4  of  the 
previous  paper.1  In  addition  the  field  dependence 
of  R„(R)  has  not  been  measured  for  temperatures 
above  36  K.  Thus  we  are  deprived  of  two  pieces 
of  experimental  information.  To  continue  the 
analysis  to  higher  temperatures  we  need  to  intro¬ 
duce  a  fourth  assumption;  (d)  the  electron  con- 
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centration  u  remains  temperature  independent 
from  36  to  58  K.  This  implies  that  the  electron 
mobility  pu  decreases  with  the  same  power  law. 

The  remaining  four  unknowns  p,,  v,.  v,.  and/) 
are  then  computed  from  the  measured  p, ,  p_,, 

RH(0),  and  p*.  With  this  fourth  assumption  we 
find  that  (Fig.  1)  the  computed  p,  (and  c,)  deviates 
only  very  slightly  from  its  previous  behavior;  and 
p  rises  dramatically  to  four  times  its  value  at  36 
K.  The  computed  behavior  of  p,  provides  justifi¬ 
cation  for  assumption  (d)  (that  the  electron  pocket 
is  not  affected  by  the  T,  transition)  while  the  ex¬ 
ponential  rise  in  p  agrees  with  the  gap  analysis  of 
Ref.  18.  It  is  easy  to  see  from  the  measured  tem¬ 
perature  behavior  of  RH/pfi,  that,  as  the  temper¬ 
ature  approaches  7\,  from  below,  either  p  rises 
rapidly  or  n  decreases  rapidly.  Since  the  conduc¬ 
tivity  in  both  directions  is  also  rising  in  this  tem¬ 
perature  range  the  first  case  is  the  obvious  choice. 
(We  have  repeated  the  computation  with  the  differ¬ 
ent  assumption  that  p  is  held  constant  between  36 
and  58  K  and  found  that  all  the  mobilities  turn 
around  and  increase  with  temperature  as  T,  is  ap¬ 
proached.  This  picture  is  certainly  harder  to 
justify.)  Thus  assumption  (d)  is  the  simplest  one 
that  provides  physically  reasonable  behavior  for 
the  mobilities  and  carrier  concentrations,  under 
the  four  constraints  provided  by  R„( 0),  p,,  p2, 
and  ptf. 

III.  RESULTS  AND  DISCUSSION 

Figure  1  shows  the  temperature  dependence  be¬ 
low  T.,  of  the  four  mobilities  obtained  by  the  solu¬ 
tion  of  Eqs.  (6)-(12).  (For  compactness  the  trans¬ 
verse  hole  mobility  ve  has  been  multiplied  by  100.) 
All  four  mobilities  show  power-law  behavior  in 
this  range  before  saturating  at  low  temperatures. 
The  longitudinal  mobilities  are  surprisingly  large 
at  liquid-helium  temperatures  (270000  cm"  /V  sec). 
However,  we  note  that  at  40  K  where  the  non-Ohmic 
effects  are  the  most  dramatic18* 19  p*  and  v,  are 
only  of  the  order  of  3  000  cm2/  V  sec.  At  4  K  the 
electron  pocket  anisotropy  is  15  compared  to  23.8 
for  the  hole  pocket.  The  carrier  concentration  n 
and  p  obtained  from  this  analysis  is  shown  in  Fig. 

2.  As  mentioned  above  both  n  and  p  are  tempera¬ 
ture  independent  below  40  K.  Above  40  K ,  p  rises 
rapidly,  reflecting  the  thermal  excitation  of  quasi¬ 
particles  across  a  diminishing  CDW  gap  on  the 
hole  surface.  The  solid  line  is  the  least-squares 
fit  which  provides  the  values  of  p  used  in  all  sub¬ 
sequent  computations.  The  broken  line  indicates 
the  assumed  temperature- independent  value  of  n 
above  40  K. 

To  check  the  validity  of  the  solution  we  have 
made  a  least-squares  fit  to  the  calculated  values 


in  Fig.  1  and  then  used  these  values  to  recompute 
p,,  p.,  Rh,  and  p v.  The  electron  concentration  is 
taken  to  be  fixed  at  1.088  ✓  1018  cm"1  and  />  is  taken 
from  Fig.  2.  Figure  3  shows  the  comparison  of 
the  computed  p,,  p,.  and  R  „( 0)  with  the  experi¬ 
mental  data  while  Fig.  4  shows  the  comparison  for 
p  M.  The  largest  discrepancy  occurs  for  li H( 0), 
especially  at  temperatures  below  12  K.  This  is  to 
be  expected  for  two  reasons.  Whereas  p, ,  pJt  and 
p*  measure  the  additive  effect  of  the  two  pockets 
the  Hall  constant  measures  the  competition  be¬ 
tween  the  mobilities  of  the  two  pockets.  Also  Rh 
cannot  be  measured  with  the  same  precision  as 
that  of  p,  and  p2.  At  low  temperatures  where  all 


FIG.  1 .  Temperature  dependence  of  the  electron  (solid 
circles)  and  hole  (open  circles)  mobilities  below  58  K, 
The  data  points  are  calculated  from  the  experimental 
data  using  Eqs.  (6)— (12) .  The  lines  are  a  least- squares 
tit  to  the  data.  In  all  cases  the  mobilities  obey  a  power- 
law  behavior  except  at  low  temperatures  where  they  sat¬ 
urate.  At  4  K  the  electron  mobility  ratio  equals  15, 
and  the  hole  mobility  ratio  is  23.8.  Note  that  the  trans¬ 
verse  hole  mobility  vc  has  been  multiplied  by  100. 
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FIG.  2.  Temperature  dependence  of  the  electron  (solid 
circles)  and  hole  (open  circles)  concentration.  The  data 
points  below  40  K  are  computed  from  Eqs.  (6)- (12)  of 
text,  and  the  error  bars  reflect  the  uncertainty  in  the 
measurements.  Above  40  K  the  broken  line  indicates  the 
assumed  temperature-independent  value  of  the  electron 
concentration  (1 .09  *  1018  cm"3).  The  open  circles  above 
40  K  are  the  computed  values  of  p  using  this  assumption. 
The  solid  line  is  the  least-squares  fit  to  the  computed 
data  on  p . 


the  mobilities  are  large  R„  becomes  very  sensi¬ 
tive  to  the  values  of  r,  t),  and  At  the  same  time, 
the  zero -field  value  of  R„  cannot  be  measured  ac¬ 
curately  due  to  the  narrow  range  of  fields  over 
which  R„  is  field  independent.  Thus  the  computed 
values  of  r,  tj,  and  £  below  12  K  are  accurate  to 
only  20*? .  We  have  found  also  that  the  values  of 
Rh(Q)  and  Rx  at  2  K  have  a  slight  variation  from 
sample  to  sample.30 

The  computed  field  dependence  of  R  „  is  shown 
in  Figs.  5  (low  fields)  and  6  (high  fields).  (The 
experimental  data  in  Fig.  6  are  from  FPC.)  While 
the  agreement  in  Fig.  5  is  excellent  except  at  very 
low  temperatures,  in  Fig.  6  we  have  had  to  adjust 
the  temperature  of  the  sample  to  improve  the  fit 
at  high  fields.  At  low  fields  (~20  kG)  the  disagree¬ 
ment  is  larger.  The  value  of  RK(0 )  (the  slope  in 
Fig.  6)  at  30  K  is  approximately  equal  and  opposite 
in  sign  to  Rx  at  4.2  K  in  FPC’s  data  whereas  in 
the  data  used  for  this  analysis  the  same  quantity 
is  only  40%.  This  accounts  for  the  large  discre¬ 
pancy  in  Fig.  6  at  low  fields.  The  disagreement 
may  be  sample  dependent. 

We  now  turn  to  the  SdH  data  and  see  to  what  ex¬ 
tent  the  model  described  here  is  supported  by  the 
FS  mapping  obtained  by  means  of  the  Shubnlkov 
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FIG.  3.  Comparison  of  the  experimental  data  for  the 
resistivities  pb  and  pc  (circles),  and  the  zero-field  Hall 
constant  firt(0)  (triangles)  with  the  two-band  model  cal¬ 
culation.  The  lines  represent  the  values  of  RH( 0),  and 
pb  and  pc  calculated  using  Eqs.  (9)  and  (12).  The  values 
of  the  mobilities  used  are  given  by  the  full  lines  in  Fig. 

1.  The  electron  concentration  is  assumed  fixed  at 
1.09  xlO18  cm"3  and  the  hole  concentration  is  given  by  the 
full  line  in  Fig.  2.  Values  of  RH{ 0)  at  low  temperatures 
appear  to  be  sample  dependent.  (The  open  and  solid  tri¬ 
angles  are  data  from  two  samples  normalized  to  the  same 
value  at  infinite  field.) 


oscillations.  The  most  extensive  measurements 
on  NbSe3  are  those  of  MB  whose  measurements  in¬ 
dicate  the  presence  of  two  branches  and  possibly 
a  third.  Along  directions  of  the  magnetic  field 
where  the  oscillations  are  strong  the  data  of  MB 
agree  with  the  preliminary  data  of  M  and  FPC. 

The  first  branch  which  is  the  most  well  defined  in 
MB’s  data  corresponds  to  an  approximately  ellip¬ 
soidal  surface  with  its  principal  axes  along  a*,  b. 
and  c.  The  shortest  dimension  is  along  b  and  the 
longest  is  along  c.  If  the  semiaxis  dimensions  are 
k,*,  kb,  and  kc  the  volume  of  the  ellipsoid  is 
V=  {Scx  (S,/4ne),/2  where  S,=  4nkbkc,  Sc  =  4»fe,«fe>. 
and  e  -kb/kc.  From  the  SdH  data  of  M,  FPC.  and 
MB,  Se  =  0.3  MG  (2.86  x  10ls  cm"2),  S,~-  1.05  MG 
(1.00  x  10“  cm"2).  The  carrier  density  within  this 
FS  pocket  is  thus  equal  to  u  =  2.17x  10I7/VT  cm"3. 
FPC  find  that  e  =*  £  gives  the  best  fit  to  the  angular 
dependence  of  the  frequency  as  the  magnetic  field 
is  rotated  towards  the  b  axis.  MB’s  data  favor  a 
smaller  e.  In  any  case  the  carrier  concentration 
is  of  the  order  of  6  x  1017cm"3  which  is  in  order- 
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FIG.  4.  Comparison  of  the  experimental  data  for  the 
magnetoreststance  mobility  with  the  two-band  model  cal¬ 
culation.  The  full  line  is  calculated  from  Eq .  (11)  with 
the  model  parameters  given  by  the  full  lines  in  Figs .  I 
and  2.  Data  from  three  samples  (Ref.  1)  are  displayed. 


FIG.  5.  Field  dependence  below  14  kG  of  the  Hall  con¬ 
stant  at  various  temperatures  below  40  K.  The  full  lines 
are  the  experimental  curves  from  Ref.  6  and  the  broken 
lines  are  calculated  from  Eqs.  (6)- (8)  using  model  para¬ 
meters  given  by  the  full  lines  in  Figs.  1  and  2.  The  com¬ 
puted  Rh  has  been  adjusted  by  a  scale  factor  for  com¬ 
parison  with  the  measurements  (see  Ref.  20).  Agreement 
is  excellent  except  at  low  temperatures  where  R H  is  very 
sensitive  to  model  parameters. 


of-magnitude  agreement  with  the  electron  concen¬ 
tration  computed  here  (1.09  x  lola  cm"3).  Hence 
the  ellipsoidal  surface  measured  in  the  SdH  data 
is  probably  the  electron  pocket.  However,  there 
are  several  difficulties  in  a  careful  comparison  of 
the  detailed  SdH  measurements  of  MB  and  FPC 
with  our  conclusions.  First,  the  value  of  s  for  e 
would  imply  a  mobility  ratio  of  64  (assuming  a 
parabolic  band)  in  the  electron  pocket.  This  is 
four  times  larger  than  the  computed  anisotropy  in 
Fig.  1.  This  disagreement  is  probably  due  to  the 
fact  that  this  FS  pocket  resembles  more  a  brief¬ 
case  than  an  ellipsoid.  MB's  data  agree  with  a  cy¬ 
lindrical  trajectory  (t  =  0)  with  the  magnetic  filed 
tilted  up  to  an  angle  of  70  from  the  a*c  plane. 

Secondly,  MB’s  data  show  other  branches  of  the 
same  order  of  magnitude  in  frequency  which  merge 
abruptly  with  the  main  branch.  Over  a  certain 
angular  range  when  the  magnetic  field  is  in  the 
a*c  plane  the  main  branch  fades  dramatically, 
leaving  a  ghost  image  of  its  original  sec  8  trajec¬ 
tory,  while  continuing  with  undiminished  intensity 
into  the  new  branch.  These  complications  suggest 
that  a  simply  connected  ellipsoidal  pocket  is  prob- 


FIG.  6.  Field  dependence  at  high  fields  of  the  Hall  re¬ 
sistance  at  various  temperatures.  The  full  lines  are 
from  Fleming,  Polo,  and  Coleman  (Ref.  8).  The  broken 
lines  are  calculated  from  the  two-band  model  using  para¬ 
meters  given  by  the  solid  lines  in  Figs.  1  and  2.  The 
calculated  Hall  resistance  has  been  normalized  to  agree 
with  the  experimental  data  at  4.2  K.  The  temperature  of 
the  sample  has  also  been  adjusted  to  give  a  better  fit  at 
high  fields.  There  is  a  disagreement  in  the  measured 
Hall  polarity  between  Ref.  6  and  Fleming  et  al.  However, 
for  the  purpose  of  comparison  we  have  changed  the  sign 
of  the  carriers  in  this  figure. 
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ably  too  naive  a  picture.  A  possible  explanation  of 
these  intriguing  patterns  is  the  intersection  of  the 
original  picket  by  higher-order  superlattice  Bragg 
planes.  Indeed  it  may  be  shown  1  that  several  se¬ 
cond-order  images  of  the  /'.-CDW  Bragg  plane 
oriented  at  45  to  the  A- axis  intersect  the  center 
of  the  Brillouin  zone  where  the  electron  pocket 
would  presumably  be  located.  The  faint  ghost 
spots  may  be  the  result  of  magnetic  breakdown  ' 
across  these  small  gaps  in  analogy  with  the  case23 
in  Cr.  while  the  abrupt  appearance  of  a  new  branch 
may  be  the  deflection  of  the  electron  trajectory  in¬ 
to  a  new  external  orbit  by  the  second-order  gap. 

The  separation  of  the  contributions  of  the  two 
bands  to  the  transport  properties  enables  us  to 
justify  the  argument  used  in  the  preceding  paper 
with  respect  to  RK{pj>e)~'-  There  it  was  argued 
that  the  quantity  A  =RH(0 )  lpbpe)  had  a  temperature 
dependence  which  was  dominated  by  the  lifetime 
of  the  carriers.  Dividing  out  the  power-law  be¬ 
havior  of  the  lifetimes  enables  us  to  observe  the 
relative  change  in  carrier  concentration.  From 
Fig.  1  the  temperature  dependence  of  the  mobilities 
above  20  K  is  T*1,2(m,),  7"2-8(p,),  r'3-5(v,),  and 
T*'‘,n(v2).  Thus  the  quantity  A  (T)  =  -(np.,Pj -pvtv2) 
has  the  temperature  dependence  -(nT*8 •“  - c/>7'*5-5) 
where  c  is  a  constant.  Within  the  calculated  pa¬ 
rameters  of  the  model  the  second  term  is  21%  of 
the  first  at  50  K  and  increases  to  61%  at  58  K. 

Thus  the  first  term  dominates  the  behavior  of  A 
in  this  temperature  range.  We  would  then  expect 
the  quantity 

A(T  T,)3-875  ocr0-1  (n -cfiT-1-4) , 

to  measure  the  relative  change  in  carrier  concen¬ 
trations.  The  residual  power  law  dependence 
T*1,8  multiplying/)  unfortunately  distorts  the  actu¬ 
al  temperature  dependence  of  />.  However  this  is 
not  serious  near  the  transition  because  of  the  ex¬ 
ponential  increase  in  p.  (A  comparison  of  the  data 
in  Fig.  3  of  the  preceding  paper  with  that  in  Fig. 

2  here  shows  the  slight  difference  between  the 
temperature  behavior  of  A(T/Ti)3‘37s  and  p.) 

IV.  CONCLUSION 

The  galvanomagnetic,  conductivity  and  Shubnikov 
data  obtained  to  date  on  NbSe3  provide  a  rather 
complete  picture  of  the  low-temperature  phase. 

To  sort  out  the  different  components  which  contri¬ 
bute  to  the  observed  quantities  we  have  adopted 
the  simplest  two-band  model  and  used  the  experi¬ 
mental  data  to  calculate  the  six  unknowns  in  the 
model.  Enough  parameters  are  fixed  by  the  mea¬ 
surements  to  leave  no  free  parameters  below  40 
K.  Above  40  K  we  have  had  to  assume  that  the 
electron  concentration  remains  temperature  inde¬ 


pendent  and  shown  that  this  is  borne  out  self-ron- 
sistently  by  the  results.  The  picture  that  we  ob¬ 
tain  is  that  just  below  7',  a  CDW  gap  occurs  on  the 
hole  surface.  As  the  hole  carriers  are  frozen  out 
across  the  developing  CDW  gap.  the  total  hole  con¬ 
centration  decreases  by  a  factor  of  4.  This  re¬ 
sults  in  only  a  factor  of  2  rise  in  the  longitudinal 
resistivity  and  an  even  smaller  increase  in  the 
transverse  resistivity  because  the  electron  mobili¬ 
ties  are  higher  at  this  temperature.  Below  40  K 
the  carrier  freeze  out  is  complete,  leaving  6.0 
xtO18  cm*3  holes  and  1.09  x  10la  cm*3  electrons. 

The  remaining  holes  are  from  portions  of  the 
same  FS  which  are  not  destroyed  by  the  gap  or 
may  possibly  belong  to  an  entirely  separate  poc¬ 
ket.  Below  36  K  the  Hall  constant  at  zero  field  in¬ 
creases  to  more  positive  values  and  changes  sign 
at  15  K  because  of  the  rapidly  increasing  hole 
mobilities.  At  very  low  temperatures  the  longi¬ 
tudinal  mobilities  become  so  large  ('250  000  cm' 

V  sec)  that  the  Hall  signal  is  rapidly  driven  to  its 
infinite  field  value  with  the  relatively  modest  field 
of  10  kG.  In  contrast  to  the  Hall  and  resistivity 
data  which  show  peak  structures  the  transverse 
magnetoresistance  obeys  a  smooth  power-law  be¬ 
havior.  This  model  also  accounts  for  the  strong 
field  dependence  of  fi„  (Fig.  6).  The  presence  of 
the  electron  pocket  is  strongly  indicated  by  the 
SdH  data  although  the  possibility  of  superlattice 
Bragg  planes  intersecting  the  electron  FS  may 
complicate  the  full  interpretation  of  the  SdH  data. 
Evidence  for  the  hole  pocket  is  not  apparent  in  the 
SdH  data.  Based  on  the  parameters  of  this  model, 
the  resistivities.  Hall  constant,  and  magnetore¬ 
sistance  are  calculated  and  good  agreement  is  ob¬ 
tained  with  the  experimental  data. 

The  number  of  carriers  participating  in  the  T, 
transition  is  approximately  2x  1019  cm*3  (16  times 
the  electron  concentration).  This  number,  based 
on  Fig.  2  is  not  meaningful  if  the  gapping  portion 
of  the  hole  surface  is  not  closed  (since  Eq.  (5)  is 
then  no  longer  valid].  In  fact  the  value  of  the 
transverse  components  of  the  spanning  vector  at 
T,  (0.5a ',0.5c*)28  suggests  that  the  FS  that  nests 
may  be  a  cosine  band25  which  intersects  the  Bril¬ 
louin  zone  face.  In  this  case  Eq.  (5)  is  clearly  in¬ 
adequate  to  describe  contributions  from  all  three 
sources.  A  third  matrix  obtained  from  a  current 
term  such  as 

,  2  2e2r  f  dSF  , 

3 = (T^ny  (2ir)3  J„p  i n ' AVk 

(where  the  integral  is  over  the  gapped  portion  of 
the  FS)  should  be  included.  In  spite  of  this  quali¬ 
fication  it  is  remarkable  that  the  temperature  de¬ 
pendence  of  p  computed  from  Eq.  (5)  shows  expli¬ 
citly  the  exponential  rise  in  quasiparticles  as  the 
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T,  transition  is  approached.  The  number  2  x  1019 
cm*5  should  serve  adequately  as  an  order-of- 
magnitude  estimate  of  the  FS  destroyed  by  the  gap 
at  T 
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APPENDIX 

Equations  (6)  and  (9)- (12)  may  be  reduced  to  a 
single  equation  involving  only  r  as  the  unknown  as 
follows.  First,  five  of  the  equations  are  rewritten 
in  terms  of  the  dimensionless  parameters  [see  Eq. 
(15)1  r.  5.  and  q,  as 

1/p,  =/>iq(ri/  +  l) ,  l/p.^pv.JrZ  +  l),  R<*  =  l/p(l  -  r) 

(Al) 


Mjr/lPiPj1  /2=- /><',«'-  frnd  +  D  ,  (A2) 

(R«  -Ru)/PiP2  =  |rpt/|v2/(l  -r)J({  +  l)(q  +  l).  (A3) 

By  forming  the  dimensionless  constants  C,  and  C2 
we  eliminate/*,  ix, ,  and  im  Thus  we  have 


P  -Hj/(PiP2)i/2  (1-r)  Sri 

(A4) 

‘  (R»-flo)  fT  (ltij)’ 

p  -Pm(PiPJi/2  ^<l-r)({  +  l) 

(A5) 

(r*  +  l)(rq  +  l)  " 

With  r  as  the  independent  variable  q  and  (;  can  be 
expressed  in  terms  of  r  through  Eqs.  (A4)  and 
(A5).  The  mobilities  v,  and  ic  can  then  be  found 
from  Eq.  (Al)  in  terms  of  r.  Finally  Rh(H),  R* 
can  be  computed  as  a  function  of  H  by  Eq.  (6)  with 
r  as  an  adjustable  parameter.  The  value  of  r 
which  provides  the  best  fit  of  RH/RX  to  the  experi¬ 
mental  data  (Fig.  5)  is  taken  to  be  the  correct  so¬ 
lution. 
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Studies  on  die  effect  of  charged  (Ti)  and  unchar  aged  (Ta)  Impurities  on  die  non- Ohmic  be¬ 
havior  of  and  microwave  absorption  in  NbSe3  show  that  (1)  the  characteristic  field  E  0  increas¬ 
es  aa  c2  (c  =  Ta  concentration) ,  in  agreement  with  a  recent  calculation  by  Lee  and  Rice;  (2)  the 
microwave  absorption  is  suppressed  by  impurities;  and  (3)  charged  impurities  are  strikingly 
more  effective  in  supressing  the  anomalous  properties.  These  conclusions  provide  very 
strong  evidence  for  Frohllch  sliding-mode  conductivity  in  the  nominally  pure  compound. 

The  central  question1  in  the  study  of  charge-  much  Interest  in  this  regard.  Recent  electron4 
density-wave  (CDW)  or  Pelerls  systems  may  and  x-ray5  diffraction  measurements  have  con- 

very  well  be  whether  collective  modes  are  sig-  firmed  superlattice  formation  at  both  Tl  and  T*, 

niflcant  In  the  conductivity  in  the  CDW  phase.  and  it  now  appears  certain  that  these  anomalous 

The  striking,  anomalous  transport  properties*'1  properties  are  due  to  excitations  in  the  CDW  con- 

discovered  in  NbSe,  below  the  two  phase  transi-  densate.  Below  Tl  and  T,  NbSe,  is  non-Ohmic. 

tions  at  T,  (142  K)  and  Tt  (58  K)  have  aroused  It  also  shows  a  strong  absorption  mode  at  micro- 
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wave  frequencies.3  The  small  dc  fields  £  (50  mV/ 
cm  near  T2;  1  v/cm  near  T,)  required  to  sup¬ 
press  the  anomalies  as  well  as  the  small  photon 
energies  (4X10'5  eV)  absorbed  in  the  microwave 
measurements  have  defied  all  theoretical  analy¬ 
ses3,8  unless  it  is  postulated  that  a  macroscopic 
number  of  charges  couples  coherently  to  the  ap¬ 
plied  fields. 

As  first  discussed  by  Frohlich,7  the  condensa¬ 
tion  of  an  incommensurate  CDW  breaks  the  trans¬ 
lational  symmetry.  The  zero-energy  excitation 
associated  with  this  broken  symmetry  is  rigid 
motion  of  the  CDW.  In  a  real  crystal,  however, 
defects  will  pin8  the  CDW.  The  application  of  a 
sufficiently  high  electric  field  will  depin  the  con¬ 
densate  leading  to  non-Ohmic  conductivity.  Also 
weak  ac  fields  will  cause  segments  of  the  CDW  to 
oscillate  about  the  pinning  centers,  leading  to  an 
absorption  mode  near  the  resonance  frequency. 
This  model  predicts  that  increasing  the  impurity 
concentration  c  will  increase  the  characteristic 
depinning  field  E0,  as  well  as  the  pinning  frequen¬ 
cy  (by  stiffening  the  spring  constant).  Charged 
impurities  are  also  expected  to  pin  the  conden¬ 
sate  more  strongly  than  uncharged  ones.  Recent¬ 
ly  Lee  and  Rice  (LR)9  have  applied  this  model  in 
detail  to  NbSea,  and  predict  that  E0  increases  as 
c 2  in  the  weak-pinning  (uncharged  impurities) 
case.  The  first  test  of  these  ideas  is  the  study 
of  the  effects  of  impurities  on  the  anomalous 
properties.  Our  data  confirm  all  the  qualitative 
arguments  described  above,  and  provide  experi¬ 
mental  verification  of  the  c2  dependence  of  E0  as 
well. 

Impurities  were  introduced  by  first  forming 
the  master  alloy  Nbj .%MX  (M  =  Ta  or  Ti).  Sever¬ 
al  remelting  and  cutting  cycles  were  performed 
to  achieve  good  homogeneity.  The  alloyed  com¬ 
pound  Nbj.,  M, Sea  was  then  grown  as  described 
by  Meerschaut  and  Rouxel.10  The  starting  Nb 
metal  had  a  nominal  concentration  of  180  ppm 
Ta.  We  have  used  the  residual  resistivity  ratio 
(RRR,  denoted  also  as  (ft)  to  monitor  sample  pur¬ 
ity.  [A  plot  of  (ft  versus  the  nominal  concentra¬ 
tion  of  Ta  impurities  cl  shows  that  1/<R  obeys  the 
equation  1/(R  =  (5.2±  0.2)*  10'3+  (19±  2)Cj(Cj<  3 
*10‘3).  In  general11  we  may  write  1/<R <*c,  where 
c  (=c0+  r,)  is  the  total  concentration  of  impuri¬ 
ties  and  defects,  and  c0  the  concentration  of  de¬ 
fects  excluding  Ta  impurities.]  Details  of  the 
non-Ohmic3  and  microwave  experiments12  have 
previously  appeared. 

Below  both  7\  and  T2  the  dc  conductivity  obeys 
the  empirical  behavior 3  o(E,T)  =oa(T)  +  cb (T) 


*exp [- E0(T)/e],  This  equation  is  valid  in  all 
samples  studied.  The  characteristic  field  E0  is 
strongly  temperature  dependent  (Fig.  1,  inset), 
attaining  a  minimum  value  £0min  a  few  degrees 
below  T2.  This  minimum  provides  a  convenient 
and  unambiguous  quantity  for  comparison  from 
sample  to  sample.13  For  each  sample  £0mln  is 
determined  by  measuring  E0  at  four  tempera¬ 
tures  bracketing  the  minimum.  We  have  plotted 
in  Fig.  1,  £0I,1,n  vs  (ft"1  for  each  Ta-doped  sam¬ 
ple.  Since  (ft1  is  linear  in  c,  £0min  varies  as  c2 
(solid  line).  A  least-squares  fit  to  the  data  gives 
£0m,n  =  44?<R-1-87  V/cm.  Also  aparent  in  Fig.  1  is 
that  in  the  low-c  limit  the  data  for  £0  show  no 
sign  of  saturation.  Annealing  of  the  purest  sam¬ 
ples  may  produce  a  further  decrease  in  Ev  (Our 
results  indicate  that  a  RRR  much  larger  than  200 
cannot  be  obtained  by  further  removal  of  Ta  im¬ 
purities.)  It  should  be  very  interesting  to  see 
how  low  a  depinning  field  one  can  attain  in  a  real 
solid.  The  effect  of  charged  impurities  (Ti)  has 
also  been  studied.  A  concentration  of  1000  ppm  Ti 


10* 


FIG.  1.  The  minimum  characteristic  field  E  0  for 
each  sample  vs  the  reciprocal  of  Its  residual  resistivity 
ratio  (RRR).  The  nominal  Ta  impurity  concentration 
c,  is  1900  ppm  (open  triangles),  1300  ppm  (closed  cir¬ 
cles),  300  ppm  (open  circles),  and  180  ppm  (closed  tri¬ 
angles).  Values  for  a  sample  quoted  in  Ref.  5  are 
shown  as  plus  sign.  The  solid  line  has  slope  2.0.  Since 
tfC 1  is  linear  In  c,  K0rr,|n  varies  as  c2,  in  agreement 
with  calculations  in  Ref.  9.  The  inset  shows  the  tem¬ 
perature  dependence  of  F0  below  T}  (58  K>. 
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nails  down  the  condensate  to  the  extent  that  no 
Ohmic  breakdown  could  be  observed  (although  a 
lower  bound  of  6  V/cm  may  be  inferred.)  This 
is  to  be  compared  to  a  typical  value  of  300  mV/ 
cm  for  Ta  impurities  in  the  same  concentration. 
The  variation  of  E0  with  the  RRR  in  two  samples 
has  also  been  reported  by  Fleming,  Moncton, 
and  McWhan.5 

The  effects  of  both  Ta  and  Ti  impurities  on  the 
microwave  absorption  mode  are  shown  in  Fig.  2, 
The  absorption  strength  is  proportional  to  the 
difference  between  the  dc  and  ac  data  for  each 
sample.  As  the  Ta  concentration  increases  (mov¬ 
ing  up  in  the  figure)  one  sees  the  suppression  of 
the  anomalous  absorption  mode.  The  data  are 
consistent  with  the  pinning  frequency  changing 
monotonically  with  the  increasing  impurity  con¬ 
centration.  Again  it  is  also  clear  that  Ti  impur¬ 
ities  have  a  more  drastic  effect  on  the  micro- 
wave  absorption  than  a  similar  concentration  of 
Ta  (samples  B  and  D ).  No  dc  Ohmic  breakdown 
was  observed  (up  to  a  field  of  6  V/cm)  in  sam¬ 
ples  A  (50000  ppm  Ta),  B  (1000  ppm  Ti),  and  C 


(5000  ppm  Ta). 

Since  claims  on  anomalous  effects  which  are 
highly  sensitive  to  sample  purity  have  been  con¬ 
troversial  in  the  past,  it  appears  essential  that 
an  independent  monitor  of  sample  purity  be  used. 
For  each  sample  studied  we  have  measured  the 
width  of  the  critical  divergence  of  dlnp/dT  (p  is 
the  low-field  resistivity)  at  T,  as  well  as  Tt  it¬ 
self.  In  Fig.  3  data  from  five  samples  show  the 
monotonic  correlation  between  the  RRR  on  the 
one  hand,  and  the  width  and  T2  on  the  other.  The 
latter  thus  provide  an  independent  measure  of 
the  sample  purity. 

The  failure  of  the  usual  Drude  model  (to  explain 
the  anomalous  transport  properties  in  the  CDW 
phase)  can  be  quantitatively  discussed  using  ener¬ 
gy  considerations.  Equating  the  minimum  field 
energy  (£0~  40  mV/cm)  picked  up  by  an  electron 
to  the  thermal  fluctuation  energy  at  50  K  we  find 
that  the  mean  free  path  must  be  of  the  order  of 
0.1  cm — a  truly  macroscopic  distance.  (This  is 
indefensible,  particularly  when  a  two-band  analy¬ 
sis14  of  NbSes  based  on  all  available  galvanomag- 


FIG.  2.  The  temperature  dependence  of  the  resistivi¬ 
ty  at  9.8  GHz  for  various  samples  of  Nb,.IAfISe3  (Af  =  Ta 
or  Ti) .  The  solid  line  is  the  dc  restsitivty  for  each 
sample.  As  discussed  in  the  text  the  absorption  mode 
is  suppressed  with  increasing  impurity  concentration 
(moving  up  in  figure) .  The  nominal  concentratlons/RRR 
are  (50000  ppm  Ta)/1.14  (curve  A),  (1000  ppm  T0/1.5S 
(curveJB),  (5000  ppm  Ta)/2.99  (curve  C),  (1900  ppm 
Ta)/  18.2  (curve  D),  and  (300  ppm  Ta)/115  (curve  F). 

Ti  impurities  are  more  effective  In  suppressing  the  ab¬ 
sorption  mode  than  Ta.  For  clarity  the  data  for  curves 
C,  D,  and  E  have  been  displaced  downwards. 


T  ( K ) 

3.  Critical  divergence  of  dlnp/dT  at  T ,  for  five 
samples  showing  the  decrease  in  T2  and  the  broadening 
of  the  transition  with  increasing  c,.  The  Ta  concentra- 
tions/RRR  are  (5000  ppm)/3  (circles  within  circles), 
(1900  ppm)/25  (Inverse  triangles),  (1300  ppm)/42  (open 
circles),  (300  ppm)/91  (closed  circles),  and  (180  ppm)/ 
111  (triangles). 
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netic  data  shows  that  the  uncondensed  carriers 
have  a  rather  ordinary  mobility  of  3000  cm*/V  s 
at  SO  K  comparable  to  that  in  Ge  at  300  K.)  Simi¬ 
lar  attempts  using  single -particle  theory  to  ex¬ 
plain  the  extra  microwave  absorption  (for  exam¬ 
ple,  by  the  photogeneration  of  carriers  across 
minigaps  or  out  of  traps)  lead  to  gap' or  trap  sizes 
which  are  (l/100)hT,  which  is  clearly  untenable. 
Models  based  on  the  partial  destruction  of  the 
CDW  state  by  the  applied  E  field  can  also  be 
ruled  out  experimentally.  The  transition  temper¬ 
atures  7\  and  Xa  are  not  affected  by  the  highest 
applied  fields.2  More  significantly,  direct  obser¬ 
vation  of  the  super  lattice  spot  intensity  by  Flem¬ 
ing,  Moncton,  and  McWhan5  shows  no  change  in 
intensity  under  fields  strong  enough  to  suppress 
the  T2  anomaly  by  half. 

In  Ref.  9,  one  of  the  mechanisms  for  non-Oh- 
mic  behavior  discussed  by  LR  is  the  depinning  of 
a  rigid  CDW  domain  of  length  L  by  an  applied 
field.  Equating  the  energy  gained  from  the  field 
to  the  pinning  energy  of  a  domain,  LR  get  the  es¬ 
timate  (in  the  weak -pinning  case) 

„  1  1  (P.V)4^ 

0=  epP'ff2ir/Q  1*1*  /0SM,I .?  * 

Here  ip  is  the  order  parameter,  p^V  the  impurity 
potential,  £,  the  coherence  length  along  i,  and  Q 
the  spanning  vector.  The  c*  dependence  is  veri¬ 
fied  by  our  data.  Using  reasonable  numbers,  LR 
calculate  E0  to  be  10  mV/cm  for  c  *  10' 9,  which 
corresponds  to  (ft'1-  5.2*10'3.  This  is  in  order- 
of-magnitude  agreement  with  the  data  in  Fig.  1. 
However,  neither  the  nonmonotonic  temperature 
dependence  of  E0  nor  the  field  dependence  of 
o(E,  T)  is  accounted  for  in  LR’s  zero-tempera¬ 
ture  calculation.  These  remain  serious  problems 
in  the  model. 

In  conclusion  our  results  show  that  the  non-Oh- 
mic  behavior  at  dc  and  the  strong  absorption 
mode  at  9.8  GHz  are  intrinsic  to  the  nominally 
pure  compound,  and  easily  suppressed  by  weak 
concentrations  of  impurities.  Charged  impurities 
(Ti)  are  strikingly  more  effective  in  suppressing 
these  anomalous  properties  than  uncharged  ones 
(Ta).  The  measured  c*  dependence  of  £0n,ln 
agrees  with  LR's  calculations  based  on  the  de¬ 
pinning  of  a  CDW  domain.  These  results  provide 
the  strongest  evidence  to  date  for  the  FrShlich 


mechanism  in  a  CDW  metal. 

We  have  benefitted  from  discussions  with  P.  A. 
Lee,  T.  M.  Rice,  S.  E.  Trullinger,  and  K.  Maki. 
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fice  of  Naval  Research  (N00014-77-C0473)  and 
in  part  by  the  National  Aeronautics  and  Space  Ad¬ 
ministration  (NAS7-100). 
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Electric- Field -Dependent  Thermopower  of  NbSe3 
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We  have  measured  the  thermoelectric  power  of  NbSe3  in  the  temperature  range,  10-70 
K,  where  charge-density-wave  formation  and  strongly  nonlinear  conductivity  has  been  re¬ 
ported.  We  report  the  first  observation  of  an  electrlc-fleld-dependent  thermopower. 

Our  results  for  NbSe,  show  that  In  the  presence  of  an  electric  field  there  is  an  additional 
negative  contribution  to  the  thermopower.  Possible  mechanisms  which  may  account  for 
this  are  discussed. 

There  has  been  a  great  deal  of  recent  interest 
in  the  non-Ohmic  behavior  observed  in  systems 
which  undergo  charge-density-wave  (CDW)  in¬ 
stabilities.1'1  Upon  application  of  an  electric 
field,  the  electrical  conductivity  is  found  to  rise 
dramatically  when  a  certain  threshold  is  sur¬ 
passed.  This  effect  is  often  found  in  quasi -one  - 
dimensional  conductors  with  tetrathiafulvalene 
tetracyanoquinodlmethane1  (TTF-TCNQ)  and 
NbSe3  (Ref.  1)  as  primary  examples.  A  variety 
of  models  have  been  proposed  to  explain  the  non¬ 
linear  characteristics.1*9  Since  the  models  pre¬ 
dict  different  types  of  charge  carriers,  thermo¬ 
electric  power  (TEP)  measurements  can  contrib¬ 
ute  to  the  understanding  of  the  nature  of  the  car¬ 
riers. 

In  this  Letter  we  report  on  TEP  measurements 
performed  on  NbSe,  in  the  presence  of  an  electric 

field.  This  is  the  first  time  that  thermopower  1  p(E)+S{E)VT, 


has  been  used  to  probe  nonlinear  effects  in  any 
material  and  consequently  the  first  time  that  an 
electric-field-dependent  TEP  has  been  observed. 
Below  the  59-K  CDW  transition  the  measure¬ 
ments  show  an  increase  in  the  magnitude  of  the 
TEP  as  electric  field  is  increased,  the  limiting 
value  approaching  the  TEP  observed  above  the 
59-K  transition.  Since  a  collective  mode  (i.e., 
sliding  CDW)  by  itself  would  produce  a  reduced 
(or  zero)  TEP,  the  present  work  implies  that 
some  other  mechanism  is  responsible  for  the  ob¬ 
served  increase  in  TEP  and  conductivity. 

In  the  experiments  we  have  measured  both  the 
electrical  resistance  and  the  TEP  of  NbSes  as  a 
function  of  both  temperature  and  electric  field. 
The  experimental  situation  is  governed  by  the 
transport  equation 
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where  ®  is  the  electric  field,  3  is  the  electric 
current  density,  VT  is  the  temperature  gradient, 
and  for  NbSe3  both  p(£')  and  S(E)  are  dependent  on 
the  electric  field  at  temperatures  below  142  K.1 
In  order  to  measure  the  TEP  in  a  high  electric 
field  we  have  to  superimpose  a  heat  current  onto 
an  electric  current  passing  through  the  sample 
and  measure  the  incremental  voltage  caused  by 
the  subsequent  temperature  difference.  Use  of  a 
steady  current  3  to  produce  the  electric  field 
proved  to  be  difficult,  as  a  large  dc  voltage  is 
produced  if  the  mean  temperature  of  the  sample 
changes.  This  is  due  to  the  strong  temperature 
dependence  of  the  resistance  in  NbSe3.  On  apply¬ 
ing  heat  to  produce  a  temperature  gradient  in  the 
sample  even  small  changes  in  temperature  cause 
this  voltage  to  be  much  larger  than  the  thermo¬ 
electric  voltage.  To  eliminate  this  problem  a 
low-frequency  ac  current  was  applied  to  the  sam¬ 
ple  via  a  dc  blocking  capacitor.  The  dc  voltages 
across  the  sample  were  measured  after  the  ac 
background  was  eliminated  by  passing  the  signal 
through  a  PAR-113  amplifier  set  at  dc  with  a 
suitable  rolloff  filter.  Use  of  ac  electric  fields 
give  similar  experimental  conditions  to  dc  as  the 
same  breakdown  effects  are  observed.  However, 
measured  values  of  p  (E)  would  differ  slightly  as 
the  ac  measurement  would  take  a  different  aver¬ 
age  over  the  nonlinear  region.  Any  incremental 
dc  effect  superimposed  on  the  ac  field  would  be 
made  in  the  nonlinear  region  although  the  actual 
value  of  the  effective  applied  electric  field  would 
depend  on  the  characteristics  of  the  current-ver- 
sus-voltage  (/-  V)  curve.  These  conditions  were 
checked  by  superimposing  a  small  dc  current  (10 
pA)  upon  the  applied  ac  current  and  the  observed 
dc  voltage  was  reduced  above  the  threshold  cur¬ 
rent  as  expected.  Thus  for  a  thermoelectric 
measurement  we  can  be  confident  that  we  are  in 
the  nonlinear  region  and  the  observed  effects  con¬ 
firm  this. 


Cu  cold  plot* 
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NbSe^  sample 
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FIG.  1.  The  experimental  arrangement  used. 


The  experimental  configuration  is  shown  in 
Fig.  1.  A  needle  of  NbSe3  was  thermally  anchored 
to  two  sapphire  blocks  with  indium  solder  and  sil¬ 
ver  paint  which  acted  as  current  and  voltage  con¬ 
tacts,  respectively.  The  thermoelectric  measure¬ 
ments  were  made  differentially  by  switching  heat 
from  one  sapphire  block  to  the  other  alternately 
in  the  presence  of  a  constant  ac  current.  This 
enabled  the  thermoelectric  voltage  to  be  distin¬ 
guished  from  any  dc  offsets.  The  resulting  vol¬ 
tages  and  temperature  differences  were  meas¬ 
ured  using  conventional  nanovoltmeters  and  ther¬ 
mocouples.  The  low  temperatures  were  achieved 
by  attaching  the  assembly  to  an  Air  Products  Dix- 
plex  refrigeration  system. 

Initially,  measurements  of  the  non-Ohmic  re¬ 
sistance  were  made  in  the  region  of  the  low-tem¬ 
perature  CDW  transition.  This  was  documented 
in  three  forms  (i)  by  measuring  the  /-  V  charac¬ 
teristic  at  various  temperatures;  (ii)  by  meas¬ 
uring  the  dc  resistance  using  a  small  constant 
current  (10  pA)  at  various  values  of  the  applied 
ac  current;  and  (iii)  by  sweeping  the  temperature 
with  various  ac  sample  currents  and  monitoring 
the  ac  voltage.  The  resulting  R-vb-T  data  [using 
method  (iii)]  are  shown  in  Fig.  2(a).  As  can  be 


FIG.  2.  (a)  The  electrical  resistance  and  (b)  the 
thermoelectric  power  of  NbSe3  as  a  function  of  tem¬ 
perature  between  10  and  80  K  at  several  values  of  die 
applied  ac  current. 
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FIG.  3.  The  electrical  conductivity  and  thermoelec¬ 
tric  power  of  NbSe,  as  a  function  of  applied  current  at 
T  =47  K,  normalized  to  values  In  die  limit  of  zero  cur¬ 
rent. 


seen  the  electric  field  breakdown  is  similar  to 
that  observed  previously,1  and  for  the  current 
levels  used  sample  heating  is  only  evident  in  the 
highest  current  level  used  (i.e.,  10  mA). 

The  TEP  data  in  this  non-Ohmic  region  is  shown 
in  Fig.  2(b).  The  zero-applied-electric -field  data 
shown  as  the  closed  circles  are  similar  to  those 
obtained  by  others.*  As  the  sample  current  (i.e., 
the  applied  electric  field)  is  increased,  the  TEP 
increases  for  electric  fields  above  some  thresh¬ 
old  value  in  similar  fashion  to  the  electrical  con¬ 
ductance.  This  is  illustrated  in  Fig.  3  which 
shows  the  dependence  of  both  R  and  S  on  the  ap¬ 
plied  current  at  T =47  K  and  the  threshold  value 
is  the  same  in  each  case.  The  resistance  curves 
in  Figs.  2(a)  and  3  are  similar  to  the  data  of  Mon- 
ceau  et  al.,x  and  both  R  and  S  tend  to  saturate  to 
values  consistent  with  an  extrapolation  of  data 
above  the  CDW  transition  to  lower  temperatures. 

Attempts  were  also  made  to  observe  similar 
effects  at  the  higher-temperature  CDW  transition 
in  this  material,  but  due  to  the  higher  sample 
currents  required  to  see  the  breakdown  effects, 
sample  heating  was  a  greater  problem.  However, 
a  few  measurements  in  this  region  indicated  a  de¬ 
crease  in  TEP  toward  the  value  obtained  above 
145  K  although  no  accurate  data  were  obtained. 
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Several  models  have  been  proposed  to  explain 
the  non-Ohmic  conductivity  in  NbSe,.  In  the  orig¬ 
inal  experimental  papers1,5  the  nonlinear  charac¬ 
teristics  were  fitted  to  a  Zener  tunneling  model 
with  the  unphysically  small  gap  value  inferred 
from  the  breakdown  field.  The  possibility  that 
hot-electron  effects  may  be  responsible  for  such 
effects  has  been  discussed  for  TTF-TCNQ,1  but 
since  these  materials  have  undergone  a  charge - 
density-wave  transition,  several  groups1'*17  have 
suggested  that  the  electric  field  effects  are  re¬ 
lated  to  excitations  of  the  CDW.  If  the  incommen¬ 
surate  CDW  is  weakly  pinned,  an  electric  field 
should  be  able  to  depin  it,  and  a  contribution  to 
the  conductivity  would  result  from  the  sliding 
CDW  in  a  manner  Bimilar  to  that  proposed  by 
Frohlich.*  If  the  CDW  is  more  strongly  pinned, 
it  may  still  be  possible  to  observe  nonlinear  ef¬ 
fects  due  to  the  creation  of  low-energy  charged 
excitations  such  as  solitons.5 

The  TEP  rules  out  some  of  these  models.  In 
the  case  of  Zener  tunneling  or  hot-electron  ef¬ 
fects,  the  added  carriers  are  in  energy  states 
considerably  above  the  Fermi  energy.  They 
would  therefore  transport  a  good  deal  of  heat  and 
produce  a  large  contribution  to  the  TEP  similar 
to  that  observed  in  semiconductors.  The  magni¬ 
tude  of  the  TEP  increases  with  electric  field  only 
to  the  value  of  -  30  p  V/K  which  was  found  above 
the  CDW  transition  and  no  higher.  It  is  then  un¬ 
likely  that  the  Zener  tunneling  or  hot  electron 
models  are  appropriate. 

Recent  x-ray  measurements*  show  that  the  lat¬ 
tice  distortion  in  NbSe,  remains  unchanged  in  the 
presence  of  electric  fields  high  enough  to  show 
non-Ohmic  behavior,  which  implies  that  the  CDW 
remains  intact  under  such  coixiitions.  This 
means  that  models  involving  solitons  or  depinning 
charge  density  waves  can  be  discussed. 

In  a  soliton  model,  the  electric  field  creates 
charged  soliton-antisoliton  pairs  which  are  mo¬ 
bile  and  increase  the  conductivity.  While  the 
TEP  of  a  soliton  gas  has  not  yet  been  calculated, 
the  probable  result  would  be  similar  to  that  of  a 
semiconductor  whose  gap  is  the  sum  of  the  soli¬ 
ton-antisoliton  rest  mass  and  whose  mobilities 
are  given  by  the  soliton  dispersion  relations  and 
scattering  rates.  Since  these  calculations  have 
not  been  performed,  it  is  possible  that  a  soliton 
model  may  explain  our  results.  It  is,  however, 
unusual  that  in  an  electric  field  the  TEP  tends  to 
return  to  the  value  obtained  above  the  CDW  transi¬ 
tion.  This  would  imply  that  the  mobility  and  ener¬ 
gy  distribution  of  the  solitons  is  similar  to  that 
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of  the  electrons  above  the  transition,  or  that  some 
fortuitous  cancelation  has  occurred. 

If  the  additional  conductivity  in  the  non-Ohmic 
region  is  due  to  a  depinned  sliding  charge -density 
wave,  the  carriers  involved  are  in  a  condensate. 
This  condensed  state  of  electron-hole  pairs  re¬ 
sults  from  a  situation  in  which  a  finite  fraction 
of  conduction  electrons  occupy  their  ground- 
state  configuration.  The  condensate  has  no  en¬ 
tropy  and  can  transport  no  heat  and  therefore 
produces  no  TEP.  The  situation  is  similar  to 
that  which  arises  in  conventional  superconductors 
where  the  TEP  drops  to  zero  at  the  transition 
temperature,  the  supercurrent  arising  from  a 
condensate  with  no  entropy.  Therefore  a  sliding 
charge-density  wave  by  itsel f  will  not  contribute 
to  the  measured  TEP  which  means  that  some  oth¬ 
er  mechanism  must  exist  in  order  to  explain  the 
observed  contribution  to  the  TEP. 

The  most  direct  interpretation  of  the  TEP  data 
is  that  the  applied  field  is  freeing  the  single-par¬ 
ticle  carriers  which  have  been  frozen  at  the  CDW 
transition.  This  would  result  in  the  TEP  return¬ 
ing  to  a  value  consistent  with  the  pretransition 
value  (which  is  observed)  and  also  be  consistent 
with  the  observed  saturation  of  the  conductivity. 
However,  this  implies  that  the  CDW  is  being  de¬ 
pleted  by  the  electric  field  which  appears  incom¬ 
patible  with  the  interpretation  of  the  x-ray  data 
mentioned  earlier. 

hi  conclusion,  we  have  reported  the  first  ob¬ 
servation  of  a  nonlinear  or  electric -field -depen¬ 
dent  thermopower.  Our  results  suggest  that  slid¬ 
ing  charge-density  waves  alone,  Zener  tunneling, 
or  hot-electron  effects  are  not  responsible  for 


the  observed  thermopower.  The  simplest  inter¬ 
pretation  involves  the  electric -field  breakdown 
of  the  charge -density -wave  state.  In  the  light  of 
recent  x-ray  studies  on  NbSe}  it  is  clear  that 
further  work  is  required  to  resolve  this. 
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thus  the  addition  of  Ta  must  do  more  than  just 
squeeze  the  NbSej  lattice. 

In  conclusion,  we  have  seen  that  by  adding  5% 
Ta  to  NbSej  we  obtain  a  three  dimensional,  ani¬ 
sotropic  superconductor  with  a  Tc  *  1.55  i  0.2K. 
The  addition  of  Ta  suppresses  the  low  tempera¬ 


ture  CDW  by  slightly  lowering  Its  transition 
temperature  and  greatly  reducing  its  amplitude, 
thus  allowing  superconductivity. 
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the  coherence  lengths  from  figure  3.  One  finds 

that  C|((0)  -  230&  and  £j_(0)  -  98X.  This  leads 

to  an  anisotropy,  £~*  »  £|/£|  •  of  2.37. 

The  anisotropic  effect  mass  theory  predicts  that 

is  independent  of  temperature  as  long  as 

4(T)  is  larger  than  the  interchain  distance.  In 

NbSej  the  interchain  distance  is  ~  48  from  X-ray 

data.  3  TaSe-i  has  a  similar  interchain  dis- 
a  j 

tance°  therefore  one  would  expect  the  same  for 
Nb  g^Ta  gjSej.  Thus  the  condition  is  fulfilled 
for  all'T.  As  can  be  seen  in  Table  1,  £-l  is 


and  we  see  a  decrease  in  the  CDW  in  the  Ta 
doped  samples.  The  superconducting  transition 
is  sharp,  ~  103!  Tc  or  ~  303!  He 2,  which  is 
another  characteristic  of  three  dimensional 
superconductors.  Thus  the  difference  in  aniso¬ 
tropy  between  NbSej  and  Nb . 9jTa. QjSej  could  be 
due  to  the  fact  that  NbSej  is  quasi  one-dimen¬ 
sional  while  the  Ta  doped  material  is  three 
dimensional. 

In  pure  TaSej  Haen  et  al1  report  a  Tc  of 
less  than  1.5K.  Fleming  et  al'’  and  Sambongi  et 


Table  1  The  anisotropy  of  the  critical  fields  of  Nb  gjTa  a®  a  function 

of  temperature.  The  points  below  0.6K  are  from  the  extrapolated  data. 

H||  and  HJ_are  the  parallel  and  perpendicular  critical  fields  respectively. 


T 

Hll 

“J. 

e'1  -  H,/Hj_ 

1.2 

4.55 

1.88 

2.42 

1.0 

8.78 

3.69 

2.38 

0.8 

12.75 

5.35 

2.38 

0.6 

16.25 

6.82 

2.38 

0.4 

19.34 

8.05 

2.40 

0.2 

21.50 

9.06 

2.37 

0. 

•  22.30 

9.38 

2.38 

Independent  of  temperature  with  a  value  of  2.39 
t  0.02.  This  is  consistent  with  the  anisotropic 
effective  mass  theory.  Another  check  of  the 
applicability  of  this  model  is  to  check  the 
angular  dependence  of  the  critical  field.  This 
study  is  currently  underway. 

Nb_ gsTaQjSej  is  not  like  the  anisotropic 
superconductor  (SN)X  where  the  anisotropy  is  due 
to  (SN)x's  filamentary  nature.^  In 
Nk.95^a.05Se3  the  ^c2  versus  T  curve  does  not 
show  the  anomalous  concave  behaviour  near  T  =  Tc 
that  was  observed  in  (SN)X.^  The  HC2  versus  T 
curve  comes  in  linearly  to  the  axis  at  Tc  for 
both  parallel  and  perpendicular  fields  in 
Nb.gjTa.osSes  (figure  3).  It  has  been  seen  in 
(SN)X  that  the  perpendicular  critical  field 
loses  its  concave  behaviour  near  Tc  when  bromine 
is  added.  This  has  been  interpreted  as  being 
due  to  more  coupling  between  the  fibers. 15  The 
anisotropic  effective  mass  model  is  not  applica¬ 
ble  to  (SN)X  ^  since  £~1  is  not  independent  of 
temperature . 

The  addition  of  Ta  appears  to  cause  the 
NbSej  to  become  more  three  dimensional.  Ong  and 
Brillih  have  found  an  anisotropy  in  the  conduc¬ 
tivity  of  NbSej  of  20  t  5;  while  from 

our  critical  field  measurements  we  find  an 
anisotropy  of  2.39  t  0.02  for  Nb.gsTa.gjSej. 

This  shows  the  Ta  doped  samples  are  less  aniso¬ 
tropic  than  the  pure  NbSe3.  Another  sign  of  the 
Ta  making  NbSej  more  three  dimensional  is  the 
fact  that  in  three  dimensional  materials  CDW's 
do  not  have  as  pronounced  an  effect  on  the 
resistance  as  they  do  in  one  dimensional  ones. 


al*  report  a  Tc  of  2. IK  in  TaSej.  Thus,  the  Tc 
of  Nb_ 95Ta_o5Se3  is  of  the  same  order  as  that  of 
TaSej.  Sambongi  et  al®  also  say  that  at  1.25K 
or  0.6  Tc,  HC2  is  1.3  kOe  while  10  kOe  parallel 
to  the  chain  only  recovers  one  fourth  of  the 
normal  state.  This  leads  to  an  anisotropy  of 
greater  than  7.7.  Thus  It  appears  that 
Nb.pjTa.QjSej  is  less  anisotropic  than  TaSe3. 
From  looking  at  HC2  for  0.6  Tc  »  .92K  one  sees 
that  Hcj^in  Nh.pjTa  gjSej  is  three  times  the 
He 2 j_ in  TaSej,  whereas  Hc2||  is  smaller.  From 
this  one  can  say  that  we  are  not  seeing  strands 
of  TaSej  undergoing  the  superconducting  transi¬ 
tion. 

It  is  possible  that  the  Ta  in  Nb.gjTagjSej 
is  squeezing  the  NbSej  lattice.  From  pressure 
studies^  it  is  known  that  NbSej  goes  supercon¬ 
ducting  under  pressure;  the  CDW  are  also  sup¬ 
pressed  under  pressure.  Haen  et  al*®  find  that 
dTcuyj/dP  «-4K/kbar.  We  see  a  change  in  the  tem¬ 
perature  of  the  higher  temperature  CDW  of  7.5  i 
2K  in  Nb.95Ta>o5Se3  (figure  1).  This  would  cor¬ 
respond  to  a  pressure  of  1.9  ±  0.5  kbar.  Ac¬ 
cording  to  Monceau  et  al5  NbSej  goes  supercon¬ 
ducting  at  -  1.1K  under  a  pressure  of  -  2  kbar. 
This  is  in  rough  agreement  with  our  Tc  for 
N®.95^a.05®e3  1.5K.  Monceau  et  al5  measured 

the  magnetic  susceptibility  to  see  the  supercon¬ 
ducting  transition  whereas  we  have  measured  the 
transition  resistively,  and  this  may  account  for 
the  difference. 18  There  can  be  a  difference  in 
Tc  due  to  the  two  methods  of  measureing  Tc. 
However  a  pressure  of  2  kbar  is  not  enough  to 
completely  smear  out  the  lower  CDW  in  NbSej, 
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Figure  1  The  resistivity  normalized  to  the  room 
temperature  value  as  a  function  of  temperature 
for  representative  NbSe^  and  Nb  ^Ta  Q^Se^  samples. 

-  125  in  the  pure  NbSe3  to  1  -  1.7  in  the  Ta 
doped  samples.  It  is  expected  that  the  actual 
Ta  doping  is  5  ±  1%.  This  can  account  for  the 
differences  in  the  Ta  doped  samples.  The  ther¬ 
mopower  studies^  also  indicate  slight  varia¬ 
tions  in  Ta  concentration  in  samples  from  the 
nominal  5%  Ta  batch. 

The  resistivity  of  the  pure  and  doped  sam¬ 
ples  for  T  <  4.2K,  normalized  to  the  4.2K  val¬ 
ues,  is  shown  in  figure  2.  The  pure  NbSe^  dis- 


Figure  2  The  low  temperature  resistivity  normal¬ 
ized  to  the  4.2K  value  as  a  function  of  temper* 
ature  in  zero  magnetic  field.  The  samples  are  the 
same  as  those  in  figure  1. 


plays  no  transition  (to  0.5K)  just  a  low  temper¬ 
ature  tail  which  levels  off  at  about  3.5K.  The 
Ta  doped  samples  show  a  superconducting  transi¬ 
tion  with  Tc  =  1.55  ±  .2K,  where  Tc  is  the  mid¬ 
point  of  the  resistive  transition  and  the  error 
represents  10  to  90%  of  the  transition  as  well 
as  deviations  from  sample  to  sample.  Two  of  our 
samples  did  not  have  their  resistances  go  com¬ 
pletely  to  zero,  although  their  residual  resist¬ 
ance  was  less  than  10%  of  their  resistance  at 
4.2K.  This  is  probably  due  to  sample  inhomoge¬ 
neities.  There  was  no  observable  hysteresis  in 
the  transition  when  the  samples  were  cooled  and 
warmed  through  Tc. 

In  figure  3  we  have  plotted  the  upper  crit¬ 
ical  field,  Hc2,  versus  T  for  two  Ta  doped  sam¬ 
ples.  Hc2  was  taken  as  the  midpoint  in  the  re¬ 
sistive  transition.  The  fields  were  approxi¬ 
mately  (±  5°)  perpendicular  and  parallel  to  the 
chain  direction.  The  bars  represent  the  10  - 
90%  transition  width.  The  sharpness  of  the 
transition  can  be  seen  in  the  insert  in  figure  3. 

From  the  data  one  can  estimate  the  critical 
field  at  zero  temperature  Hc2(0).  One  finds 


Figure  3  The  upper  critical  field  parallel 

and  perpendicular,  as  a  function  of  temperature 
tor  Nb  ^Ta  ^^Se^*  The  dashed  lines  are  a  least 
squares  fit  of  the  data  to  the  Abrikosov-Gor'kov 
universal  function.  The  resistive  transitions  as 
a  function  of  magnetic  field  at  T  *  0.54K  are  shown 
in  the  insert. 

Hc2| (0)  -  22  kOe  while  Hc2  j.  -  9  kOe.  For  an 
anisotropic  superconductor,  the  Ginsburg-Landau 
(GL)  theory  of  superconductivity  gives 

dHlS  *o  .  dHi  ♦o  . 

- ”  - 5 —  and  -  -  -  ' 

dT  27t£/Tc  dT  2tt5||  SjTc 

■  7  2 

where  <(>„  *  2  x  10  Oe-cm  ,  £||  is  the  GL  coher¬ 
ence  length  along  the  chain  direction  and  is 
the  coherence  length  in  the  transverse  direc¬ 
tion.  Using  these  expressions  one  can  calculate 
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We  have  measured  the  resistivity  of  NbSej  doped  with  5%  Ta  from  room 
temperature  down  to  0. 5K  and  compared  our  results  with  similar  measure¬ 
ments  on  pure  NbSej.  The  pure  sample  remains  normal  to  the  lowest  tem¬ 
perature  (0.5K),  whereas  the  doped  sample  has  a  sharp  transition  to  the 
superconducting  state  with  Tc  *  1.5  ±  0.2K.  Measurements  of  the  criti¬ 
cal  magnetic  field  indicate  that  the  Ta  doped  samples  are  homogeneous, 
anisotropic  three  dimensional  superconductors. 


The  transition  metal  trichalcogenides, 

NbSe3  and  TaSe^,  are  highly  anisotropic  metals. 
NbSe3  has  two  anomalies  in  its  resistivity  as  a 
function  of  temperature,  one  at  145K  and  the 
other  at  59K.  The  anomalies  have  been  iden¬ 
tified  as  charge  density  wave  (CDW)  transi¬ 
tions.  Below  these  transitions,  NbSe^  shows 
striking  nonohmic  behaviour. 3»^  From  magnetiza¬ 
tion  measurements  Monceau  et  al^  do  not  see  any 
sign  of  superconductivity  down  to  50  mK  under 
atmospheric  pressure  in  NbSe3>  however  a  drop 
in  its  resistance  between  2.2K  and  1 . 5K  to  a 
plateau  and  then  a  linear  drop  below  0.4K  is 
seen  by  Haen  et  al.® 

In  TaSe^  there  are  no  resistive  anoma¬ 
lies.^*®  Two  groups,  Fleming  et  al^  and 
Sambongi  et  al,®  report  that  TaSe3  is  supercon¬ 
ducting  at  atmospheric  pressure  with  a  transi¬ 
tion  temperature,  Tc  of  ~  2,0  -  2. IK. 

It  has  previously  been  suggested^  that  CDW 
formation  and  superconductivity  are  antagonis¬ 
tic,  both  competing  for  states  near  the  Fermi 
surface.  The  NbSe3-TaSe3  system  is  therefore 
interesting  for  studying  this  competition. 
Pressure  measurements3  on  pure  NbSe3  have  shown 
an  increase  in  superconductivity  from  Tc  -  0  to 
Tc  ~  2.5K  with  a  pressure  of  7  kbar.  Pressure 
also  suppresses  the  CDW. 10  jt  both  lowers  the 
temperature  at  which  the  resistive  anomalies 
occur  and  diminishes  the  size  of  the  anomalies. 
Thus  it  appears,  from  the  measurements  on  NbSe3» 
that  by  suppressing  the  CDW  one  is  able  to 
obtain  superconductivity. 

In  the  present  study  we  have  used  the 
alloying  of  NbSe3  with  TaSe3  to  suppress  the  CDW 
transition  tn  a  quite  separate  way.  The  addi¬ 


tion  of  5%  TaSe3  to  NbSe3  is  sufficient  to  smear 
the  CDW  transition  and  produce  a  Tc  comparable 
to  that  of  TaSej. 

We  have  measured  the  resistivity,  p,  versus 
temperature,  T,  for  several  samples  of 
Nb ,95Ta.05®e3  anc*  The  samples  were  rib¬ 

bon  like  with  typical  dimensions  being  4.5  mm 
long,  0.01  -  0.02  mm  wide  for  those  samples  with 
Ta,  0.02  -  0.04  mm  wide  for  NbSe3»  and  '  0.006 
mm  thick.  The  resistance  was  measured  along  the 
chain  direction  using  the  usual  four  probe  meth¬ 
od  with  1  mil  gold  leads  silver  painted  to  the 
sample.  The  low  temperatures  were  achieved  in  a 
He3  cryostat,  the  samples  being  either  in  gaseous 
He3  or  submerged  in  liquid  He^  depending  on  the 
temperature.  Both  DC  and  low  frequency  (<  100 
Hz)  AC  measurements  were  taken  with  a  current 
through  the  samples  of  1  pA.  The  AC  measure¬ 
ments  were  performed  with  a  lock-in  voltmeter. 
Critical  field  measurements  were  made  on  the 
superconducting  samples  by  monitoring  the  resis¬ 
tance  while  sweeping  a  magnetic  field  at  con¬ 
stant  temperature  (±  2  mk) . 

The  resistivity,  normalized  to  its  room 
temperature  value,  is  shown  in  figure  1  for  both 
the  pure  NbSe3  and  Ta  doped  samples.  Here  one 
can  see  that  the  CDW's  smear  out  with  the  addi¬ 
tion  of  the  Ta  impurity.  The  lower  temperature 
CDW  anomaly  is  almost  nonexistent.  Ong  et  al** 
have  done  a  much  more  detailed  study  of  this 
transition.  In  thermopower  measurements  Kwak  et 

have  also  seen  that  the  145K  transition  is 
less  affected  than  the  59K  transition  for  the  5% 
doping.  Those  samples  with  Ta  are  much  more 
disordered  than  the  pure  NbSe^  samples.  The  re¬ 
sistivity  ratio,  p(300)/p(4.2) ,  changes  from 
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